BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	732,534	T→C	10.8%	L691L (TTG→CTG)	clpB →	chaperone protein ClpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	732,534	0	T	C	10.8%	517.7 / 39.2	184	L691L (TTG→CTG)	clpB	chaperone protein ClpB
Reads supporting (aligned to +/- strand): ref base T (88/76); new base C (5/15); total (93/91)
Fisher's exact test for biased strand distribution p-value = 1.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 187 total aligned reads displayed.

CTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATG  >  NC_002947/732386‑732683

 tCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTg                                                                                                                                                      >  2:2180639/1‑150 (MQ=255)

 tCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGAATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTg                                                                                                                                                      >  1:4994429/1‑150 (MQ=255)

      cgcgcgCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGt                                                                                                                                                 <  2:2694433/150‑1 (MQ=255)

          cgtcTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCtt                                                                                                                                             <  2:1006418/147‑2 (MQ=255)

           gcCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTgg                                                                                                                                            >  2:4397567/1‑150 (MQ=255)

            tcTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGa                                                                                                                                           <  2:2857147/149‑1 (MQ=255)

             cTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGat                                                                                                                                          <  1:2354814/150‑2 (MQ=255)

             cTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCGGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGat                                                                                                                                          <  1:318518/150‑2 (MQ=255)

             cTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGaa                                                                                                                                          >  2:3059935/1‑150 (MQ=255)

             cTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGaa                                                                                                                                          >  2:504895/1‑150 (MQ=255)

               gATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGa                                                                                                                                        >  1:3025363/1‑150 (MQ=255)

                 tCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACg                                                                                                                                      <  2:4856873/150‑1 (MQ=255)

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

                   ggTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgc                                                                                                                                    >  2:148460/1‑150 (MQ=255)

                    gTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgcc                                                                                                                                   <  1:5041745/150‑1 (MQ=255)

                      gCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgccgc                                                                                                                                 <  2:4337189/150‑1 (MQ=255)

                       caccacCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgccgcc                                                                                                                                >  2:708257/1‑150 (MQ=255)

                          caccaGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGa                                                                                                                             >  1:3386690/1‑150 (MQ=255)

                          caccaGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGa                                                                                                                             >  2:4945908/1‑150 (MQ=255)

                          caccaGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGa                                                                                                                             >  2:1567972/1‑150 (MQ=255)

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

                             ccggttATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCg                                                                                                                          >  1:3994997/6‑150 (MQ=255)

                                 cTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAg                                                                                                                      <  2:3839826/150‑1 (MQ=255)

                                  tATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCAGATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGAcggt                                                                                                                     <  2:4564983/150‑4 (MQ=255)

                                    tGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGAcggtca                                                                                                                   <  1:3067712/150‑6 (MQ=255)

                                       aGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACgg                                                                                                                <  2:2527804/150‑1 (MQ=255)

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

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

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

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

                                                 gagGGCTGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTg                                                                                                     <  1:1720666/150‑1 (MQ=255)

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

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

                                                    gtcgtttACCTGACCGAGGCAGTGCGTCGCAAGCCTTTCTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACTGCCGCCTGACCGACAGCCTCGGTGCGC‑CCGTggct                                                                                                  <  2:3419683/145‑3 (MQ=37)

                                                        gcgttcTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGG‑CCGTACGGTCGATTTcc                                                                                              <  2:1928968/145‑1 (MQ=255)

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

                                                         cgttcTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGG‑CCGTACGGTCGATTTCCg                                                                                             <  1:4710728/146‑1 (MQ=255)

                                                           aCCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCa                                                                                           >  2:2079503/1‑150 (MQ=255)

                                                           aCCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCa                                                                                           >  2:2079224/1‑150 (MQ=255)

                                                            ccTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCaa                                                                                          >  1:4714913/1‑150 (MQ=255)

                                                             cTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAac                                                                                         >  1:3920785/1‑150 (MQ=255)

                                                             cTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAac                                                                                         >  2:1796831/1‑150 (MQ=255)

                                                             cTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAac                                                                                         >  1:259370/1‑150 (MQ=255)

                                                               gACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAacac                                                                                       >  1:3154591/1‑150 (MQ=255)

                                                                 ccGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACcg                                                                                     <  2:1463012/150‑1 (MQ=255)

                                                                 ccGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACcg                                                                                     >  1:823569/1‑150 (MQ=255)

                                                                  cGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACcgt                                                                                    <  2:1343606/150‑1 (MQ=255)

                                                                   gAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACcgtg                                                                                   >  2:1247745/1‑150 (MQ=255)

                                                                    aGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACcgtga                                                                                  >  1:1129830/1‑150 (MQ=255)

                                                                              gcgcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGAc                                                                        <  1:285762/150‑1 (MQ=255)

                                                                              gcgcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGAc                                                                        >  2:823378/1‑150 (MQ=255)

                                                                                gcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCt                                                                      >  1:2698313/1‑150 (MQ=255)

                                                                                gcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCt                                                                      >  1:2698472/1‑150 (MQ=255)

                                                                                  aaGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTcc                                                                    >  2:2274709/1‑150 (MQ=255)

                                                                                    gCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCaa                                                                  >  1:2271628/1‑150 (MQ=255)

                                                                                     ccTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAAc                                                                 >  1:254600/1‑150 (MQ=255)

                                                                                          cTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTggg                                                            <  1:3777210/150‑1 (MQ=255)

                                                                                            aggtgatCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGG‑CCGTACGGTCGATTTCCGCAACACGGTCATCATCATGACCTCGAACCTCGGcg                                                          >  2:866040/2‑149 (MQ=255)

                                                                                             ggtggtGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTc                                                         <  2:4936704/150‑1 (MQ=255)

                                                                                                  tgctgcTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGgcgc                                                    <  1:571949/150‑1 (MQ=255)

                                                                                                   gctgctGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCa                                                   <  2:2971601/150‑1 (MQ=255)

                                                                                                       ctgGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGAtt                                               >  1:4911611/1‑150 (MQ=255)

                                                                                                        tgGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTc                                              <  1:4399040/150‑1 (MQ=255)

                                                                                                          gACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAg                                            >  2:3974035/1‑150 (MQ=255)

                                                                                                          gACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAg                                            <  1:139721/150‑1 (MQ=255)

                                                                                                          gACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAg                                            <  1:1963749/150‑1 (MQ=255)

                                                                                                           aCGAGGTGGAAAAAGCCCACCCGGATGTGTTCACCGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAgg                                           <  2:3393983/150‑1 (MQ=255)

                                                                                                              aGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAAc                                        >  1:4885391/1‑150 (MQ=255)

                                                                                                                  ggAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACtggt                                    >  1:3723069/1‑150 (MQ=255)

                                                                                                                  ggAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACtggt                                    <  2:2265685/150‑1 (MQ=255)

                                                                                                                  ggAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACtggt                                    >  1:1605247/1‑150 (MQ=255)

                                                                                                                   gAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACtggtg                                   <  1:148460/150‑1 (MQ=255)

                                                                                                                    aaaaaGCCCTCCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACtggtgg                                  <  1:1554488/150‑1 (MQ=255)

                                                                                                                    aaaaaGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACtggtgg                                  <  1:2506504/150‑1 (MQ=255)

                                                                                                                     aaaaGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTggtggg                                 >  2:1690901/1‑150 (MQ=255)

                                                                                                                      aaaGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTggtgggt                                <  2:5057831/150‑1 (MQ=255)

                                                                                                                      aaaGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTggtgggt                                <  1:2535266/150‑1 (MQ=255)

                                                                                                                       aaGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTggtgggtg                               >  1:1490662/1‑150 (MQ=255)

                                                                                                                           ccACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACcg                           >  2:2199853/1‑150 (MQ=255)

                                                                                                                            cACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACcgc                          >  2:2362847/1‑150 (MQ=255)

                                                                                                                              cccGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGa                        <  1:473241/150‑1 (MQ=255)

                                                                                                                              cccGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGa                        <  1:923043/150‑1 (MQ=255)

                                                                                                                              cccGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGa                        <  1:3625296/150‑1 (MQ=255)

                                                                                                                               ccGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAg                       <  1:116913/150‑1 (MQ=255)

                                                                                                                               ccGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAg                       >  2:4260044/1‑150 (MQ=255)

                                                                                                                                cGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAAAACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAgg                      <  1:3729258/150‑1 (MQ=255)

                                                                                                                                 ggATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGc                     >  1:3639696/1‑150 (MQ=255)

                                                                                                                                       gtTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCg               >  2:4712697/1‑150 (MQ=255)

                                                                                                                                         tCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTg             <  2:4994429/150‑1 (MQ=255)

                                                                                                                                         tCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTg             <  2:2424574/150‑1 (MQ=255)

                                                                                                                                          cAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGc            <  2:4327562/150‑1 (MQ=255)

                                                                                                                                           aaCGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCg           >  1:3120736/1‑150 (MQ=255)

                                                                                                                                           aaCGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCg           >  2:1534186/1‑150 (MQ=255)

                                                                                                                                              gTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCa        <  1:2180639/150‑1 (MQ=255)

                                                                                                                                              gTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCa        <  1:4394061/150‑1 (MQ=255)

                                                                                                                                                  tGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTga    >  1:77020/1‑150 (MQ=255)

                                                                                                                                                  tGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGAGTGACCGCGAGGCACAGCGTGCGGCAGTga    >  1:4868372/1‑150 (MQ=255)

                                                                                                                                                   gTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTgat   <  2:3091231/150‑1 (MQ=255)

                                                                                                                                                   gTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTgat   <  2:1511557/150‑1 (MQ=255)

                                                                                                                                                    ttGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTgatg  <  2:2520963/150‑1 (MQ=255)

CTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATG  >  NC_002947/732386‑732683

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: `‑`