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

CGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941564‑941862                                                                                                                                                       |                                                                                                                                                      cGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGT‑GGTGGGGTg                                                                                                                                                      <  1:906282/151‑1 (MQ=255)  gACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGa                                                                                                                                                     <  1:619256/151‑1 (MQ=255)      actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACt                                                                                                                                                   <  2:237810/149‑1 (MQ=255)      actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTa                                                                                                                                                 <  2:61892/151‑1 (MQ=255)      actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACtt                                                                                                                                                  >  2:1189547/1‑150 (MQ=255)         cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACt                                                                                                                                               <  2:176123/150‑1 (MQ=255)            ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGATGGACAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCt                                                                                                                                           <  2:424540/151‑1 (MQ=255)             cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCtt                                                                                                                                          <  1:279535/151‑1 (MQ=255)             cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGGTTTCTATCtt                                                                                                                                          >  2:1226414/1‑151 (MQ=255)               cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGa                                                                                                                                        <  1:1067849/151‑1 (MQ=255)                  aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACa                                                                                                                                     <  1:1182484/151‑1 (MQ=255)                  aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGGCTTTCTAttttttac                                                                                                                                      >  1:252122/1‑146 (MQ=255)                       gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCAc                                                                                                                                <  1:1231664/151‑1 (MQ=255)                        cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACg                                                                                                                               <  2:537666/151‑1 (MQ=255)                         ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGc                                                                                                                              <  2:6220/151‑1 (MQ=255)                            aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGAGGGGTGACTTAATATTTTTaaaatg                                                                                                                                   >  1:136429/1‑139 (MQ=255)                              gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaa                                                                                                                         <  2:657399/151‑1 (MQ=255)                                  gAACTATACCTACGGCGGCTGGGCGGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGGGTTATTATCTTTAAAAAGGCCGCCACaaaa                                                                                                                       >  1:1260295/1‑148 (MQ=255)                                  gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTTCTATCTTTAAAAAGCACGCCAGaaa                                                                                                                        >  2:931934/1‑148 (MQ=255)                                      tataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTTCTATCTTTAACAGGCACGCCAGAAAACTaa                                                                                                                   >  2:540815/1‑149 (MQ=255)                                       ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGTACTTTCAATCTTTAAAAGGCACCGCACAAAACTAAttt                                                                                                                >  2:1093326/1‑149 (MQ=255)                                        taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGAGGTACTATCTTTAAAAGGGACCCCCAAAATCTAAtatc                                                                                                               >  1:193636/1‑148 (MQ=255)                                            tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGGCTTTCTATCTTTAAAAAGGACGCCCGAAAACTAATGccccc                                                                                                            >  2:966258/1‑148 (MQ=255)                                              cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATg                                                                                                         <  2:1231855/151‑1 (MQ=255)                                                 cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGCGGGGGGGCTTTCTATCTTTAAAAAGCACGCCAGAAAACTAATTCCCGcccg                                                                                                         >  1:1045485/1‑145 (MQ=255)                                                     tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGGTTCCTATCTTCAACAGGCACCCCCGAAAACTAATTCCCGCGTTatcat                                                                                                    >  1:10657/1‑149 (MQ=255)                                                        gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGAGTTATTATCTTTAAAAGGCACGCCAGAAAACTAAAGCCCCCATGatcatc                                                                                                   >  1:493536/1‑147 (MQ=255)                                                           gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTACTATCTTGAACAGGCACGCCCGAAAACTAATGCCCGCATGATCATCATCGATc                                                                                            >  2:413390/1‑151 (MQ=255)                                                              ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTTCTATCTTTAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCcgc                                                                                         >  2:425865/1‑151 (MQ=255)                                                                  aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTTCTATCTTTAACAAGCACCGCAGAAATCTAATTCCCGCATGATTATCATCGATCcg                                                                                          >  1:180095/1‑146 (MQ=255)                                                                           cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGAATTACTAACCTGAAAAAGCCCGCCAGAAATCTAATGCCCGCCTGAACATTATTGAACCGCGCcataac                                                                                 >  2:418594/1‑144 (MQ=255)                                                                             gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCg                                                                          <  2:594125/151‑1 (MQ=255)                                                                              atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGGTTACTATCTTGAAAAGGCACGCCCGAAAACTAAAGCCCCCCCGATCCTCATCGATCCGCGCTATACCGACAcc                                                                           >  1:964141/1‑149 (MQ=255)                                                                                atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAAACTAAATCCCGCCAGATCATCATCGATCCGCGCTATAACGACACCgg                                                                         >  2:281704/1‑149 (MQ=255)                                                                                  cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTACTATCTTTTAAAGGCACGCCAGAAAACTAATGCCCGCCTGATTCTCATCGATCCGCGCTATACCCACACCCGTGcc                                                                     >  1:973182/1‑151 (MQ=255)                                                                                    aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGGCTTTCTATCTTGAACCAGGACGCCAGGAACCTAATACCCCCCAGATTCTCATCGATCCGCGCTATTCATCCACCGGGGcccg                                                                   >  2:473134/1‑149 (MQ=255)                                                                                    aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGCGCCgg                                                                   >  2:95156/1‑151 (MQ=255)                                                                                    aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGGTTACTATCTTTAAAAAGCACGCCCGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCCCTACCGCCACCGGTg                                                                       >  1:818345/1‑147 (MQ=255)                                                                                    aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTACTATCTTTAACAGGCCCGCCCGAAAACTAATGCCCGCATGATTATCATCCATCCGCGCCATAACGACACCGGGGCCgg                                                                   >  2:1054027/1‑151 (MQ=255)                                                                                    aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGAGTTACTATTTTTAAAAAGCACCCCCGAAAACTAATGCCCGCCAGATCATCATCGACTCGCGGTAAACCCACACCCGGGcccg                                                                   >  1:1254542/1‑149 (MQ=255)                                                                                        cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcg                                                                <  1:1142809/150‑1 (MQ=255)                                                                                             aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGa                                                          <  2:1246994/151‑1 (MQ=255)                                                                                             aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTTATATTTTTAACAGGCCCCCCAGAAAACTTATGCCCGCCGGATCATTTTCGATTCGCGCCTTACCGATCCccgggccggg                                                                  >  1:157437/1‑143 (MQ=255)                                                                                                    aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTg                                                    <  2:816872/150‑1 (MQ=255)                                                                                                      tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGa                                                  <  1:544764/150‑1 (MQ=255)                                                                                                       gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCCATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATc                                                >  2:1220621/1‑151 (MQ=255)                                                                                                          gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTACTATCCTTAACAAGCACGCCAGAAAACTAATGGCCGCATGATCATCATCCATCCCCGCTATAACGACACCGGTGGCGGGCGCGAAGATGAGTGGGTccc                                              >  1:1246761/1‑150 (MQ=255)                                                                                                              ggTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGCGGGGGGGCCTAATAACTTGAAAACGCGCGCCACGAAAATAAAAACCGCCCGACTATTATCGATCCCCGCTATACCGACACCCGGGGCGGGGGGGGAGACGCGGGGATCCCTaat                                          >  1:1189137/1‑148 (MQ=255)                                                                                                                taataaCCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATTTTTAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCCACACCGGGGCCGGGGGGGAAGATGGGTGGGTCCCCAtt                                          >  2:328999/1‑148 (MQ=255)                                                                                                                taataaCCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTTAAAAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCCGTGCCGGGCGCGAAGATGGGGGGGTCCCCAtttg                                        >  2:1053556/1‑148 (MQ=255)                                                                                                                     aCCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTAATATCTTGAACAGGCAAGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGAAGATGAGTGGATCCCTATTCGTCCggg                                  >  2:228606/1‑151 (MQ=255)                                                                                                                      cccTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                                  <  1:450198/150‑1 (MQ=255)                                                                                                                       ccTGGCGAAACGCGAATGAGTGGCGGGTGGGGGGA‑TTACTATCTTGAACAAGGGCGCCCGAAAACTAATACCCCCATGATCATCATCGAACCGCGCCATAACGACACCGGGGCCCGGGGCGGAAAAGTGAGGGTCCCCCTTCGTcc                                     >  1:1141003/1‑146 (MQ=255)                                                                                                                                aCGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCg                       <  1:1118199/151‑1 (MQ=255)                                                                                                                                aCGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCg                       <  1:32194/151‑1 (MQ=255)                                                                                                                                aCGCGAATGAGTGGC‑GGTGGGGGGACTTTCTATTTTTAAAAAGGACGCCCGAAAAATAATGGCCGCCTGATCATCATCGATCCGCGCTATAACGACCCCCGGGGCGGGGGCGGAGATGGGTGGGTCCCTATTTGTCCGGGGACAGAAGcc                        >  1:927558/1‑150 (MQ=255)                                                                                                                                  gcgAATGAGTGGC‑GGTGGGGGGGCTTACTCTCTTGAACAAGCACGCCAGAAAAATAATGCCCGCACGGTTATTATCGAACCGCGCTATACCGACACCCGTGCCGGGGGCGAAGATGAGTGGGTCCCTATTCGTCCGGGGACAGATGcc                        >  1:707458/1‑148 (MQ=255)                                                                                                                                     aaTGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTg                  <  2:291231/151‑1 (MQ=255)                                                                                                                                     aaTGAGTGGC‑GGTGGGGGGGCTTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCGGTGCCGGGCGCGAAGATGAGGGGATCCCTATTCGTCCGGGAACAGATGCCGCACTg                  >  1:400332/1‑151 (MQ=255)                                                                                                                                            ggC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg           <  1:986508/151‑1 (MQ=255)                                                                                                                                              c‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt         >  1:1243767/1‑151 (MQ=255)                                                                                                                                                ggTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTc        >  2:258641/1‑151 (MQ=255)                                                                                                                                                ggTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTc        <  1:281698/151‑1 (MQ=255)                                                                                                                                                 gTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt       >  1:453743/1‑151 (MQ=255)                                                                                                                                                  tGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg      >  1:372393/1‑151 (MQ=255)                                                                                                                                                   ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg     >  2:1193137/1‑151 (MQ=255)                                                                                                                                                   ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg     >  2:3142/1‑151 (MQ=255)                                                                                                                                                    gggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc    <  1:205376/151‑1 (MQ=255)                                                                                                                                                      gTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  2:364310/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941564‑941862

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