BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	732,445	A→T	9.3%	Y661F (TAC→TTC) ‡	clpB →	chaperone protein ClpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	732,445	0	A	T	9.3%	741.7 / 49.2	236	Y661F (TAC→TTC) ‡	clpB	chaperone protein ClpB
Reads supporting (aligned to +/- strand): ref base A (106/108); new base T (9/13); total (115/121)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 239 total aligned reads displayed.

GACCGAACTGTGCAAGGCGCTGGCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCA  >  NC_002947/732296‑732594

gACCGAACTGTGCAAGGCGCTGGCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTa                                                                                                                                                        >  2:4354373/1‑150 (MQ=255)

gACCGAACTGTGCAAGGCGCTGGCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTa                                                                                                                                                        >  2:6805359/1‑150 (MQ=255)

 aCCGAACTGTGCAAGGCGCTGGCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGTGTATGAAGAGGGCGGTTAc                                                                                                                                                       <  1:6189296/150‑1 (MQ=255)

   cGAACTGTGCAAGGCGCTGGCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCt                                                                                                                                                     >  2:6263853/1‑150 (MQ=255)

         gtgCAAGGCGCTGGCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGa                                                                                                                                               <  2:1747115/150‑1 (MQ=255)

                 cgcTGGCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTgc                                                                                                                                       <  1:415703/150‑1 (MQ=255)

                  gcTGGCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTgcg                                                                                                                                      >  1:6724865/1‑150 (MQ=255)

                    tGGCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTgcggc                                                                                                                                    <  1:3526019/150‑1 (MQ=255)

                     ggCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTgcggcg                                                                                                                                   >  1:2432918/1‑150 (MQ=255)

                      gCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGgcgc                                                                                                                                  >  2:7518381/1‑150 (MQ=255)

                        cGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTCTGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCaa                                                                                                                                >  2:5621505/1‑150 (MQ=255)

                        cGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCaa                                                                                                                                >  2:5825801/1‑150 (MQ=255)

                          aGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGc                                                                                                                              <  2:1836003/150‑1 (MQ=255)

                              ccTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTa                                                                                                                          <  2:783763/150‑1 (MQ=255)

                                    cGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCggt                                                                                                                    <  1:4046656/150‑1 (MQ=255)

                                     gACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCggtg                                                                                                                   <  1:4659799/150‑1 (MQ=255)

                                     gACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCggtg                                                                                                                   <  2:2170430/150‑1 (MQ=255)

                                        acCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGgtggtg                                                                                                                >  2:5527061/1‑150 (MQ=255)

                                        acCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGgtggtg                                                                                                                >  2:1668960/1‑150 (MQ=255)

                                           gAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGtgctg                                                                                                             >  1:1441427/1‑150 (MQ=255)

                                             agagGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTgctgct                                                                                                           <  2:2264429/150‑1 (MQ=255)

                                                     tGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAgg                                                                                                   <  2:8894103/150‑1 (MQ=255)

                                                     tGGTGCGCAACGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAgg                                                                                                   <  2:8347083/150‑1 (MQ=255)

                                                       gTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTg                                                                                                 <  2:5317150/150‑1 (MQ=255)

                                                       gTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTg                                                                                                 <  2:8830477/150‑1 (MQ=255)

                                                        tGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTgg                                                                                                >  2:7629523/1‑150 (MQ=255)

                                                           acATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGaaa                                                                                             <  2:8637448/149‑1 (MQ=255)

                                                             aTCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGaaaaa                                                                                           <  2:8154561/150‑1 (MQ=255)

                                                              tCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAg                                                                                          <  1:7629523/150‑1 (MQ=255)

                                                               cGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGc                                                                                         <  1:1909250/150‑1 (MQ=255)

                                                                 aCATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGccc                                                                                       >  2:8596529/1‑150 (MQ=255)

                                                                 aCATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGccc                                                                                       <  1:6425320/150‑1 (MQ=255)

                                                                   aTGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCAc                                                                                     >  2:2279442/1‑150 (MQ=255)

                                                                     gTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCAccc                                                                                   <  1:4102836/150‑1 (MQ=255)

                                                                      tCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCg                                                                                  >  2:6715528/1‑150 (MQ=255)

                                                                      tCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCg                                                                                  >  2:1409313/1‑150 (MQ=255)

                                                                      tCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCg                                                                                  >  1:357332/1‑150 (MQ=255)

                                                                         gAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGAt                                                                               <  2:4014887/150‑1 (MQ=255)

                                                                         gAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGAt                                                                               <  1:3362127/150‑1 (MQ=255)

                                                                         gAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGAt                                                                               <  1:1076046/150‑1 (MQ=255)

                                                                             tCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATgtgt                                                                           >  2:2233501/1‑150 (MQ=255)

                                                                              cATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGtt                                                                          <  2:6061784/150‑1 (MQ=255)

                                                                               aTGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTc                                                                                                                       <  1:361252/104‑1 (MQ=255)

                                                                               aTGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTc                                                                                                                       >  2:361252/1‑104 (MQ=255)

                                                                               aTGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTc                                                                         >  1:1742086/1‑150 (MQ=255)

                                                                                tGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCa                                                                        >  2:6434251/1‑150 (MQ=255)

                                                                                  gagaAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAAc                                                                      >  2:495342/1‑150 (MQ=255)

                                                                                  gagaAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAAc                                                                      <  1:6470458/150‑1 (MQ=255)

                                                                                   agaAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACg                                                                     >  1:2301782/1‑150 (MQ=255)

                                                                                     aaGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTg                                                                   >  1:2614323/1‑150 (MQ=255)

                                                                                       gCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCt                                                                 >  1:7049506/1‑150 (MQ=255)

                                                                                       gCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCt                                                                 >  2:3006119/1‑150 (MQ=255)

                                                                                             cGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCa                                                           >  1:8805439/1‑150 (MQ=255)

                                                                                             cGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCa                                                           <  1:5157677/150‑1 (MQ=255)

                                                                                             cGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCa                                                           <  2:7191334/150‑1 (MQ=255)

                                                                                              gTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAg                                                          <  1:2642052/150‑1 (MQ=255)

                                                                                                gagccgtcTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGt                                                        <  1:9162723/143‑1 (MQ=255)

                                                                                                 gcgcgcCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTg                                                       <  2:3750873/150‑1 (MQ=255)

                                                                                                   gcgcCTGATCGGTGCACCACCAGGCAATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCATACACGGTGGTGCTGCTGGACGAGGTGGAAAAAGACCACCCCGATGTGTTAAAAGTGCTGTTTCATGTGCt                                                     >  1:8596929/1‑150 (MQ=255)

                                                                                                      tcTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGa                                                  >  2:3416529/2‑150 (MQ=255)

                                                                                                      ccTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAATCCTTACTCGGTGGTGCTGCTGGACGAGGTGTAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGa                                                  <  1:1288983/150‑1 (MQ=255)

                                                                                                      ccTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGa                                                  <  1:5388745/150‑1 (MQ=255)

                                                                                                       cTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGaa                                                 <  1:2369024/150‑1 (MQ=255)

                                                                                                        tGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAg                                                >  2:386966/1‑150 (MQ=255)

                                                                                                          aTCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGAc                                              <  1:5560604/150‑1 (MQ=255)

                                                                                                           ttggTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACg                                             <  1:7392211/148‑1 (MQ=255)

                                                                                                            cGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACgg                                            <  2:5338503/150‑1 (MQ=255)

                                                                                                             ggagCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTAAACGTGCTGTTGCAGGTGCTGGAAGACGgc                                           >  1:2783514/4‑150 (MQ=255)

                                                                                                              gTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgcc                                          >  2:4048676/1‑150 (MQ=255)

                                                                                                                  accaccAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGc                                                      <  1:3640302/134‑1 (MQ=255)

                                                                                                                  accaccAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCt                                      >  2:182672/1‑150 (MQ=255)

                                                                                                                     tcccggttATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGAc                                   >  1:5040665/8‑150 (MQ=255)

                                                                                                                     tcccggttATGTCGGCTATGAAGAGGGTGGCGTTCTGAAGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGAc                                   >  2:7455129/8‑150 (MQ=255)

                                                                                                                     accaGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGAc                                   >  2:2872114/1‑150 (MQ=255)

                                                                                                                     accaGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGAc                                   <  1:2094474/150‑1 (MQ=255)

                                                                                                                       caGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCg                                 >  1:8079580/1‑150 (MQ=255)

                                                                                                                        aGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGa                                <  1:8373731/150‑1 (MQ=255)

                                                                                                                           cTATGTAGGGTATGAAGAGGGCGGTGACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAg                             <  1:8740818/150‑1 (MQ=255)

                                                                                                                                 aGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCTGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACgg                       <  2:342868/150‑1 (MQ=255)

                                                                                                                                  gggTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACgg‑g                     >  2:4413221/1‑150 (MQ=255)

                                                                                                                                  ggctATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGAc                                   <  1:7825273/134‑1 (MQ=255)

                                                                                                                                    gTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑gcg                   >  2:4725646/1‑150 (MQ=255)

                                                                                                                                      aTGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCa                 <  1:4354373/150‑1 (MQ=255)

                                                                                                                                      aTGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCa                 <  1:6805359/150‑1 (MQ=255)

                                                                                                                                      aTGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCa                 >  1:8534463/1‑150 (MQ=255)

                                                                                                                                       tGAAGAGGTCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCAc                >  1:8197970/1‑150 (MQ=255)

                                                                                                                                       tGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCAc                >  1:1728307/1‑150 (MQ=255)

                                                                                                                                        gAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCAcc               >  1:3048433/1‑150 (MQ=255)

                                                                                                                                         aaGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCCCCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCg              >  2:2962864/1‑150 (MQ=255)

                                                                                                                                         aaGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCg              >  2:8734498/1‑150 (MQ=255)

                                                                                                                                         aaGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCg              >  2:4671342/1‑150 (MQ=255)

                                                                                                                                           gagGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTg            >  1:4707677/1‑150 (MQ=255)

                                                                                                                                           gagGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTg            >  1:3684047/1‑150 (MQ=255)

                                                                                                                                           gagGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTg            <  2:1122795/150‑1 (MQ=255)

                                                                                                                                            agGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGG‑AAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGTGCGC‑CCGTggc          <  1:1520511/150‑2 (MQ=255)

                                                                                                                                             gggCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGa          <  2:2679872/150‑1 (MQ=255)

                                                                                                                                              ggCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGTGCGC‑CCGTggct         <  2:5647422/150‑3 (MQ=255)

                                                                                                                                                  gTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTcc     >  2:2480633/1‑150 (MQ=255)

                                                                                                                                                  gTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTcc     <  1:5915628/150‑1 (MQ=255)

                                                                                                                                                    tACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGc   <  1:4420156/150‑1 (MQ=255)

GACCGAACTGTGCAAGGCGCTGGCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCA  >  NC_002947/732296‑732594

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