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

ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCT  >  NC_000913/941569‑941857
                                                                                                                                                 |                                                                                                                                                  
actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACt                                                                                                                                                <  1:952580/149‑1 (MQ=255)
      ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAt                                                                                                                                          <  1:1527708/149‑1 (MQ=255)
      ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑Ct                                                                                                                                       <  2:394691/151‑1 (MQ=255)
         cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGATTTCTAT‑Ctt                                                                                                                                      >  2:1369816/1‑149 (MQ=255)
         cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑GGGGGGGGCTTTCTAT‑Ctt                                                                                                                                      >  1:94424/1‑149 (MQ=255)
          gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑GGGGGGGGCTTACTAt‑ttt                                                                                                                                      >  2:354929/1‑148 (MQ=255)
           cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTTCTAT‑CTTTaa                                                                                                                                   >  2:1430227/1‑150 (MQ=255)
              aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTTCTAT‑CTTTaaaa                                                                                                                                 >  2:602754/1‑147 (MQ=255)
                 gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCa                                                                                                                             <  1:340888/150‑1 (MQ=255)
                   ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGc                                                                                                                          <  1:1197528/151‑1 (MQ=255)
                   ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTAATAT‑CTTTaaaaag                                                                                                                               >  2:743040/1‑142 (MQ=255)
                    gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑GGGGGGGGCGTTCTAT‑CTTGAAAAAGGACGc                                                                                                                          >  1:1418827/1‑150 (MQ=255)
                     aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCa                                                                                                                        <  2:825777/151‑1 (MQ=255)
                     aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCa                                                                                                                        <  1:274460/151‑1 (MQ=255)
                     aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑GGGGGGGACTTACTAT‑CTTTAACAGGGACGcc                                                                                                                         >  2:259058/1‑150 (MQ=255)
                           ttaaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACCTTATTT‑CTTTAACCAGCACGCCCGaa                                                                                                                     >  2:214027/3‑148 (MQ=255)
                           tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTTCTAT‑CTTTAAAAAGCACGCCCGaaa                                                                                                                    >  2:1249515/1‑149 (MQ=255)
                               cTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGGTTTCTAT‑CTTTAAAAAGCACGCCCGAAAACTaa                                                                                                               >  1:486332/1‑150 (MQ=255)
                                 ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAAt                                                                                                              <  1:763570/149‑1 (MQ=255)
                                  taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑GGGGGGGGGTTTTTAT‑CTTTTACAAGCACGCCACAAAACTAATgg                                                                                                            >  1:400882/1‑149 (MQ=255)
                                  taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑GGGGGGGGGTTTCTAT‑CTTTAAAAAGGACGCCCGAAATCTaaatgc                                                                                                           >  1:1267368/1‑147 (MQ=255)
                                   aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGGTTACTTT‑TTTTAACAACGGCGCCAGAAAACTTAAGcc                                                                                                           >  1:1585909/1‑150 (MQ=255)
                                       aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCa                                                                                                       <  2:951278/150‑1 (MQ=255)
                                       aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCa                                                                                                       <  2:1331352/150‑1 (MQ=255)
                                       aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                      <  1:1467394/151‑1 (MQ=255)
                                       aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                      <  2:889770/151‑1 (MQ=255)
                                         ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTACTAT‑TTTGAACAGGCACGCCCGAAATCTAATGCCCGCCTg                                                                                                     >  2:906047/1‑150 (MQ=255)
                                           cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TAGGGTGGCTTATTTT‑CTTGAAACGGCACGCCCGAAATCTAATGcaccca                                                                                                       >  1:1097823/1‑141 (MQ=255)
                                            ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatca                                                                                                 <  1:437795/151‑1 (MQ=255)
                                                gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcat                                                                                             <  2:1247868/151‑1 (MQ=255)
                                                gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatca                                                                                              <  2:1106747/150‑1 (MQ=255)
                                                   cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTTTTAT‑CTTCAACAGGCCCGCCCCAAAACTAATGCCCGCCTGatcatcttc                                                                                            >  2:1403128/1‑149 (MQ=255)
                                                   cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑GGGGGGGTCGTAATATATTTGAACA‑GCCCGGCCGAAATTCAATGCCCGCTCCatcatcatc                                                                                            >  1:1402744/1‑149 (MQ=255)
                                                             cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCtat                                                                                <  2:529521/151‑1 (MQ=255)
                                                              aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCtata                                                                               <  1:235328/151‑1 (MQ=255)
                                                               gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAc                                                                              <  1:259134/151‑1 (MQ=255)
                                                                   gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGTCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGac                                                                          <  2:1510553/151‑1 (MQ=255)
                                                                     cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGaca                                                                         <  2:104310/150‑1 (MQ=255)
                                                                     cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTACTAT‑TTTGAACAGGGACGCCAGAAATCTAATGCCCGCCTGAATATCATCGATCCGCGCTATACCGaca‑c                                                                       >  2:1082347/1‑151 (MQ=255)
                                                                     cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTACTAT‑CTTTAACAGGCACCCCCGAAAACTAATGCCCGCATGATCATCATCGATCCCCGCTATACCCaca‑c                                                                       >  1:277963/1‑151 (MQ=255)
                                                                      ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑cc                                                                      <  2:146865/151‑1 (MQ=255)
                                                                            cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGc                                                                 <  1:1449724/150‑1 (MQ=255)
                                                                              aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCgg                                                              <  1:1645177/151‑1 (MQ=255)
                                                                              aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTACTAT‑CTTGAACAGGCACGCCCGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCCACA‑CCCGGGCCgg                                                              >  2:1495131/1‑151 (MQ=255)
                                                                               aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTACTAT‑CTTGAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGGTCCCCGCTATACCGCCA‑Cccgggccggg                                                             >  2:408904/1‑151 (MQ=255)
                                                                                  cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGgcgc                                                          <  2:34555/151‑1 (MQ=255)
                                                                                    gTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGAATTACTAT‑TTTTAACAAGCACGCCAGAAAACTAATGCCCCCAAGATCATCATCGATCCCCGCCATAACGGCA‑CCGGGGGCGGgcgc                                                          >  2:441453/1‑149 (MQ=255)
                                                                                      aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGaa                                                       <  2:758498/150‑1 (MQ=255)
                                                                                            gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTACTAT‑CTTTAAAAGGCACGCCCGAAATCTAATGGCCGCAAGATCATCATCGATCCGCGCTCTACCCACA‑CCCGTGCCGGGCGCGAAGATGAg                                                 >  1:1372460/1‑150 (MQ=255)
                                                                                            gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTTCTAT‑CTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCCATACCGACA‑CCGGTGCCGGGggcg                                                         >  2:640263/1‑142 (MQ=255)
                                                                                              aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTg                                               <  2:1262806/150‑1 (MQ=255)
                                                                                                 gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGAt                                            <  2:473461/150‑1 (MQ=255)
                                                                                                   tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑GGGGGGGGATTACTAT‑CTTTAAAAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGGTCCGGCCCATACCGGCACCCGGGGCCGGG‑GCGGAGATGAGGGGGTccc                                         >  2:954647/1‑151 (MQ=255)
                                                                                                    gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                        <  2:1090408/151‑1 (MQ=255)
                                                                                                       tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt                                      <  2:1485039/150‑1 (MQ=255)
                                                                                                        ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                    <  2:1374198/151‑1 (MQ=255)
                                                                                                         gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGATGGGGGGA‑TTAACAA‑CTTTAACAGGGCAGCCACAAAAATAATGGGCGCAGGGGCCTCATCGAACCCCGTCTTAACGAAC‑CCCGGGGCGGGGGCGAAGATCAGTgggacc                                          >  1:1077376/1‑140 (MQ=255)
                                                                                                          taataaCCCTGGCTAAACGCGAATGAGTGTCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGGCA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCg                                   <  1:397530/150‑1 (MQ=255)
                                                                                                              aaCCCTGGCGAAACGCGAATGAGTGGCGG‑AGGGGTGGCTTACTTT‑ATTGAAAACGCCCGCCCGAAAATTAAAGCCCGCCTGAGAATCCACGGTGCGCGCCATATCCACC‑CCCGGGGCGGGGGCGAAGATGAGCGGAGGCCCATTCGTCCgg                              >  2:1673644/1‑151 (MQ=255)
                                                                                                                cccTGGCGAAACGCGAATGAGTGGCGG‑TTGGGGGTCCTACTTT‑CTTTAACAGGCAAGCCAGAAAACTAATGCCCGCCTGATCATCATCGATCCGCGCCATAACAACA‑CCGGGGCCGGGCGCGGAGATGAGTGGGTCCCCATTCGTCCggg                             >  1:946070/1‑150 (MQ=255)
                                                                                                                 ccTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTACTAT‑CTTGAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGGGCCGGGCGCGAAGATGAGTGGATCCCCCTTCGTCCggg                             >  2:1436675/1‑149 (MQ=255)
                                                                                                                 ccTGGCGAAACGCGAATGAGTGGCGG‑CGGGGGGATTTATTAT‑CCTCAACAGGCACGCCAGAAATCTAATGCCCGCCCGATTATCATCGCACCGCGCTATACCCCCA‑CCGCGGCCCGGCGCGAAGATGAGTGGGTCCCTATTCGTCCggg                             >  1:8387/1‑149 (MQ=255)
                                                                                                                  cTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGTTATTACTCT‑CTTTAACAAGCACGCCAGCAAACTAATGGCCGCCACATCCTCATCCATTCCCGCCATACCGACC‑ACCGCGCCGGGGGCCAAGATTGGGGGATCCCTATTCGTCCgggga                           >  2:1046225/1‑148 (MQ=255)
                                                                                                                     gCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTACTAT‑CTTGAAAAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACA‑CCGGGGCCGGGGGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc                          >  2:619932/1‑148 (MQ=255)
                                                                                                                      cGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                      <  1:1354010/151‑1 (MQ=255)
                                                                                                                       gAAACGCGAATGAGTGGCGG‑TGGGGGGACTTACTAT‑CTTGAAAAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGGGCCGGGCGCGAAGATGAGTGGGTCCCTATTCGTCCGGGGACAGa                       >  2:239582/1‑149 (MQ=255)
                                                                                                                        aaaCGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGc                    <  2:1578904/151‑1 (MQ=255)
                                                                                                                         aaCGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGcc                   <  1:826266/151‑1 (MQ=255)
                                                                                                                           cgcgAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGc                 <  1:1249470/151‑1 (MQ=255)
                                                                                                                            gcgAATGAGTGGCGG‑TGGGGTGGCCTACTAT‑CTTGAAAAGGCACGCCCGAAAGCTCATGCCCGCGTGATCATAATCGATCCGCGCCATACCGACA‑ACGGTGGCGGGGGCGGAGAGGAGGGGGTCCCCCTTCGTCCGGGGACAGATGCCg                  >  1:59497/1‑149 (MQ=255)
                                                                                                                            gcgAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCa                <  1:483950/151‑1 (MQ=255)
                                                                                                                            gcgAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCa                <  1:1006136/151‑1 (MQ=255)
                                                                                                                            gcgAATGAGTGGCGG‑TGGGGGGACTTACTAT‑CTTGAACAAGGACGCCCGAAATATAATGCCCGCATGATCATCATCGATCCGCGCTATAACGGCA‑CCGGTGCCGGGGGCGAAGATGAGTGGGTCCCCATTCGTCCGGGGACAGATGcc                   >  2:94703/1‑148 (MQ=255)
                                                                                                                             cgAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCCTCATCGATCCGCGCTATACCGACA‑CCGGGGCCCGGCGCGGAGATGAGGGGGTCCCCCTTCGTCCCGGGACAGATGCCGc                 >  1:1427784/1‑149 (MQ=255)
                                                                                                                              gAATGAGTGGCGG‑TGGGGGGGCTTACTAT‑CTTGAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAACCACA‑CCGGTGCCGGGCGCGAAGATGAGTGGGTCCCTATTTGTCCGGGAACAGATTCCGCa                >  2:594166/1‑149 (MQ=255)
                                                                                                                                aTGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTgg            <  2:330669/151‑1 (MQ=255)
                                                                                                                                 tGAGTGGCGG‑TGGGGTGGATTTCTAT‑CTTGAACAAGCACGCCCGAAAATTAATGCCCGCATGCTCATCCTCCATCCGCGCGATACCCCCA‑CCCGTGCCGGGCGCGACGATGAGGGGGTCCCCATTCGTCCGGGAACAGATGCCGCCCTg             >  1:1064129/1‑149 (MQ=255)
                                                                                                                                      ggCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg      >  1:999432/1‑151 (MQ=255)
                                                                                                                                          g‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt  <  1:471608/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTAT‑CTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCT  >  NC_000913/941569‑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: ‑