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

CTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTG  >  NC_000913/941567‑941858
                                                                                                                                                  |                                                                                                                                                   
ctactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAAAGAGTGGCGGGTGGGTGACt                                                                                                                                                 >  2:276431/1‑151 (MQ=255)
        ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCt                                                                                                                                        <  2:514066/151‑1 (MQ=255)
                    cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGAGTTACTATCTTGAAAAGGCAc                                                                                                                             >  2:731379/1‑150 (MQ=255)
                     ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGCTACTATCTTTAACAAGCAc                                                                                                                             >  1:687673/1‑149 (MQ=255)
                          gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaa                                                                                                                      <  2:1123346/151‑1 (MQ=255)
                           tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACGCCCGaa                                                                                                                      >  1:851767/1‑150 (MQ=255)
                              gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTTAACAAGGACCCCCGAAAACt                                                                                                                  >  2:819946/1‑151 (MQ=255)
                               aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTa                                                                                                                 <  2:416746/151‑1 (MQ=255)
                                aCTATACATACGGCGGCTGGGCAGATGGCAACAGCCCGTCCGATGTCGAAAACCGTAAACTGGTAGTGCTGGTTGGTTATAACCCGGGCGAAACGCGAATGAGAGGCGGGGGGGTGTCTTTATCTCTTGAAAAggac                                                                                                                              >  2:646905/1‑135 (MQ=255)
                                   ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATg                                                                                                              <  1:411466/150‑1 (MQ=255)
                                    taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGTCTTACTATCTTTAAAAGGCCCCCCCGAAATTTAATGcc                                                                                                            >  1:250850/1‑151 (MQ=255)
                                    taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTGAACAGGCCCGCCCGAAATCTAATGcc                                                                                                            >  1:455221/1‑151 (MQ=255)
                                         aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGc                                                                                                         <  1:642456/149‑1 (MQ=255)
                                             cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatc                                                                                                   <  1:86581/151‑1 (MQ=255)
                                             cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAGGCACGCCAGAAAAATAATGCCCCCATGatc                                                                                                   >  1:3202/1‑151 (MQ=255)
                                               gcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGACATACTATCTTTTACAAGGGCGCCAGAAATATAATGccc                                                                                                           >  1:371969/1‑141 (MQ=255)
                                                     cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGa                                                                                           <  2:644383/151‑1 (MQ=255)
                                                         tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAATCTAATGCCCGCAAGATCATCATCGAttc                                                                                        >  1:54161/1‑148 (MQ=255)
                                                              aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCta                                                                                  <  2:687712/151‑1 (MQ=255)
                                                               cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAAACTAAAGCCCGCATGATCATCATCGATCcg                                                                                       >  1:637255/1‑145 (MQ=255)
                                                                 gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTTCTATCTTTAACCGGCACCCCAGAAATCTAATGCCCGCATGGTCATCATCGATCCGCGCtata                                                                                >  2:139093/1‑150 (MQ=255)
                                                                   ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAGGCACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCTATAccc                                                                             >  2:333944/1‑150 (MQ=255)
                                                                   ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCAGAAAATTAATGCCCCCACGATCATCATCGATCCCCCCtataac                                                                              >  2:728630/1‑148 (MQ=255)
                                                                   ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACCTATTATCTTGAACAGGGACGCCAGAAATCTAAAGCCCGCCCGATCATCCTCGATCCGCGCtat                                                                                 >  2:29325/1‑147 (MQ=255)
                                                                   ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTATTATCTTGAAAAAGCACGCCCGAAATCTAATGTCCCCATGGTCATCATCGATCCGCGCTATAcc                                                                              >  2:1071573/1‑150 (MQ=255)
                                                                     gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATCTTGAACAAGCACGCCCGAAAAATAATGCCCCCCAGATCATTATCGGTCCCCGCcataac                                                                              >  1:50286/1‑146 (MQ=255)
                                                                       cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGAACATCATAGATCCGCCCTATAACGacac                                                                        >  1:501756/1‑151 (MQ=255)
                                                                           tatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGGAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACAcccg                                                                     >  1:380431/1‑148 (MQ=255)
                                                                            atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAGGCACGCCAGAAAACTAATGCCCCCCTGATCATCATCGATCCGCGCCATACCCACAcccgcg                                                                   >  1:408604/1‑147 (MQ=255)
                                                                                aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAAACTAATACCCGCACGAACCTCATCGACTCCCGCCATACCCACACCCGTGcccg                                                               >  2:618319/1‑149 (MQ=255)
                                                                                aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGTGCCgg                                                               >  1:85404/1‑151 (MQ=255)
                                                                                  aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTCTATCTTGAAAAGGGACGCCCGAAAACTAAAGCCCGCCTGAACATCATCGATCCCCGCCATACCGACACCCGGGCCcgg                                                              >  2:1125992/1‑150 (MQ=255)
                                                                                          gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                     <  2:961102/151‑1 (MQ=255)
                                                                                            tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTATTAACTTGAAAAGGCACGCCAGAAATCTAAAGCCCGCCTGATCATCATCCGTCCGCGCCCAACCGACACCGGTGCCGGGcgcg                                                          >  2:166718/1‑144 (MQ=255)
                                                                                             ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGa                                                   <  2:236697/150‑1 (MQ=255)
                                                                                             ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATATTGAAAAAGCACGCCAGAAAACTAAAGCCCGCCTGATACTCATCGAACCGCGCTATACCGACACCCGGGCCCGGCGCGAAAATg                                                    >  2:949375/1‑149 (MQ=255)
                                                                                              gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGt                                                 <  2:1099324/151‑1 (MQ=255)
                                                                                                   gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGAt                                             <  2:371958/150‑1 (MQ=255)
                                                                                                   gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATc                                            <  1:618335/151‑1 (MQ=255)
                                                                                                      gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATccc                                          <  1:678762/150‑1 (MQ=255)
                                                                                                      gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGGACGCCAGAAAACTAAAGCCCCCATGATCATCATCCATCCGCGCCATACCCACACCCGTGCCCGGCGCGAAGATGAGTGGATccc                                          >  1:405337/1‑150 (MQ=255)
                                                                                                       tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAAGCACGCCAGAAATCTAATGCCCGCATGATCCTCATCGATCCCCGCTATACCGAAACCGGTGCCGGGGGCGAAGATGAGGGGGTccc                                          >  1:65312/1‑149 (MQ=255)
                                                                                                        ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTa                                        <  2:1021993/150‑1 (MQ=255)
                                                                                                          ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                     <  2:118458/151‑1 (MQ=255)
                                                                                                             aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCTGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGt                                   <  2:637232/150‑1 (MQ=255)
                                                                                                              ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACTGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTcc                                 <  2:526893/151‑1 (MQ=255)
                                                                                                                aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCgg                               <  1:1107919/151‑1 (MQ=255)
                                                                                                                aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGAAGATGAGTGGATCCCTATTCGTCCgg                               >  1:982298/1‑151 (MQ=255)
                                                                                                                     tGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCAATCTTGAACACGCACGCCCGCAAACTAATGCCCCCATGATCATCATCGATTCGCGCTATAACCACACCCGGGCCCGGCGCGAAGATGAGTGGATCCCTATTCGCCCgg                               >  2:292786/1‑146 (MQ=255)
                                                                                                                      ggCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAAAAAACTAATGCCCCCATGATTATCATCGGTCCGCGCTATACCCACACCGGTGCCGGGCGCGCAGATGAGTGGATCCCTATTCTTCCgggga                            >  1:8789/1‑146 (MQ=255)
                                                                                                                       gCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAACAGGACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCTATACCCACACCCGCGCCGGGCGCGAAGATGAGTGGGTCCCCATTCGTCCGGGaa                            >  1:248956/1‑147 (MQ=255)
                                                                                                                        cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                       <  1:859790/151‑1 (MQ=255)
                                                                                                                        cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                       <  2:636466/151‑1 (MQ=255)
                                                                                                                         gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATg                      <  2:920822/151‑1 (MQ=255)
                                                                                                                          aaaCGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTTAAAAGGCCCGCCAGAAATCTAATGCCCCCATGATCCTCATCGGTTCGCGCTCTACCGCCACCCGGGGCCGGCGGGGAAGAGAGTGGGTCCCCATTTGGCCGGGGACAGaag                      >  2:769722/1‑148 (MQ=255)
                                                                                                                          aaaCGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGc                     >  2:183605/1‑151 (MQ=255)
                                                                                                                           aaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGcc                    <  1:789163/151‑1 (MQ=255)
                                                                                                                           aaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGcc                    <  2:281776/151‑1 (MQ=255)
                                                                                                                           aaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGcc                    <  2:782729/151‑1 (MQ=255)
                                                                                                                            aCGCGAATGAGTGGCGGTGGGGTGACGTACTATCCTGAACAGGCACGCCAGAAAACTAATGGCCGCAAGATTATCAACGATCCGCGCCATAACCACACCCGTGCCGGGCGCGCAGATGAGGGGGTCCCTATTCGTCCGGGGACAGATgtc                    >  2:609912/1‑148 (MQ=255)
                                                                                                                               cgAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGGACGCCAGAAATCTAATGCCCGCATGGTCATCATCGCTCCGCGCTATACCGACACCCGGGCCGGGCGCGAAGATGAATGGGACCCCATTCGTTCGGGAACAGATGCCGc                  >  2:642215/1‑149 (MQ=255)
                                                                                                                                 aaTGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTg              <  1:819970/151‑1 (MQ=255)
                                                                                                                                       tGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc        >  1:216396/1‑151 (MQ=255)
                                                                                                                                        ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc        <  1:120488/150‑1 (MQ=255)
                                                                                                                                          cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt     <  2:182258/151‑1 (MQ=255)
                                                                                                                                             tGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg  >  2:900170/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTG  >  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: