Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,838,550	G→C	5.6%	R701P (CGG→CCG)	PP_3388 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,838,550	0	G	C	5.6%	276.5 / 7.0	90	R701P (CGG→CCG)	PP_3388	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/47);  new base C (2/3);  total (40/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

CTGCTGGAGGGGATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTAATATACCGCCA  >  NC_002947/3838402‑3838697                                                                                                                                                       |                                                                                                                                                    gtgctgGAGGGGATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCgg                                                                                                                                                    <  2:774844/149‑1 (MQ=255)     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:307045/150‑1 (MQ=255)            gATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGaa                                                                                                                                         >  1:1868451/1‑150 (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:359247/1‑150 (MQ=40)              tGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACa                                                                                                                                       <  1:1805444/150‑1 (MQ=255)               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                                                                                                                            >  2:2078808/1‑150 (MQ=255)                         ttCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGa                                                                                                                            <  1:716036/150‑1 (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                                                                                                                                     >  1:2263187/1‑118 (MQ=14)                                           ccTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGTT‑CGCGCCGTACCGCGCGaccacc                                                                                                                                     <  2:2263187/123‑6 (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                                                                                        <  2:2074350/148‑1 (MQ=255)                                                               gAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTa                                                                                        >  1:2074350/1‑148 (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                                                                                                                        >  2:1451603/1‑108 (MQ=255)                                                                       aaCTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACgag                                                                                                                        <  1:1451603/108‑1 (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                                                                    <  2:1620390/150‑1 (MQ=40)                                                                                 ggcCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCtgt                                                                    <  1:2193287/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:562342/1‑150 (MQ=255)                                                                                             tACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGttt                                                        >  2:844934/1‑150 (MQ=255)                                                                                             tACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGttt                                                        >  1:2148669/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                               <  2:1788659/150‑1 (MQ=255)                                                                                                                      ggAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGcc                               <  1:1727728/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                                                                      <  2:1948122/104‑1 (MQ=255)                                                                                                                             gTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCt                                                                      >  1:1948122/1‑104 (MQ=255)                                                                                                                             gTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCgg                        <  1:926873/150‑1 (MQ=255)                                                                                                                                 cGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTca                    <  2:1351009/150‑1 (MQ=255)                                                                                                                                     gATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATgg                                                                                     >  2:1923847/1‑81 (MQ=39)                                                                                                                                     gATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATgg                                                                                     >  2:1923724/1‑81 (MQ=39)                                                                                                                                     gATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATgg                                                                                     <  1:1923724/81‑1 (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)                                                                                                                                                       |                                                                                                                                                    CTGCTGGAGGGGATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAG‑GTGGT‑CGCGCGGTACGGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTAATATACCGCCA  >  NC_002947/3838402‑3838697

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