Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7130TG40.0% 91.7 / 28.2 65V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (2/37);  new base G (26/0);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 3.77e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCT  >  NC_000913/941563‑941857
                                                                                                                                                      |                                                                                                                                                  
gCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                                                                                    <  2:829791/151‑1 (MQ=255)
 cGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTg                                                                                                                                                   <  2:612680/151‑1 (MQ=255)
 cGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTg                                                                                                                                                   <  1:1004948/151‑1 (MQ=255)
      actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTa                                                                                                                                              <  1:488104/151‑1 (MQ=255)
       ctcctcCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGGAGTGTTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGatgta                                                                                                                                              >  1:1051060/1‑146 (MQ=255)
            ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATC‑t                                                                                                                                       >  2:996390/1‑151 (MQ=255)
             cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑tt                                                                                                                                      <  2:661609/151‑1 (MQ=255)
               cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTg                                                                                                                                     <  2:527887/150‑1 (MQ=255)
               cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGCCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGa                                                                                                                                    <  1:349482/151‑1 (MQ=255)
                 cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTATC‑TTTaa                                                                                                                                   >  2:507310/1‑150 (MQ=255)
                      tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGc                                                                                                                              <  2:738268/150‑1 (MQ=255)
                       gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑c                                                                                                                           <  2:683265/151‑1 (MQ=255)
                        cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑Cg                                                                                                                          <  2:477707/151‑1 (MQ=255)
                            aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATT‑TTTaaaaaggccgcc                                                                                                                         >  1:810681/1‑139 (MQ=255)
                              gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGaa                                                                                                                    <  1:670155/151‑1 (MQ=255)
                                   aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGGTTACTATCTTTGAAAAGG‑A‑CGCCCGAAAACt                                                                                                                >  1:992448/1‑150 (MQ=255)
                                   aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATT‑TTGAAAAAGGA‑CGCCAGaaa                                                                                                                   >  1:563618/1‑147 (MQ=255)
                                         aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGCTTACTATC‑TTTAAAAGGGA‑CGCCCGAAATCTaaaagcc                                                                                                         >  2:930978/1‑146 (MQ=255)
                                          ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTATC‑TTTAAAAGGCA‑CGCCCGAAAACTaaaaccc                                                                                                         >  1:1205553/1‑145 (MQ=255)
                                            tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGCTATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGGAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTCCATC‑TTGTAAAAGCA‑CGGCACAAAAATAATGcc                                                                                                          >  1:479519/1‑147 (MQ=255)
                                            tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGACCTACTATT‑TTTAAAAGGGACCCCCAAAAATTTAAtcacgcc                                                                                                       >  2:668401/1‑145 (MQ=255)
                                             aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCAt                                                                                                     <  1:194183/151‑1 (MQ=255)
                                              cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATg                                                                                                    <  1:22016/151‑1 (MQ=255)
                                               ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGTCCTACTAAC‑TTTAAAAGGCA‑CCCCACAAAAATAATGTCCGCAt                                                                                                     >  1:519234/1‑149 (MQ=255)
                                               ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTGAAAAAGGA‑CGCCCGAAATCTAAAGCCCGCCTg                                                                                                    >  2:494676/1‑150 (MQ=255)
                                               ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATC‑TTGAAAACGCA‑CGCCAGAAAATTAATGTCCCCATg                                                                                                    >  2:626185/1‑150 (MQ=255)
                                                      gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGatcatcat                                                                                            <  2:726873/151‑1 (MQ=255)
                                                      gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGatcatca                                                                                             <  2:476335/150‑1 (MQ=255)
                                                      gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGatcatca                                                                                             <  2:952560/150‑1 (MQ=255)
                                                      gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGatcatca                                                                                             <  2:50982/150‑1 (MQ=255)
                                                        gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                          <  2:466087/151‑1 (MQ=255)
                                                        gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAAGGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                          <  1:893896/151‑1 (MQ=255)
                                                              ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACGTACTATT‑TTTAACAAGGA‑CGCCCGAAAACTAATGGCCGccaggtcattcatga                                                                                         >  1:751099/1‑132 (MQ=255)
                                                                 aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCt                                                                                 <  1:53264/151‑1 (MQ=255)
                                                                    aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCtata                                                                              <  1:369736/151‑1 (MQ=255)
                                                                      cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTACTATC‑TTTAAAAAGCA‑CCCCAGAAAACTAATACCCGCATGGTCATCatccatcccccccat                                                                               >  1:202168/1‑141 (MQ=255)
                                                                       ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATC‑TTGAACAAGCA‑CGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCCATAcc                                                                            >  1:34488/1‑150 (MQ=255)
                                                                                atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACGCCGGt                                                                   <  1:206033/150‑1 (MQ=255)
                                                                                  cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGc                                                                 <  2:341385/150‑1 (MQ=255)
                                                                                       aCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAACCAGCA‑CGCCCGAAAACTAATGCCCGCACGATCCTCATCGATCCGCGCTTTACCCACAACCGTGCCggg                                                             >  2:878650/1‑149 (MQ=255)
                                                                                        cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcgc                                                          <  1:19803/151‑1 (MQ=255)
                                                                                            aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAg                                                      <  1:517336/151‑1 (MQ=255)
                                                                                            aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGCTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAAAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAg                                                      <  1:207719/151‑1 (MQ=255)
                                                                                             aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGa                                                     <  2:921715/151‑1 (MQ=255)
                                                                                                tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGCATACTATC‑TTGAAAAAGCG‑CCCCCTAAAATTAAAGGCCGCCTCAATAACATCATTCCCCCCTATTCCGCCACCCGCGGCGGGCGGGGAGAt                                                    >  1:872630/1‑149 (MQ=255)
                                                                                                 ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATg                                                   <  1:1114970/149‑1 (MQ=255)
                                                                                                   tAGTGCTGTTTGGTCATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAGC‑TTGAACAGGCG‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGACCCGCGCTATATCGACACCGGTGCCGGTCGCGAAGATgt                                                  <  1:198962/148‑2 (MQ=255)
                                                                                                    aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATT‑TTTAAAAGGCA‑CCCCAGAAAACTAATGCCCCCCTGATCATCATCGATCCGCGCCATACCCACACCCGGGCCGGGCGCGGAGATGAGTgg                                              >  2:1001614/1‑151 (MQ=255)
                                                                                                            ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAAACTAATGCCCCCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGGTCCCTAt                                      >  2:1079251/1‑151 (MQ=255)
                                                                                                            ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGCTTAACATC‑TTTAACAAGCA‑ACCCAGGAATTTAATGACCGCATGATCATCATCGATTCGCGCTTTACCCCCACCGGGGGCGGGCGCGAAGATGAGGGGaaacc                                         >  1:1044527/1‑144 (MQ=255)
                                                                                                                taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTTTC‑TTTAACAGGCA‑CGCCAGAAAACTAATGCCCGCCTGAACATCATCGATCCGCGCTATACCGACACCCGTGCCGGGCGCGGAGATGAGGGGGTCCCTAtt                                     >  2:1173760/1‑148 (MQ=255)
                                                                                                                taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTGAAAAAGCA‑CGCCAGAAAACTAATGCCCGCATGATCATCATCGAACCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGGGGATCCCCAtt                                     >  2:703929/1‑148 (MQ=255)
                                                                                                                taataaCCCTGGCGAAACGCGAATGAGTGGCGGTCGGGGGGATTAATATC‑TTGAAAAGGCA‑CGCCAGAAAACTAATGCCCGCCTGATCATCTTCCGTCCGCGCCATAACGACACCGGGGCCGGGGGCGAACATGAGTGGGTcc                                          >  1:718815/1‑143 (MQ=255)
                                                                                                                    aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATT‑TTTAAAACGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCAATGAACCGCGCCATACCCACACCCGGGCCGGGCGGGAAGATGAGGGGATCCCTATTCGTCCg                               >  2:532170/1‑150 (MQ=255)
                                                                                                                     aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                             <  1:996337/151‑1 (MQ=255)
                                                                                                                      cccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGa                            <  2:85972/151‑1 (MQ=255)
                                                                                                                         tGGCGAAACGCGAATGAGTGGCGGAGGGGGTAGTTAATATC‑TTTAACAGGGA‑CGCCCGAAATCTAATGCCCCCACGGTCCTCAACCATCCTCCCTCTACCCGCACCGGTGACCGGCGCGGAGATGAGGGGGGCCCTATTCCTCCgcggac                          >  1:1179668/1‑145 (MQ=255)
                                                                                                                          ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACa                         <  1:1025335/150‑1 (MQ=255)
                                                                                                                          ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACa                         <  2:591076/150‑1 (MQ=255)
                                                                                                                           gCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGa                       <  2:1006527/151‑1 (MQ=255)
                                                                                                                            cGAAACGCGAATGAGTGGCCGTGGGGTGACTTCCTATC‑TTGAACGGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAGCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACat                        <  1:382568/149‑2 (MQ=255)
                                                                                                                                aCGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAACAGGCA‑CGCCCGAAATCTAATGCCCGCATGATCATCCTCGATCCGCGCTATACCCACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTTGTCCGGGAACAGATTcc                   >  2:488566/1‑150 (MQ=255)
                                                                                                                                aCGCGAATGAGTGGCGGTGGGGGGAATTACTATC‑TTGAACAGGCA‑CGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGACAGATGcc                   >  2:161954/1‑150 (MQ=255)
                                                                                                                                          gTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTaa        <  2:34488/151‑1 (MQ=255)
                                                                                                                                            ggCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg      <  2:914938/151‑1 (MQ=255)
                                                                                                                                                gTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt  <  1:995883/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCA‑CGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCT  >  NC_000913/941563‑941857

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