Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	941,713	0	T	G	44.8%	125.6 / 34.5	87	V252G (GTG→GGG)	dmsA	dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (11/37);  new base G (39/0);  total (50/37)
Fisher's exact test for biased strand distribution p-value = 6.48e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941569‑941862                                                                                                                                                  |                                                                                                                                                       actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCtt                                                                                                                                                   >  1:771847/1‑150 (MQ=255)   tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGG‑GGGGGGGCtt                                                                                                                                                   >  1:610461/1‑148 (MQ=255)    cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACt                                                                                                                                                <  1:1450982/150‑1 (MQ=255)       ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATc                                                                                                                                             <  2:1023542/150‑1 (MQ=255)       ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATc                                                                                                                                             <  1:195720/150‑1 (MQ=255)        cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCCTTCTATCtt                                                                                                                                           >  1:41846/1‑151 (MQ=255)           gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGaa                                                                                                                                        <  1:1164478/151‑1 (MQ=255)           gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGaa                                                                                                                                        <  2:153750/151‑1 (MQ=255)             aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGG‑GGGGGGGCTTACTAttttttaa                                                                                                                                       >  2:1205943/1‑146 (MQ=255)             aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGG‑GGGGGGGATTACTATTTTTaa                                                                                                                                        >  2:1402339/1‑149 (MQ=255)              aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAg                                                                                                                                     <  2:1261316/151‑1 (MQ=255)                 tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGc                                                                                                                                   <  1:1232627/150‑1 (MQ=255)                 tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCTTTCCATCttta                                                                                                                                         >  1:1319185/1‑142 (MQ=255)                 tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCTTACTATCTTTaaaaag                                                                                                                                    >  1:415197/1‑145 (MQ=255)                  gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCAc                                                                                                                                 <  2:1071682/151‑1 (MQ=255)                    ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGGTTAATATCTTTaaaaag                                                                                                                                    >  2:960005/1‑142 (MQ=255)                         gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCTTTCTATCTTTAAAAGGCACGcc                                                                                                                              >  1:1246746/1‑147 (MQ=255)                          tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGG‑GGGGGGGCTTAATATCTTTAAAAAGCCCGCCAGaa                                                                                                                          >  2:1546929/1‑150 (MQ=255)                           ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGaaa                                                                                                                         <  2:251436/150‑1 (MQ=255)                           ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGAAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGG‑GGGGGGGCTTTCTTTCTTTAAAAAGGACGCCCGaaa                                                                                                                         >  2:1344002/1‑150 (MQ=255)                            tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGACTT‑CACTTTTTAACAGGCAAGCCCGaaa                                                                                                                         >  1:1357519/1‑149 (MQ=255)                                cTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGG‑GGGGGGGCTTTCTATCTTTAACAAGCACGCCAGAAATCTaa                                                                                                                    >  1:313718/1‑150 (MQ=255)                                 tataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGG‑GGGGGGGCTTTCTATCTTTAAAAGGCACGCCAGAAAACTaa                                                                                                                    >  1:1322524/1‑149 (MQ=255)                                  ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGATTACTATCTTTAAAAGGCACGCCAGAAAACTAATGc                                                                                                                 >  1:336087/1‑151 (MQ=255)                                         cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATg                                                                                                         <  2:548929/151‑1 (MQ=255)                                         cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCTTACTATCTTGAAAAAGCACGCCAGAAAACTAATGC‑CCGCATg                                                                                                         >  2:429984/1‑151 (MQ=255)                                             ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGG‑GGGGGGGCTTACTATCTTTAACAAGCCCGCCAGAAATCTAATTC‑CCGCATg                                                                                                         >  2:985708/1‑147 (MQ=255)                                               cTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAATCTAATGC‑CCCCATGatca                                                                                                     >  1:265729/1‑149 (MQ=255)                                                tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGG‑GGGGTGACTTACTATCTTGAAAAGGCACGCCAAAAATTTAAAGC‑CCGCCTGatcatc                                                                                                   >  1:1483879/1‑150 (MQ=255)                                                tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGG‑GGGGGGGCCTTGTATTTTGTACAGGGAGGACCCAAAACTAATGCGCCGCagatcctc                                                                                                    >  2:1191654/1‑143 (MQ=255)                                                  ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGatcatcatc                                                                                                <  1:1337538/151‑1 (MQ=255)                                                    cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGa                                                                                              <  2:282748/151‑1 (MQ=255)                                                     aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCTTACTATCTTGAAAAGGCACGCCCGAAAACTAATGC‑CCGCCTGATCATCATCGa                                                                                              >  1:602173/1‑150 (MQ=255)                                                      gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCTTTCTATCTTTAACAGGCACGCCAGAAAACTAATGC‑CCGCATGATCATCATCGa                                                                                              >  1:814285/1‑149 (MQ=255)                                                         ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCCCATGATCATCATCGATcc                                                                                           >  1:1176662/1‑149 (MQ=255)                                                           cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGACTTACTATTTTTAAAAGGCACGCCCGAAAACTAATGC‑CCGCATGATCATCATCGATcc                                                                                           >  1:1535341/1‑147 (MQ=255)                                                             aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGTGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCta                                                                                     <  1:1061948/151‑1 (MQ=255)                                                             aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCTTAATATCTTTAAAAAGCACGCCAGAAAACTAAATC‑CCCCATGATCATCCTCGATcccc                                                                                         >  2:1220970/1‑145 (MQ=255)                                                               aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCtata                                                                                   <  1:1553432/151‑1 (MQ=255)                                                                  ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCg                                                                                <  1:1281335/151‑1 (MQ=255)                                                                   cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGa                                                                               <  2:952684/151‑1 (MQ=255)                                                                      cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCTTTCTATTTTGAACAAGCACGCCAGAAAACTAATGC‑CCGCATGATCATCATCCATCCGCGCTATACCcaca                                                                             >  2:899195/1‑150 (MQ=255)                                                                      cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCTTACTATCTTTAACAAGCACGCCAGAAAACTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGacac                                                                            >  1:617711/1‑151 (MQ=255)                                                                      cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCCTGATCATCATCGATCCGCGCTATACCGacac                                                                            >  1:441165/1‑151 (MQ=255)                                                                        gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAAAAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCCCGCTATACCGACAcc                                                                           >  1:13788/1‑150 (MQ=255)                                                                         atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTAATATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCGTGATCATCATCGATGCGCGCTATAcc                                                                                 <  2:1189765/143‑1 (MQ=255)                                                                            tCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGC‑GGGGGGGGCTACTATCTTTTATAAGCGCCCCCGAAATTCAATTC‑CCCCATCGTCATTATTAATCCGCGCGATAACGACGccaggg                                                                       >  1:121929/1‑146 (MQ=255)                                                                             cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCCGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGcc                                                                     >  1:527136/1‑151 (MQ=255)                                                                             cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGcc                                                                     <  1:936947/151‑1 (MQ=255)                                                                              gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTTAAAAAGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGcc                                                                     >  1:1141967/1‑150 (MQ=255)                                                                               aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAAACTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGCCACCGGGGCCgg                                                                   >  2:229156/1‑151 (MQ=255)                                                                               aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGACTTACTATCTTTAACAAGGACCCCAGAAATCTAATGC‑CCGCAAGATCATCATCGATCCGCGCCATACCCACACCGGGGCCgg                                                                   >  1:621446/1‑151 (MQ=255)                                                                                aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGCGTTACTATCTTTAACAGGCACGCCGGAAATCTAATGC‑CCCCGTTATCATCATCGATCCGCGCCATACCGACCCCCGTGCCggg                                                                  >  1:659312/1‑151 (MQ=255)                                                                                  aCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGG‑GGGGGGGCTTTATTTCTTGAACAGGCAAGCCCGAAAACTAATGA‑CCGCCTGATCATCATCGATCCGCGCCATACCCACACCGGTGGCggg                                                                  >  2:813568/1‑149 (MQ=255)                                                                                   cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcgc                                                               <  1:496092/151‑1 (MQ=255)                                                                                    aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGcgcg                                                              <  2:1030262/151‑1 (MQ=255)                                                                                     gTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGACTTACTATCTTGAACAGGCACGCCAGAAAACTAATGC‑CCCCATGATCATCATCGATCCGCGCCATACCCACACCGGGGCCGGGcgcg                                                              >  1:737550/1‑150 (MQ=255)                                                                                      tAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGA‑GGGGGGGGTTACTATCTTGAAAAGGCACGCCCGAAATCTAATAC‑CCGCATGATCATCATCGATCCGCGCCATACCGACACCGGGGCCGGGcgcg                                                              >  1:922507/1‑149 (MQ=255)                                                                                        aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGACTTTCTATCTTTGAGAGGCCCGCCCGAAATCTAAAAC‑CCCCAAGAGCCTAATTCATCCGCGCTCTTCCGACACCGGTGCCggggg                                                                >  1:485126/1‑143 (MQ=255)                                                                                           tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCTTTCTATCTTTAAAAAGGACGCCAGAAATCTAATGC‑CCGCCAGATTATCATCGATCCGCCCTATACCCACCCCGGGGCCGGGCGCGAAg                                                           >  1:473168/1‑147 (MQ=255)                                                                                            ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAg                                                      <  1:1174099/151‑1 (MQ=255)                                                                                             gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGt                                                     <  1:338889/151‑1 (MQ=255)                                                                                             gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGACTTACTATCTTTAACAGGCACGCCCGAAATCTAATGC‑CCGCCTGATCATTATCGATCCGCGCTATACCGACACCCGGGCCGGGCGCGAAGATg                                                        >  1:1046998/1‑148 (MQ=255)                                                                                               aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCTTATTATCTTTCACAAGCACGCCAGAAATATAATGG‑CCGCATGATCATCATCGCTCCGCGCTATACCCACACCGGGGCCCGGGGCGAAGATGAGTgg                                                   >  2:638043/1‑151 (MQ=255)                                                                                                gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGa                                                  <  1:1477472/151‑1 (MQ=255)                                                                                                  gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGAt                                                 <  1:360225/150‑1 (MQ=255)                                                                                                    tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGGCTTATTATCTTGAACAAGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCCATCCGCGCTATACCCCCACCGGGGGCGGGGGCGGAAATGGGGGGGTccc                                              >  2:9802/1‑151 (MQ=255)                                                                                                       ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTTCTATCTTGAAAAAGCACGCCAGAAATTTAATGC‑CCGCAAGATCATCATTGCACCGCGCTATACCGACACCCGGGCCGGGCGCGGAGATGAGTGGATcccca                                            >  1:1431230/1‑148 (MQ=255)                                                                                                        tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt                                           <  1:154112/150‑1 (MQ=255)                                                                                                        tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGGCTTACTATCTTGAAAAGGCACGCCCGAAATCTAATGC‑CCGCATGATCATCATCCATCCGCGCTATACCGACAACGGGGGCGGGCGCGGAGATGAGTGGATCCCTAtt                                          >  2:1220134/1‑151 (MQ=255)                                                                                                         ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                          <  1:410939/150‑1 (MQ=255)                                                                                                             ataaCCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTcc                                     <  2:1117861/151‑1 (MQ=255)                                                                                                             ataaCCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGGGACTTACTATCTTGAAAAAGCACGCCCGAAAACTAATGC‑CCCCATGATCATCATCGATTCGCGCTATACCGACACCGGGGCCGGGGGCGGAGATGAGTGGGTCCCCATTCGTcc                                     >  1:286222/1‑151 (MQ=255)                                                                                                               aaCCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCg                                    <  1:34036/150‑1 (MQ=255)                                                                                                                     ggCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc                               <  2:1423322/149‑1 (MQ=255)                                                                                                                     ggCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACa                              <  1:347444/150‑1 (MQ=255)                                                                                                                            cgcgAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGc                      <  2:41846/151‑1 (MQ=255)                                                                                                                                   gAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGt                <  1:258912/150‑1 (MQ=255)                                                                                                                                    aGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa              <  1:492148/151‑1 (MQ=255)                                                                                                                                       ggCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg           <  1:1453087/151‑1 (MQ=255)                                                                                                                                          ggT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt         <  2:1088827/150‑1 (MQ=255)                                                                                                                                          ggT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTc        >  2:1398150/1‑151 (MQ=255)                                                                                                                                              ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg     >  2:1445169/1‑151 (MQ=255)                                                                                                                                               gggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc    <  2:302004/151‑1 (MQ=255)                                                                                                                                               gggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc    >  1:677361/1‑151 (MQ=255)                                                                                                                                                 gTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  2:127233/1‑151 (MQ=255)                                                                                                                                                 gTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCAATGGTTAACGGTCTg      <  2:711697/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                       ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGT‑GGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGC‑CCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941569‑941862

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