Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,838,555	G→C	5.4%	G703R (GGC→CGC)	PP_3388 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,838,555	0	G	C	5.4%	281.8 / 6.9	92	G703R (GGC→CGC)	PP_3388	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (36/51);  new base C (2/3);  total (38/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

TGGAGGGGATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTAATATACCGCCAGCAAAA  >  NC_002947/3838406‑3838703                                                                                                                                                        |                                                                                                                                                     tgGAGGGGATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACg                                                                                                                                                      >  1:2382748/1‑150 (MQ=255)  ggAGGGGATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACgg                                                                                                                                                     >  2:926873/1‑150 (MQ=40)   gAGGGGATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGgc                                                                                                                                                    >  2:307045/1‑150 (MQ=40)     ggggATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGgcgc                                                                                                                                                  >  1:1436995/1‑150 (MQ=255)     ggggATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGgcgc                                                                                                                                                  <  2:1436995/150‑1 (MQ=255)        gATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGaa                                                                                                                                               >  1:1868451/1‑150 (MQ=255)        gATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGaa                                                                                                                                               <  1:307045/150‑1 (MQ=255)         aTGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGTT‑CGCGCCGTACCGCGCGacc                                                                                                                                              <  1:619222/150‑3 (MQ=255)         aTGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAAc                                                                                                                                              >  1:1788659/1‑150 (MQ=255)          tGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACa                                                                                                                                             >  1:1236054/1‑150 (MQ=255)          tGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACa                                                                                                                                             <  1:1805444/150‑1 (MQ=255)          tGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACa                                                                                                                                             >  1:359247/1‑150 (MQ=40)           gCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGTT‑CGCGCCGTACCGCGCGACCac                                                                                                                                            <  2:2257887/150‑2 (MQ=14)            cTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGACCAG‑GTGGT‑CGCGCTGGACGGCGCGAACAgt                                                                                                                                           >  2:21019/1‑150 (MQ=255)              gATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAgtgt                                                                                                                                         <  1:2437596/150‑1 (MQ=255)               aTGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGtgtg                                                                                                                                        >  1:1721698/1‑150 (MQ=255)                     ttCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGa                                                                                                                                  <  1:716036/150‑1 (MQ=255)                     ttCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGa                                                                                                                                  >  2:2078808/1‑150 (MQ=255)                     ttCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGa                                                                                                                                  >  1:2187718/1‑150 (MQ=40)                        agcgcgagcgaCTTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACga                                                                                                                               <  2:1721698/147‑1 (MQ=255)                         aaGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGTGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACgag                                                                                                                              <  2:1236054/150‑1 (MQ=255)                           gcgcgcGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTACGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACgagga                                                                                                                            <  1:1949745/150‑1 (MQ=255)                                cgAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGagcagc                                                                                                                       <  1:224433/150‑1 (MQ=255)                                  aaTTGCGTGACAGTGATGCGCGGGTGAAGAGCGAACTTGTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCgg                                                                                                                     <  2:141827/150‑1 (MQ=255)                                    ttGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATGA‑‑GTGTTGCGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGcc                                                                                                                   <  1:415740/150‑1 (MQ=11)                                    ttGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGcc                                                                                                                   >  1:2321635/1‑150 (MQ=38)                                      gCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGc                                                                                                                 <  2:359247/150‑1 (MQ=37)                                       ccTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGTT‑CGCGCCGTACCGCGCGaccacc                                                                                                                                           <  2:2263187/123‑6 (MQ=14)                                       ccTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGTT‑CGCGCCGTACCGCGCGaccacc                                                                                                                                           >  1:2263187/1‑118 (MQ=14)                                       ccTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCt                                                                                                                >  1:1943613/1‑150 (MQ=38)                                         tGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCTATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAgtg                                                                                                                                          <  1:542736/122‑1 (MQ=255)                                         tGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAgtg                                                                                                                                          >  2:542736/1‑122 (MQ=40)                                          gACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTa                                                                                                             <  1:2078808/150‑1 (MQ=255)                                           aCAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAg                                                                                                            >  1:1320448/1‑150 (MQ=255)                                             aGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAggg                                                                                                          >  1:761510/1‑150 (MQ=255)                                                gATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGt                                                                                                       >  2:1986757/1‑150 (MQ=255)                                                 aTGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGtt                                                                                                      <  2:1287962/150‑1 (MQ=40)                                                  tGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGGACAGTGTGACCGGACGGGGAGCCGCGGCCGCTTTAGGGCGttt                                                                                                     >  2:349152/1‑150 (MQ=255)                                                  tGCGCGGGTGAAGACCGAACCTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGttt                                                                                                     >  1:2234745/1‑150 (MQ=255)                                                    cgcgGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGa                                                                                                   <  1:2215857/150‑1 (MQ=255)                                                      cgGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGTCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGaaa                                                                                                 <  1:2307637/150‑1 (MQ=255)                                                      cgGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGaaa                                                                                                 <  1:21019/150‑1 (MQ=255)                                                         gTGAAGACAGAACGTCTGGCGGGCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTa                                                                                              <  2:1320448/150‑1 (MQ=255)                                                          tGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGACAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTAt                                                                                             <  1:286823/150‑1 (MQ=255)                                                           gAAGCCCGAACTTCTGGCGTCGGCTATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCCGTTCCTACCTTCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATg                                                                                            <  1:1971156/150‑1 (MQ=14)                                                           gAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTa                                                                                              >  1:2074350/1‑148 (MQ=255)                                                           gAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTa                                                                                              <  2:2074350/148‑1 (MQ=255)                                                           gAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATg                                                                                            >  2:613744/1‑150 (MQ=255)                                                             aGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATggg                                                                                          <  2:1422703/150‑1 (MQ=255)                                                               aCCGAACTTCTGGCGGCCGCGATGAATACGGGAGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GGGGT‑CGCGCGGTACGGCGCGGACAGTGTGACCGGCCGAGGAGCAGCGGCCGCTTTAGGGCGTTCGAAACTATGGGgc                                                                                        >  1:62278/1‑150 (MQ=14)                                                                ccGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGgcg                                                                                       >  2:1029044/1‑150 (MQ=255)                                                                  gAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGcc                                                                                     <  1:1029044/150‑1 (MQ=255)                                                                   aaCTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACgag                                                                                                                              <  1:1451603/108‑1 (MQ=255)                                                                   aaCTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACgag                                                                                                                              >  2:1451603/1‑108 (MQ=255)                                                                      ttCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGtt                                                                                 <  2:2382748/150‑1 (MQ=40)                                                                         tGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGG‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCCGTTTAg                                                                              >  2:956942/1‑150 (MQ=255)                                                                          ggcggcCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGGCGCTTTAGGGCGTTTGAAACTATGGGGCGCTAGTTTAGc                                                                             >  1:131850/1‑150 (MQ=255)                                                                             ggcCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCtgt                                                                          <  1:2193287/150‑1 (MQ=40)                                                                             ggcCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCtgt                                                                          <  2:1620390/150‑1 (MQ=40)                                                                             ggcCGCGATGAATACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCGAGTTTAGCtgt                                                                          <  2:1114223/150‑1 (MQ=255)                                                                               ccGCGATGACTACGGGTGCGGCGGGATTCGAACTTGGCGACA‑CTACGCTCTTCCGCTTTGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTCGAAACTATGGGGCGCCAGTTTAGCTGTGa                                                                        <  2:62278/150‑1 (MQ=11)                                                                               ccGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGa                                                                        >  1:1564697/1‑150 (MQ=255)                                                                                  cgATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAg                                                                     <  1:1382037/150‑1 (MQ=255)                                                                                   gATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAgg                                                                    >  2:1827965/1‑150 (MQ=255)                                                                                         tACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGttt                                                              >  1:2148669/1‑150 (MQ=255)                                                                                         tACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGttt                                                              >  2:844934/1‑150 (MQ=255)                                                                                         tACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGttt                                                              >  1:562342/1‑150 (MQ=255)                                                                                          aCGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTc                                                             <  1:613744/150‑1 (MQ=255)                                                                                             gCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTa                                                          <  2:1344386/150‑1 (MQ=255)                                                                                                  ggcggGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGc                                                     <  2:2234745/150‑1 (MQ=255)                                                                                                       gCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCtggt                                                <  2:1868451/150‑1 (MQ=255)                                                                                                         ttgggggTTGGAGCGGGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCtggtgg                                              <  2:131850/144‑1 (MQ=255)                                                                                                         ttCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCtggtgg                                              <  2:1943613/150‑1 (MQ=255)                                                                                                         ttCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCtggtgg                                              <  1:1986757/150‑1 (MQ=255)                                                                                                                  ggAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGcc                                     <  1:1727728/150‑1 (MQ=255)                                                                                                                  ggAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGcc                                     <  2:1788659/150‑1 (MQ=255)                                                                                                                     gCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGTGGGTTTCGTACTCGCCTGGTGGGACTTCGCCGAt                                  <  1:461102/150‑1 (MQ=255)                                                                                                                     gCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGAt                                  >  2:288117/1‑150 (MQ=255)                                                                                                                        aGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCg                               >  1:1351009/1‑150 (MQ=255)                                                                                                                        aGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCg                               <  1:844934/150‑1 (MQ=255)                                                                                                                         gTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCt                                                                            >  1:1948122/1‑104 (MQ=255)                                                                                                                         gTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCt                                                                            <  2:1948122/104‑1 (MQ=255)                                                                                                                         gTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCgg                              <  1:926873/150‑1 (MQ=255)                                                                                                                             cGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTca                          <  2:1351009/150‑1 (MQ=255)                                                                                                                                 gATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATgg                                                                                           <  1:1923724/81‑1 (MQ=39)                                                                                                                                 gATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATgg                                                                                           >  2:1923847/1‑81 (MQ=39)                                                                                                                                 gATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATgg                                                                                           >  2:1923724/1‑81 (MQ=39)                                                                                                                                 gATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATgg                                                                                           <  1:1923847/81‑1 (MQ=39)                                                                                                                                         tgtt‑CGCGCCGTACCGCGCGACCAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGGTTGAAACTATGGGGCGCCAGTTCAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTAatat               >  2:1520658/4‑150 (MQ=25)                                                                                                                                         tggt‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTAatat               >  2:1705321/1‑150 (MQ=255)                                                                                                                                                cgcgGAACGGCGCGAACCGTGTGACCGGACGAGGAGCAGCGGCTGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTNCTCGNCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTAATATACCGcc         <  1:349152/150‑1 (MQ=255)                                                                                                                                                 gcgGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTAATATACCGCCa        <  2:962096/150‑1 (MQ=255)                                                                                                                                                     tACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTAATATACCGCCAGCaa    <  2:761510/150‑1 (MQ=255)                                                                                                                                                       cGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTAATATACCGCCAGCaaaa  <  1:1827965/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                     TGGAGGGGATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTAATATACCGCCAGCAAAA  >  NC_002947/3838406‑3838703

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