Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7130TG37.3% 95.9 / 16.1 59V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (6/31);  new base G (22/0);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 4.20e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941563‑941860
                                                                                                                                                      |                                                                                                                                                      
gCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                                                                                        <  2:639081/151‑1 (MQ=255)
gCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                                                                                        <  1:257429/151‑1 (MQ=255)
 cGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGTCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGTCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                                                                                        <  2:763705/150‑1 (MQ=255)
 cGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTg                                                                                                                                                       <  1:294408/151‑1 (MQ=255)
      actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                                   >  2:279286/1‑150 (MQ=255)
         cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTa                                                                                                                                               <  2:1139037/151‑1 (MQ=255)
             cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                                                                                        <  2:14370/138‑1 (MQ=255)
             cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAt                                                                                                                                              <  1:1343327/148‑1 (MQ=255)
              gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTg                                                                                                                                         <  1:1023442/151‑1 (MQ=255)
                 cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTaa                                                                                                                                       >  2:507758/1‑150 (MQ=255)
                      tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTACTATC‑TTTAACa                                                                                                                                     >  1:522085/1‑147 (MQ=255)
                              gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGaa                                                                                                                         <  2:718510/151‑1 (MQ=255)
                              gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGaa                                                                                                                         <  1:403502/151‑1 (MQ=255)
                                ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTTAACAGGCAAGCCCGaaa                                                                                                                        >  2:834741/1‑150 (MQ=255)
                                  gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATc                                                                                                                      <  1:266015/150‑1 (MQ=255)
                                  gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATc                                                                                                                      <  2:692726/150‑1 (MQ=255)
                                   aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑a                                                                                                                   <  2:294028/151‑1 (MQ=255)
                                              cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATg                                                                                                        <  2:240056/151‑1 (MQ=255)
                                              cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATg                                                                                                        <  1:949322/151‑1 (MQ=255)
                                              cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTTTC‑TTTAAAAAGGACGCCCCAAAACT‑AATCTCCGcgc                                                                                                         >  1:1127033/1‑148 (MQ=255)
                                               gccggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGCGGCGGTGGGGTGGGTTACTCTC‑TTCAAAACGCCCGGCACaaa                                                                                                                        >  1:274657/3‑135 (MQ=255)
                                                gcggcTGGGCAGATGGCAGCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATACCCCTGGCGAAACGCGAATGAGTGGGGGTGGGGTGACTTACTATC‑TAGAACAAGCACGCCAGCAATCT‑AATGCCCGCagg                                                                                                        <  1:581056/149‑3 (MQ=255)
                                                gcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATT‑TTTAACAGGCACGCCAGAAATCT‑AAAGCCCGCATGa                                                                                                       >  2:105976/1‑150 (MQ=255)
                                                  ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTGAAAAGG‑ACGCCAGAAATAT‑AATGACCGCAtta                                                                                                       >  1:44197/1‑146 (MQ=255)
                                                        gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCg                                                                                              <  1:1277414/151‑1 (MQ=255)
                                                        gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTTTC‑TTTAAAAAACACGCCTGAAATCC‑AATGCCCGCATGatccaca                                                                                                 >  1:61856/1‑144 (MQ=255)
                                                         cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTACTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTCCATC‑TTTAAAAGGCACGCCCGGAAACTAAATGCCCCCA‑Gatcat                                                                                                   >  1:1226739/1‑145 (MQ=255)
                                                          aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTTCTATC‑TTGAAAAGGCCCGGCAGAAAACT‑AATGCCCGCACGAACATCATCGAt                                                                                            >  2:569681/1‑151 (MQ=255)
                                                          aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTTAACAGGCACGCCAGAAAACT‑AATTCCCGCCTGATCATCATCGa                                                                                             >  2:497921/1‑150 (MQ=255)
                                                          aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTTCTATC‑TTTAAAAGGCACGCCAGAAAACT‑AATGTCCGGCTGATCATTATCGAt                                                                                            >  2:1333958/1‑151 (MQ=255)
                                                           gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGTCTTACTATT‑TTTAACAGGCACACCAGAAAATT‑AATGCCCTCAAAattatcatctta                                                                                            >  2:46425/1‑147 (MQ=255)
                                                           gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGGAATAATATC‑TTTAACAGACCCGGCCGCAATAT‑AATACCCGCCTGGTCATCAACCATc                                                                                           >  1:833047/1‑151 (MQ=255)
                                                               gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCcgcg                                                                                       <  1:1274342/151‑1 (MQ=255)
                                                                 aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCgcgc                                                                                      <  1:582717/150‑1 (MQ=255)
                                                                   cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCtat                                                                                   <  2:1043037/151‑1 (MQ=255)
                                                                       ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTTAAAAGGCACGCCAGAAAAAT‑AATGCCCGCATGATCATCATCGATCCGCGCtata                                                                                  >  2:819419/1‑148 (MQ=255)
                                                                              atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCgg                                                                        <  2:459660/151‑1 (MQ=255)
                                                                              atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTGAACAGGCACGCCCGAAATCT‑AATGCCCGCATTATCATCATCGATTCGCGCTATACCGCCACCgg                                                                        >  2:513287/1‑151 (MQ=255)
                                                                                   gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATC‑TTTAACAAGCACGCCAGAAATCT‑AAATCCCGCCCGGTCATCATCGATCCGCGCTATACCGACACCGGGGc                                                                     >  2:39594/1‑149 (MQ=255)
                                                                                    aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCgg                                                                  <  2:914211/151‑1 (MQ=255)
                                                                                    aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTTAACAGGCACGCCAGAAAACT‑AATGCCCCCATGATCATCATCGATCCCCGCCATACCGACACCGGGGGCgg                                                                  >  2:285121/1‑151 (MQ=255)
                                                                                     aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCggg                                                                 <  1:468075/151‑1 (MQ=255)
                                                                                          gTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGa                                                            <  1:506497/151‑1 (MQ=255)
                                                                                              gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                        <  1:975335/151‑1 (MQ=255)
                                                                                                tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTTAACAAGCACGCCCGAAATCT‑AATGCCCCCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGGAGATg                                                       >  2:487320/1‑150 (MQ=255)
                                                                                                   tAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATC‑TTTAACAAGCACGCCCGAAAACT‑AATGCCCGCATGATCATCATCGATCCGCGCTATAACGACACCGGGGCCGGGCGCGGAGATGAg                                                     >  1:1344459/1‑149 (MQ=255)
                                                                                                     gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGa                                                 <  2:627560/151‑1 (MQ=255)
                                                                                                         tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATC‑TTGCACCAGCACGCCAGAAATCT‑TATGCCCGCATGATAATCATTGATTCGCGCTATAACCCCCCCCGCGCCGGGCGGGAAGATGAGTGGATccc                                             >  2:116230/1‑151 (MQ=255)
                                                                                                           tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGCTCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGGGGATccc                                             >  2:135888/1‑149 (MQ=255)
                                                                                                             tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTTCTAGC‑TGGGACAGGCACGCGCGAAATCT‑AATGCGCGCAAGATCATCATCGATCCGCGGTATACTAACACCTGGTCCAGGCGCGAATGTGAGTGTGACCCTAt                                          <  1:1329537/150‑1 (MQ=255)
                                                                                                              ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                         <  2:365007/150‑1 (MQ=255)
                                                                                                              ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                         <  2:804394/150‑1 (MQ=255)
                                                                                                              ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                        <  1:1267428/151‑1 (MQ=255)
                                                                                                                    aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATT‑TTTAACAGGCACGCCAGAAATCT‑AATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGGTCCCCATTCGT‑CCgg                                 >  2:705484/1‑151 (MQ=255)
                                                                                                                          ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAg                           <  1:1108393/151‑1 (MQ=255)
                                                                                                                            cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCAGCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGGGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGAt                         <  1:1186063/151‑1 (MQ=255)
                                                                                                                             gAAACGCGAATGAGTGGCGGAGGGGGGGCTTTCTTTC‑TTGAACAGGAACGCCAGGAAATT‑AAAGCCCGCCAGATCATCACCGGTCCCCGCCATACCGACAACCCTGGCGGGCGCGAAGAAGAGGGGATCCCCATTCGT‑CCGGGGACAGGTg                        >  1:897366/1‑151 (MQ=255)
                                                                                                                                   cgAATGAGTGGCGGTGGGGGGACTTAATATC‑TTGAAAAGGGACGCCAAAAATCT‑TAAGCCCGCAAGATCATCATCGATCCGCGCTATAACCACACCCGGGCCGGGGGCGAAGAAGAGGGGGTTCCTTTTCGTCCCGGGaaaaatggc                       >  2:858578/1‑140 (MQ=255)
                                                                                                                                       tGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGACGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTgg               <  2:624549/150‑1 (MQ=255)
                                                                                                                                          gTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGGTTaa           <  2:1323540/151‑1 (MQ=255)
                                                                                                                                          gTGGCGGTGGGGGGGCTTACTATC‑TTGAACAAGCACGCCAGAAAACT‑AATGCCCCCATGATCATCATCGATCCGCGCTATACCCACACCGGGGCCGGGGGCGAAGATGAGTGGATCCCTCTTTGT‑CCGGGGACAGATGCCGCCCTGGGTa            >  2:920542/1‑150 (MQ=255)
                                                                                                                                               ggTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGGTTAACGGTc      >  2:1410073/1‑151 (MQ=255)
                                                                                                                                                   gggTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc  <  2:1331230/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCT‑AATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941563‑941860

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