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

CTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTG  >  NC_000913/941567‑941858
                                                                                                                                                  |                                                                                                                                                   
ctactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAAAGAGTGGCGGGTGGGTGACt                                                                                                                                                 >  2:276431/1‑151 (MQ=255)
        ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑t                                                                                                                                        <  2:514066/151‑1 (MQ=255)
                    cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGAGTTACTATC‑TTGAAAAGGCAc                                                                                                                             >  2:731379/1‑150 (MQ=255)
                     ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGCTACTATC‑TTTAACAAGCAc                                                                                                                             >  1:687673/1‑149 (MQ=255)
                          gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGaa                                                                                                                      <  2:1123346/151‑1 (MQ=255)
                           tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTGAACAAGCACGCCCGaa                                                                                                                      >  1:851767/1‑150 (MQ=255)
                              gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATC‑TTTAACAAGGACCCCCGAAAACt                                                                                                                  >  2:819946/1‑151 (MQ=255)
                               aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTa                                                                                                                 <  2:416746/151‑1 (MQ=255)
                                aCTATACATACGGCGGCTGGGCAGATGGCAACAGCCCGTCCGATGTCGAAAACCGTAAACTGGTAGTGCTGGTTGGTTATAACCCGGGCGAAACGCGAATGAGAGGCGGGGGGGTGTCTTTATCTC‑TTGAAAAggac                                                                                                                              >  2:646905/1‑135 (MQ=255)
                                   ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATg                                                                                                              <  1:411466/150‑1 (MQ=255)
                                    taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGTCTTACTATC‑TTTAAAAGGCCCCCCCGAAATTTAATGcc                                                                                                            >  1:250850/1‑151 (MQ=255)
                                    taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATC‑TTGAACAGGCCCGCCCGAAATCTAATGcc                                                                                                            >  1:455221/1‑151 (MQ=255)
                                         aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGc                                                                                                         <  1:642456/149‑1 (MQ=255)
                                             cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatc                                                                                                   <  1:86581/151‑1 (MQ=255)
                                             cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAAAAGGCACGCCAGAAAAATAATGCCCCCATGatc                                                                                                   >  1:3202/1‑151 (MQ=255)
                                               gcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGACATACTATC‑TTTTACAAGGGCGCCAGAAATATAATGccc                                                                                                           >  1:371969/1‑141 (MQ=255)
                                                     cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGa                                                                                           <  2:644383/151‑1 (MQ=255)
                                                         tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAACAAGCACGCCAGAAATCTAATGCCCGCAAGATCATCATCGAttc                                                                                        >  1:54161/1‑148 (MQ=255)
                                                              aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCta                                                                                  <  2:687712/151‑1 (MQ=255)
                                                               cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTGAACAGGCACGCCAGAAAACTAAAGCCCGCATGATCATCATCGATCcg                                                                                       >  1:637255/1‑145 (MQ=255)
                                                                 gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTTCTATC‑TTTAACCGGCACCCCAGAAATCTAATGCCCGCATGGTCATCATCGATCCGCGCtata                                                                                >  2:139093/1‑150 (MQ=255)
                                                                   ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAAAAGGCACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCTATAccc                                                                             >  2:333944/1‑150 (MQ=255)
                                                                   ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTGAACAAGCACGCCAGAAAATTAATGCCCCCACGATCATCATCGATCCCCCCtataac                                                                              >  2:728630/1‑148 (MQ=255)
                                                                   ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACCTATTATC‑TTGAACAGGGACGCCAGAAATCTAAAGCCCGCCCGATCATCCTCGATCCGCGCtat                                                                                 >  2:29325/1‑147 (MQ=255)
                                                                   ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTATTATC‑TTGAAAAAGCACGCCCGAAATCTAATGTCCCCATGGTCATCATCGATCCGCGCTATAcc                                                                              >  2:1071573/1‑150 (MQ=255)
                                                                     gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATC‑TTGAACAAGCACGCCCGAAAAATAATGCCCCCCAGATCATTATCGGTCCCCGCcataac                                                                              >  1:50286/1‑146 (MQ=255)
                                                                       cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATT‑TTTAACAGGCACGCCAGAAATCTAATGCCCGCATGAACATCATAGATCCGCCCTATAACGaca‑c                                                                        >  1:501756/1‑151 (MQ=255)
                                                                           tatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAACAGGCACGCCAGGAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑cccg                                                                     >  1:380431/1‑148 (MQ=255)
                                                                            atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAAAAGGCACGCCAGAAAACTAATGCCCCCCTGATCATCATCGATCCGCGCCATACCCACA‑cccgcg                                                                   >  1:408604/1‑147 (MQ=255)
                                                                                aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAACAGGCACGCCAGAAAACTAATACCCGCACGAACCTCATCGACTCCCGCCATACCCACA‑CCCGTGcccg                                                               >  2:618319/1‑149 (MQ=255)
                                                                                aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTGAACAGGCACGCCCGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCgg                                                               >  1:85404/1‑151 (MQ=255)
                                                                                  aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTCTATC‑TTGAAAAGGGACGCCCGAAAACTAAAGCCCGCCTGAACATCATCGATCCCCGCCATACCGACA‑CCCGGGCCcgg                                                              >  2:1125992/1‑150 (MQ=255)
                                                                                          gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGAt                                                     <  2:961102/151‑1 (MQ=255)
                                                                                            tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTATTAAC‑TTGAAAAGGCACGCCAGAAATCTAAAGCCCGCCTGATCATCATCCGTCCGCGCCCAACCGACA‑CCGGTGCCGGGcgcg                                                          >  2:166718/1‑144 (MQ=255)
                                                                                             ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGa                                                   <  2:236697/150‑1 (MQ=255)
                                                                                             ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATA‑TTGAAAAAGCACGCCAGAAAACTAAAGCCCGCCTGATACTCATCGAACCGCGCTATACCGACA‑CCCGGGCCCGGCGCGAAAATg                                                    >  2:949375/1‑149 (MQ=255)
                                                                                              gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGt                                                 <  2:1099324/151‑1 (MQ=255)
                                                                                                   gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGAt                                             <  2:371958/150‑1 (MQ=255)
                                                                                                   gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATc                                            <  1:618335/151‑1 (MQ=255)
                                                                                                      gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATccc                                          <  1:678762/150‑1 (MQ=255)
                                                                                                      gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTGAACAAGGACGCCAGAAAACTAAAGCCCCCATGATCATCATCCATCCGCGCCATACCCACA‑CCCGTGCCCGGCGCGAAGATGAGTGGATccc                                          >  1:405337/1‑150 (MQ=255)
                                                                                                       tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTGAAAAAGCACGCCAGAAATCTAATGCCCGCATGATCCTCATCGATCCCCGCTATACCGAAA‑CCGGTGCCGGGGGCGAAGATGAGGGGGTccc                                          >  1:65312/1‑149 (MQ=255)
                                                                                                        ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTa                                        <  2:1021993/150‑1 (MQ=255)
                                                                                                          ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                     <  2:118458/151‑1 (MQ=255)
                                                                                                             aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCTGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGt                                   <  2:637232/150‑1 (MQ=255)
                                                                                                              ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CTGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTcc                                 <  2:526893/151‑1 (MQ=255)
                                                                                                                aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCgg                               <  1:1107919/151‑1 (MQ=255)
                                                                                                                aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACA‑CCGGGGCCGGGGGCGAAGATGAGTGGATCCCTATTCGTCCgg                               >  1:982298/1‑151 (MQ=255)
                                                                                                                     tGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCAATC‑TTGAACACGCACGCCCGCAAACTAATGCCCCCATGATCATCATCGATTCGCGCTATAACCACA‑CCCGGGCCCGGCGCGAAGATGAGTGGATCCCTATTCGCCCgg                               >  2:292786/1‑146 (MQ=255)
                                                                                                                      ggCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTTAACAGGCACGCCAAAAAACTAATGCCCCCATGATTATCATCGGTCCGCGCTATACCCACA‑CCGGTGCCGGGCGCGCAGATGAGTGGATCCCTATTCTTCCgggga                            >  1:8789/1‑146 (MQ=255)
                                                                                                                       gCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAACAGG‑ACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCTATACCCACA‑CCCGCGCCGGGCGCGAAGATGAGTGGGTCCCCATTCGTCCGGGaa                            >  1:248956/1‑147 (MQ=255)
                                                                                                                        cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                       <  1:859790/151‑1 (MQ=255)
                                                                                                                        cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                       <  2:636466/151‑1 (MQ=255)
                                                                                                                         gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATg                      <  2:920822/151‑1 (MQ=255)
                                                                                                                          aaaCGCGAATGAGTGGCGGTGGGGGGACTTTCTATC‑TTTAAAAGGCCCGCCAGAAATCTAATGCCCCCATGATCCTCATCGGTTCGCGCTCTACCGCCACCCGGGGCCGGCGGGGAAGA‑GAGTGGGTCCCCATTTGGCCGGGGACAGaag                      >  2:769722/1‑148 (MQ=255)
                                                                                                                          aaaCGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGc                     >  2:183605/1‑151 (MQ=255)
                                                                                                                           aaCGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGcc                    <  1:789163/151‑1 (MQ=255)
                                                                                                                           aaCGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGcc                    <  2:281776/151‑1 (MQ=255)
                                                                                                                           aaCGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGcc                    <  2:782729/151‑1 (MQ=255)
                                                                                                                            aCGCGAATGAGTGGCGGTGGGGTGACGTACTATC‑CTGAACAGGCACGCCAGAAAACTAATGGCCGCAAGATTATCAACGATCCGCGCCATAACCACA‑CCCGTGCCGGGCGCGCAGATGAGGGGGTCCCTATTCGTCCGGGGACAGATgtc                    >  2:609912/1‑148 (MQ=255)
                                                                                                                               cgAATGAGTGGCGGTGGGGGGACTTACTATC‑TTGAACAGGGACGCCAGAAATCTAATGCCCGCATGGTCATCATCGCTCCGCGCTATACCGACA‑CCCGGGCCGGGCGCGAAGATGAATGGGACCCCATTCGTTCGGGAACAGATGCCGc                  >  2:642215/1‑149 (MQ=255)
                                                                                                                                 aaTGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTg              <  1:819970/151‑1 (MQ=255)
                                                                                                                                       tGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc        >  1:216396/1‑151 (MQ=255)
                                                                                                                                        ggCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc        <  1:120488/150‑1 (MQ=255)
                                                                                                                                          cGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt     <  2:182258/151‑1 (MQ=255)
                                                                                                                                             tGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg  >  2:900170/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTG  >  NC_000913/941567‑941858

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