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

CGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGG  >  NC_000913/941564‑941859                                                                                                                                                      |                                                                                                                                                      cGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTg                                                                                                                                                      <  2:1198278/151‑1 (MQ=255)  gACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGa                                                                                                                                                     <  1:456321/151‑1 (MQ=255)         cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCt                                                                                                                                               >  2:903636/1‑150 (MQ=255)             cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑tt                                                                                                                                         >  1:134074/1‑151 (MQ=255)              gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtttccattttt                                                                                                                                         >  1:343279/1‑142 (MQ=255)                 cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAAc                                                                                                                                     <  2:1249491/151‑1 (MQ=255)                  aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTtctttc‑ttttac                                                                                                                                     >  1:827588/1‑146 (MQ=255)                       gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATC‑TTTAACAGGCAc                                                                                                                               >  2:465271/1‑151 (MQ=255)                          gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGcc                                                                                                                            <  2:288069/151‑1 (MQ=255)                           aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCa                                                                                                                           <  1:83402/151‑1 (MQ=255)                           aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATC‑TTGAACAAGCACGccc                                                                                                                           >  1:1171384/1‑150 (MQ=255)                            aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTTTC‑TTTAAAAAGCACGcccg                                                                                                                          >  2:919175/1‑149 (MQ=255)                              gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGaa                                                                                                                        <  1:1097120/151‑1 (MQ=255)                              gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGaa                                                                                                                        <  2:1366206/151‑1 (MQ=255)                                  gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCt                                                                                                                    <  2:182399/151‑1 (MQ=255)                                  gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCt                                                                                                                    <  2:1012841/151‑1 (MQ=255)                                  gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTATC‑TTTAAAAGGCACGCCAGAAAACt                                                                                                                    >  1:968268/1‑151 (MQ=255)                                   aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTTAACAAGCACCCCAGAAAACt                                                                                                                    >  1:329486/1‑150 (MQ=255)                                    aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑aa                                                                                                                 <  1:1022548/151‑1 (MQ=255)                                    aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTGAAAAGG‑ACGCCCGAAAACT‑a                                                                                                                  >  1:1283916/1‑150 (MQ=255)                                       ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AAt                                                                                                                <  1:640900/149‑1 (MQ=255)                                       ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGACTACTATC‑TTTAACAGGCAAGCCAGAAAACT‑AATGc                                                                                                              >  1:733331/1‑151 (MQ=255)                                           cTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTTCTATT‑TTTAAAAGGCACGCCCGAAAACT‑AATGCCCGc                                                                                                          >  2:1285115/1‑151 (MQ=255)                                            tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATC‑TTTTAAAAGCACGCCCGAAAAAT‑AATTcccccc                                                                                                         >  1:676353/1‑148 (MQ=255)                                              cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATg                                                                                                       <  2:5100/151‑1 (MQ=255)                                               ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTTAACAGGCACGCCAGAAATCT‑AATTCCCGCATg                                                                                                       >  1:326669/1‑150 (MQ=255)                                                 cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATT‑TTTAAAAAGCCCCCCAGAAAACTAAATGccccccggac                                                                                                     >  2:1204432/1‑144 (MQ=255)                                                  ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGatca                                                                                                   <  2:888378/151‑1 (MQ=255)                                                  ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTGAACAAGCACGCCCGAAAACT‑AATGCCCCCATGatc                                                                                                    >  1:425321/1‑150 (MQ=255)                                                   gcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGatcat                                                                                                  <  1:894539/151‑1 (MQ=255)                                                        gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCg                                                                                             <  2:1232645/151‑1 (MQ=255)                                                         cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTGAACAGGGACGCCCGAAAACT‑AAAGCCCCCCTGATCATCATCg                                                                                             >  2:326194/1‑150 (MQ=255)                                                              ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCcg‑c                                                                                      <  2:654718/151‑1 (MQ=255)                                                                cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCg‑cgc                                                                                    <  1:645798/151‑1 (MQ=255)                                                                cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAACAAGCACGCCCGAAATCT‑AATGCCCCCCTGATCATCATCCATCCg‑cgc                                                                                    >  1:309378/1‑151 (MQ=255)                                                                 aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATACCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCCCGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCt                                                                                   <  2:438367/151‑1 (MQ=255)                                                                    aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCtata                                                                                <  2:1121855/151‑1 (MQ=255)                                                                    aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTATTATC‑TTTAAAAGGCACGCCCGAAAACT‑AATGCCCGCCTGATCATCATCGATCCG‑CGCtat                                                                                 >  2:447189/1‑150 (MQ=255)                                                                         gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAAAAAGGACGCCAGAAAACT‑AATGCCCGCCTGATCATCATAGATCCG‑CGCTATAACGac                                                                           >  2:869588/1‑151 (MQ=255)                                                                              atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGATTTACTATC‑TTGAACAGGCACGACAGAAATAA‑AATTGCCGCCGGATAAACATACATCCG‑CGCCCTACAGacacacg                                                                      >  2:413082/1‑148 (MQ=255)                                                                              atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTGAACAAGCACGCCAGAAAACT‑AATGCCCGCATGATCATCATCGATCCG‑CGCCATACCGACAcccg                                                                      >  1:56123/1‑149 (MQ=255)                                                                                   gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCg                                                                 <  2:1062524/151‑1 (MQ=255)                                                                                   gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCg                                                                 <  1:1353124/151‑1 (MQ=255)                                                                                    aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCCTTCTATC‑CTTAACAAGCACGCCACAAAACA‑AATGCCCCCATTAACAACATCGAACCG‑CGCTTTACCCACACCGGTGCCgg                                                                >  2:1160817/1‑151 (MQ=255)                                                                                      aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGgc                                                              <  2:558573/151‑1 (MQ=255)                                                                                             aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAAATTCTCTC‑TTGAAAAGGCACGCCAGAAATAT‑AATGCCCGCCTGATCCTCCTTGGTCCGCCGCATTACC‑ACACccgggccgggggggg                                                          >  2:873157/1‑143 (MQ=255)                                                                                              gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                      <  2:333317/151‑1 (MQ=255)                                                                                               cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATC‑TTTAAAAAGCACGCCAGAAAAAT‑AATGGCCGCATGATCATTATCGATCCG‑CGCTATACCCGCACCCGGGCCGGGGGCGAAGATg                                                     >  2:1354576/1‑151 (MQ=255)                                                                                                 ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGa                                                    <  2:745929/150‑1 (MQ=255)                                                                                                 ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAg                                                   <  2:1164776/151‑1 (MQ=255)                                                                                                 ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTAGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAg                                                   <  2:392341/151‑1 (MQ=255)                                                                                                    aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTg                                                 <  1:556338/150‑1 (MQ=255)                                                                                                            ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTa                                         <  2:858686/150‑1 (MQ=255)                                                                                                              ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                       <  2:275146/150‑1 (MQ=255)                                                                                                              ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                      <  1:1164951/151‑1 (MQ=255)                                                                                                                 aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAAAAGGCACGCCCGAAATCT‑AATTCCCGCCTGATCATCATCGATCCG‑CGCTATACCCACACCGGGGCCGGGGGCGAAGATGAGTGGATCCCTATTCGTc                                   >  1:1222805/1‑151 (MQ=255)                                                                                                                  ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTGAAAAAGCACGCCAGAAAACT‑AATGCCCGCAAGATCATCATCGATCCG‑CGCTATACCGACAACCGGGCCGGGCGCGGAAATGAGTGGGTCCCCATTCGTcc                                  >  2:249605/1‑151 (MQ=255)                                                                                                                   taaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAACAAGCACGCCAGAAATCT‑AATGCCCGCCTGATCATCATCGATCCG‑CGCTATACCCACACCGGGGCCGGGGGCGAAGATGAGTGGGTCCCTATTCGTCCg                                 >  1:1105717/1‑151 (MQ=255)                                                                                                                        cTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTTAAAAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGGTCCCTATTCGTCCGGGAAc                            >  2:1011298/1‑151 (MQ=255)                                                                                                                          ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc                            <  1:151695/149‑1 (MQ=255)                                                                                                                          ggCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGAc                            >  1:87597/1‑149 (MQ=255)                                                                                                                             gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATg                       <  2:306580/151‑1 (MQ=255)                                                                                                                              aaaCGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTTAAAAAGCACGCCAGAAAACT‑TATGGCCGCATGATCATCATCGATCCC‑CGCTATAACGACACCCGTGCCGGGCGCGAAGATGAGTGGATCCCTATTTGTCCGGGGACAGATGc                      >  2:51734/1‑151 (MQ=255)                                                                                                                                aCGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCg                    <  2:42858/151‑1 (MQ=255)                                                                                                                                aCGCGAATGAGTGGCGGTGGGGGGAGTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑TATGGCCGCATGATTATCATCGATCCC‑CGCTATACCCACACCGGGGCCGGGCGCGGAGATGAGGGGATCCCCATTCGTCCGGGGACAGATGGCg                    >  1:242546/1‑151 (MQ=255)                                                                                                                                   cgAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCAc                 <  1:591334/151‑1 (MQ=255)                                                                                                                                     aaTGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTg               <  1:1079237/151‑1 (MQ=255)                                                                                                                                            ggCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg        >  1:766933/1‑151 (MQ=255)                                                                                                                                              cGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt      >  1:1083582/1‑151 (MQ=255)                                                                                                                                               ggTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTc     >  1:778177/1‑151 (MQ=255)                                                                                                                                                gTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt    >  1:266479/1‑151 (MQ=255)                                                                                                                                                  ggggTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg  <  1:874183/151‑1 (MQ=255)                                                                                                                                                  ggggTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg   <  2:329473/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCG‑CGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGG  >  NC_000913/941564‑941859

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