BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,501	G→A	6.5%	A17A (GCC→GCT)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,501	0	G	A	6.5%	390.3 / 9.8	165	A17A (GCC→GCT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base G (66/86); new base A (4/7); total (70/95)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 167 total aligned reads displayed.

CTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATA  >  CP000730/36363‑36641

cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                                            >  2:804942/1‑141 (MQ=35)

 tttttttttttATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTTTTTTTTTTTTTTTATCACGATAACTCAATGTATATCTTATATAGTAGCGAACGGCt                                                                                                                                           >  2:187431/8‑141 (MQ=11)

  tcCAAATTTGATGCAAAATTGGGGCATATTCTTGAACCCAACGGTACACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGACATTAGCCAACGGCTa                                                                                                                                          >  1:642584/1‑141 (MQ=255)

   ccAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTAc                                                                                                                                         >  1:799020/1‑141 (MQ=35)

        tttGATACAAAATTGGGGTATAGGCTGGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGa                                                                                                                                    <  2:269987/141‑1 (MQ=255)

        tttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGa                                                                                                                                    <  2:316745/141‑1 (MQ=34)

         ttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                                                   <  1:628033/141‑1 (MQ=35)

           gATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATaa                                                                                                                                 >  1:76439/1‑141 (MQ=35)

             tACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACa                                                                                                                               <  2:159597/141‑1 (MQ=35)

                 aaaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCt                                                                                                                           <  2:409808/141‑1 (MQ=32)

                    ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                                                        <  1:84229/141‑1 (MQ=32)

                    tcggtgcctatTCTTGAACCCAACTGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                                                        <  2:635442/133‑1 (MQ=12)

                         gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                                                   <  2:746652/141‑1 (MQ=32)

                         gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                                                   >  1:799230/1‑141 (MQ=32)

                         gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAAAGGCTACAGTGATAACACCCTTGTTAAAt                                                                                                                   >  2:258617/1‑141 (MQ=17)

                                   gAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATc                                                                                                         >  1:635475/1‑141 (MQ=32)

                                     aCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCTATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTg                                                                                                       <  2:11494/141‑1 (MQ=17)

                                      cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                                                      <  2:766076/141‑1 (MQ=32)

                                        cAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                                                    >  2:483561/1‑141 (MQ=32)

                                          aCGGTAGACCGTTTACTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATa                                                                                                  <  2:140983/141‑1 (MQ=17)

                                            ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                                                >  1:326928/1‑141 (MQ=32)

                                               ctaCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCa                                                                                             >  2:426471/3‑141 (MQ=17)

                                               aGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCa                                                                                             <  2:835493/141‑1 (MQ=32)

                                                 aCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATACTTGTTAAATTGTTTATATCTGAAATAGTTCATa                                                                                           >  1:236753/1‑141 (MQ=17)

                                                        aaTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaaga                                                                                    >  1:603272/1‑141 (MQ=32)

                                                        aaTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaaga                                                                                    >  1:839384/1‑141 (MQ=32)

                                                          tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACt                                                                                  >  2:792175/1‑141 (MQ=32)

                                                             agaACGTTTACACCACGTGCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCt                                                                               <  2:128590/140‑1 (MQ=17)

                                                              gaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtt                                                                              >  2:198487/1‑141 (MQ=32)

                                                                aCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtttt                                                                            >  2:165490/1‑141 (MQ=32)

                                                                 cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttttt                                                                           >  1:904890/1‑141 (MQ=32)

                                                                   ttgacacCACGTTCCCTTAATAGTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGt                                                                         <  2:171591/138‑1 (MQ=17)

                                                                      acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                                                      <  1:763666/141‑1 (MQ=32)

                                                                            cGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATtat                                                                <  1:1026/141‑1 (MQ=32)

                                                                                    tAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaaa                                                        <  1:827325/141‑1 (MQ=34)

                                                                                        atTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                                                    >  2:180952/1‑141 (MQ=35)

                                                                                              gATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTg                                              >  1:845402/1‑141 (MQ=35)

                                                                                              gATATATCACGATAACTCAATGTATATATTATATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTTAAATTATACTATAAATTCAACtttt                                              >  2:838110/1‑140 (MQ=255)

                                                                                               atatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGc                                             <  1:603254/141‑1 (MQ=35)

                                                                                                 atatCACGATAACTCAATGTTTATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATCAATTCAACTTTGCaa                                           <  1:673840/141‑1 (MQ=21)

                                                                                                 atatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCaa                                           <  1:14002/141‑1 (MQ=35)

                                                                                                 atatCACGATAACTCAATGTATATCTTAGATAGAAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCaa                                           >  1:128668/1‑141 (MQ=25)

                                                                                                  tatCACGATAACTCCATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAc                                          <  1:700923/141‑1 (MQ=25)

                                                                                                  tatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCaaa                                          >  1:613374/1‑140 (MQ=25)

                                                                                                  tatCACGATAACTCAATGTATATCTTAGATAGGAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAc                                          <  1:246474/141‑1 (MQ=25)

                                                                                                  taccACGATAACTCAATGCATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAc                                          <  2:358165/138‑1 (MQ=0)

                                                                                                   atCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACa                                         <  1:316275/141‑1 (MQ=35)

                                                                                                    tCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATCCAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAg                                        <  1:820189/141‑1 (MQ=255)

                                                                                                    tCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTTTTATTATACTATAAATTCAACTTTGCAACAg                                        >  2:628344/1‑141 (MQ=255)

                                                                                                     cACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGa                                       >  1:291372/1‑141 (MQ=35)

                                                                                                      aCGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGaa                                      >  1:622792/1‑141 (MQ=35)

                                                                                                         aTAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCg                                   <  2:72555/141‑1 (MQ=9)

                                                                                                          tatcTCAACGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTTAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTACATTATAAATTCAACTTTGCAACAGAACCGt                                  <  1:42087/138‑1 (MQ=11)

                                                                                                          tagttCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                                  <  2:152066/137‑1 (MQ=9)

                                                                                                          tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                                  <  1:191032/138‑1 (MQ=12)

                                                                                                          tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGCGAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                                  >  2:878546/4‑141 (MQ=11)

                                                                                                          tAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt                                  >  1:675169/1‑141 (MQ=37)

                                                                                                          tAACTCAATGAATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTTTAATTCTTCTATAAATTCACCTTTTCACCAGaccgtc                                  >  2:421477/1‑136 (MQ=11)

                                                                                                           agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAAATTTGCAACAGAACCGta                                 >  2:507103/3‑141 (MQ=9)

                                                                                                           agcTCAAAGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta                                 <  1:33037/139‑1 (MQ=12)

                                                                                                           agatCAATGCATTGCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta                                 <  1:162961/138‑1 (MQ=11)

                                                                                                            gcttaaTGCATATCTTAGCTAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat                                <  1:148706/137‑1 (MQ=255)

                                                                                                            gcttaaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat                                <  2:869792/137‑1 (MQ=255)

                                                                                                            gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat                                >  1:357550/2‑141 (MQ=17)

                                                                                                            aCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat                                >  1:225804/1‑141 (MQ=39)

                                                                                                             cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt                               <  1:193331/141‑1 (MQ=21)

                                                                                                              tCAATGCATATCTTAGGTAGTGGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt                              >  1:465728/1‑140 (MQ=11)

                                                                                                              tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt                              >  1:496056/1‑140 (MQ=16)

                                                                                                              tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt                              <  1:634326/141‑2 (MQ=19)

                                                                                                              tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt                              >  2:326846/1‑140 (MQ=12)

                                                                                                               cAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtattat                             <  2:177978/141‑1 (MQ=40)

                                                                                                               cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt                             <  2:37666/141‑3 (MQ=255)

                                                                                                                  ttcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatg                          <  2:395794/138‑6 (MQ=255)

                                                                                                                  agcattttttagttagTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCCTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatg                          <  1:611092/133‑6 (MQ=11)

                                                                                                                   tcatatCTTAGATAGTAGCCAACGGCTACAGTGATATCCTCTTTGGTATGCTTTTTATATCTGAAATAGTTCATACAGAAGATTCCTTTTTGTTAAAATTATATTATCACTTCAACTTTGCACCAGAACCGtattttacga                         <  2:874800/139‑7 (MQ=11)

                                                                                                                   gcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgta                         <  2:332651/139‑12 (MQ=11)

                                                                                                                   gcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgta                         <  2:560298/139‑12 (MQ=11)

                                                                                                                   gcatatCATAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgta                         <  2:736995/139‑12 (MQ=11)

                                                                                                                    catatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtat                        >  2:60833/2‑129 (MQ=11)

                                                                                                                    catatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGAATCCTTTTTGTTAAAGTTGTATTATAAATTCAAATTTGCAACAGAAcctcctttaatgtt                        >  2:546936/2‑129 (MQ=11)

                                                                                                                     atatGTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTATATAACTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAAGTCAACTTTGCAACAGAAcatcctaaccgtatt                       >  2:796574/1‑127 (MQ=11)

                                                                                                                     atatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtatt                       >  1:572117/1‑128 (MQ=11)

                                                                                                                     atatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAAAGAACCGtattttatgagt                       >  2:269609/1‑133 (MQ=11)

                                                                                                                      tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagta                      >  2:97809/1‑132 (MQ=11)

                                                                                                                      tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtatta                      <  1:392040/141‑15 (MQ=11)

                                                                                                                      tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagta                      <  2:233923/141‑10 (MQ=255)

                                                                                                                       atCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATaga                     >  1:849932/1‑141 (MQ=255)

                                                                                                                         cgtAGATAGTAGCCAACGGCTACAGTGATAAAATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAgaga                   <  1:414035/139‑1 (MQ=255)

                                                                                                                         cTTAGATAGTAACCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATTCACACCATGGAACTTTGCAACAGAACCGTATTATGGAATAgaga                   <  1:436656/141‑1 (MQ=14)

                                                                                                                          ttAGATAGGAGCCAACGGCTACAGTGATAACATCCTTGTTAAATGGTTTATATATGGTGGAGCCCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacact                  <  1:432210/141‑14 (MQ=11)

                                                                                                                           tAGATAGTAGCCCACCGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc                 <  1:304388/141‑15 (MQ=11)

                                                                                                                              aTAGTAGTCAACTGCAACAGTCATAACATCCGTGTTAAATTGTTTATATCTGAAATAGTTCATACAGCAGACTCCCTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGCTGTTg              <  2:588108/141‑1 (MQ=17)

                                                                                                                                      cAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGATCCGTAGTATGGACTAGAGATGTTGGTAACAtt      <  2:602391/141‑1 (MQ=255)

                                                                                                                                        aCGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTta    <  2:90200/141‑1 (MQ=255)

                                                                                                                                        aCGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTta    <  1:781071/141‑1 (MQ=255)

                                                                                                                                        aCGGCTAAAGGGAAAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTAATACGGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCCACAGACCCGTATTATGGAAGAGAGATGTTGGTacctctac    <  2:236753/141‑8 (MQ=17)

                                                                                                                                         cGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAAACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCACCTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTtat   <  1:496599/141‑1 (MQ=255)

                                                                                                                                         cGGCTACAGTGATAACATCCTTGTTAAATTGTTTATAACTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTtat   <  1:337642/141‑1 (MQ=255)

                                                                                                                                          ggCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTtata  >  1:147273/1‑141 (MQ=255)

                                                                                                                                          ggCTACAGTGATAACATCCTTGTAAAATTGCTTATATCTGAAATAGTTCATCCAGAAGACTACTTTTTGTTAAAATTATACTATATATTCAACTTTGCAAAAGAACCGTATTATGGAATAGAGATGTTGATAACATTtcta  >  2:136814/1‑141 (MQ=17)

CTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATA  >  CP000730/36363‑36641

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`