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

GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCT  >  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)
            ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATCt                                                                                                                                       >  2:996390/1‑151 (MQ=255)
             cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCtt                                                                                                                                      <  2:661609/151‑1 (MQ=255)
               cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTg                                                                                                                                     <  2:527887/150‑1 (MQ=255)
               cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGCCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGa                                                                                                                                    <  1:349482/151‑1 (MQ=255)
                 cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTATCTTTaa                                                                                                                                   >  2:507310/1‑150 (MQ=255)
                      tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGc                                                                                                                              <  2:738268/150‑1 (MQ=255)
                       gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCAc                                                                                                                           <  2:683265/151‑1 (MQ=255)
                        cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACg                                                                                                                          <  2:477707/151‑1 (MQ=255)
                            aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATTTTTaaaaaggccgcc                                                                                                                         >  1:810681/1‑139 (MQ=255)
                              gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaa                                                                                                                    <  1:670155/151‑1 (MQ=255)
                                   aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGGTTACTATCTTTGAAAAGGACGCCCGAAAACt                                                                                                                >  1:992448/1‑150 (MQ=255)
                                   aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATTTTGAAAAAGGACGCCAGaaa                                                                                                                   >  1:563618/1‑147 (MQ=255)
                                         aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGCTTACTATCTTTAAAAGGGACGCCCGAAATCTaaaagcc                                                                                                         >  2:930978/1‑146 (MQ=255)
                                          ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTATCTTTAAAAGGCACGCCCGAAAACTaaaaccc                                                                                                         >  1:1205553/1‑145 (MQ=255)
                                            tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGCTATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGGAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTCCATCTTGTAAAAGCACGGCACAAAAATAATGcc                                                                                                          >  1:479519/1‑147 (MQ=255)
                                            tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGACCTACTATTTTTAAAAGGGACCCCCAAAAATTTAAtcacgcc                                                                                                       >  2:668401/1‑145 (MQ=255)
                                             aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                     <  1:194183/151‑1 (MQ=255)
                                              cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATg                                                                                                    <  1:22016/151‑1 (MQ=255)
                                               ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGTCCTACTAACTTTAAAAGGCACCCCACAAAAATAATGTCCGCAt                                                                                                     >  1:519234/1‑149 (MQ=255)
                                               ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAAGGACGCCCGAAATCTAAAGCCCGCCTg                                                                                                    >  2:494676/1‑150 (MQ=255)
                                               ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATCTTGAAAACGCACGCCAGAAAATTAATGTCCCCATg                                                                                                    >  2:626185/1‑150 (MQ=255)
                                                      gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcat                                                                                            <  2:726873/151‑1 (MQ=255)
                                                      gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatca                                                                                             <  2:476335/150‑1 (MQ=255)
                                                      gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatca                                                                                             <  2:952560/150‑1 (MQ=255)
                                                      gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatca                                                                                             <  2:50982/150‑1 (MQ=255)
                                                        gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                          <  2:466087/151‑1 (MQ=255)
                                                        gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAAGGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                          <  1:893896/151‑1 (MQ=255)
                                                              ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACGTACTATTTTTAACAAGGACGCCCGAAAACTAATGGCCGccaggtcattcatga                                                                                         >  1:751099/1‑132 (MQ=255)
                                                                 aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCt                                                                                 <  1:53264/151‑1 (MQ=255)
                                                                    aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCtata                                                                              <  1:369736/151‑1 (MQ=255)
                                                                      cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTACTATCTTTAAAAAGCACCCCAGAAAACTAATACCCGCATGGTCATCatccatcccccccat                                                                               >  1:202168/1‑141 (MQ=255)
                                                                       ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTGAACAAGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCCATAcc                                                                            >  1:34488/1‑150 (MQ=255)
                                                                                atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACGCCGGt                                                                   <  1:206033/150‑1 (MQ=255)
                                                                                  cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGc                                                                 <  2:341385/150‑1 (MQ=255)
                                                                                       aCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACCAGCACGCCCGAAAACTAATGCCCGCACGATCCTCATCGATCCGCGCTTTACCCACAACCGTGCCggg                                                             >  2:878650/1‑149 (MQ=255)
                                                                                        cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcgc                                                          <  1:19803/151‑1 (MQ=255)
                                                                                            aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAg                                                      <  1:517336/151‑1 (MQ=255)
                                                                                            aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGCTGGGGTGACTTACTATCTTGAACAGGCACGCCAAAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAg                                                      <  1:207719/151‑1 (MQ=255)
                                                                                             aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGa                                                     <  2:921715/151‑1 (MQ=255)
                                                                                                tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGCATACTATCTTGAAAAAGCGCCCCCTAAAATTAAAGGCCGCCTCAATAACATCATTCCCCCCTATTCCGCCACCCGCGGCGGGCGGGGAGAt                                                    >  1:872630/1‑149 (MQ=255)
                                                                                                 ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATg                                                   <  1:1114970/149‑1 (MQ=255)
                                                                                                   tAGTGCTGTTTGGTCATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAGCTTGAACAGGCGCGCCAGAAATCTAATGCCCGCATGATCATCATCGACCCGCGCTATATCGACACCGGTGCCGGTCGCGAAGATgt                                                  <  1:198962/148‑2 (MQ=255)
                                                                                                    aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATTTTTAAAAGGCACCCCAGAAAACTAATGCCCCCCTGATCATCATCGATCCGCGCCATACCCACACCCGGGCCGGGCGCGGAGATGAGTgg                                              >  2:1001614/1‑151 (MQ=255)
                                                                                                            ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCCCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGGTCCCTAt                                      >  2:1079251/1‑151 (MQ=255)
                                                                                                            ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGCTTAACATCTTTAACAAGCAACCCAGGAATTTAATGACCGCATGATCATCATCGATTCGCGCTTTACCCCCACCGGGGGCGGGCGCGAAGATGAGGGGaaacc                                         >  1:1044527/1‑144 (MQ=255)
                                                                                                                taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTTTCTTTAACAGGCACGCCAGAAAACTAATGCCCGCCTGAACATCATCGATCCGCGCTATACCGACACCCGTGCCGGGCGCGGAGATGAGGGGGTCCCTAtt                                     >  2:1173760/1‑148 (MQ=255)
                                                                                                                taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAAGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGAACCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGGGGATCCCCAtt                                     >  2:703929/1‑148 (MQ=255)
                                                                                                                taataaCCCTGGCGAAACGCGAATGAGTGGCGGTCGGGGGGATTAATATCTTGAAAAGGCACGCCAGAAAACTAATGCCCGCCTGATCATCTTCCGTCCGCGCCATAACGACACCGGGGCCGGGGGCGAACATGAGTGGGTcc                                          >  1:718815/1‑143 (MQ=255)
                                                                                                                    aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTTAAAACGCACGCCAGAAATCTAATGCCCGCATGATCATCAATGAACCGCGCCATACCCACACCCGGGCCGGGCGGGAAGATGAGGGGATCCCTATTCGTCCg                               >  2:532170/1‑150 (MQ=255)
                                                                                                                     aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                             <  1:996337/151‑1 (MQ=255)
                                                                                                                      cccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGa                            <  2:85972/151‑1 (MQ=255)
                                                                                                                         tGGCGAAACGCGAATGAGTGGCGGAGGGGGTAGTTAATATCTTTAACAGGGACGCCCGAAATCTAATGCCCCCACGGTCCTCAACCATCCTCCCTCTACCCGCACCGGTGACCGGCGCGGAGATGAGGGGGGCCCTATTCCTCCgcggac                          >  1:1179668/1‑145 (MQ=255)
                                                                                                                          ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACa                         <  1:1025335/150‑1 (MQ=255)
                                                                                                                          ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACa                         <  2:591076/150‑1 (MQ=255)
                                                                                                                           gCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGa                       <  2:1006527/151‑1 (MQ=255)
                                                                                                                            cGAAACGCGAATGAGTGGCCGTGGGGTGACTTCCTATCTTGAACGGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAGCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACat                        <  1:382568/149‑2 (MQ=255)
                                                                                                                                aCGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCCTCGATCCGCGCTATACCCACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTTGTCCGGGAACAGATTcc                   >  2:488566/1‑150 (MQ=255)
                                                                                                                                aCGCGAATGAGTGGCGGTGGGGGGAATTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGACAGATGcc                   >  2:161954/1‑150 (MQ=255)
                                                                                                                                          gTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTaa        <  2:34488/151‑1 (MQ=255)
                                                                                                                                            ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg      <  2:914938/151‑1 (MQ=255)
                                                                                                                                                gTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt  <  1:995883/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCT  >  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: