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

Only 100 of 136 total aligned reads displayed.

GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941563‑941862
                                                                                                                                                      |                                                                                                                                                     
gCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                                                                                       <  2:783900/151‑1 (MQ=255)
  gACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGa                                                                                                                                                     <  2:567942/151‑1 (MQ=255)
     tactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                                  >  2:959654/1‑151 (MQ=255)
     tactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAAtt                                                                                                                                                  >  2:124316/1‑151 (MQ=255)
      actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTa                                                                                                                                                 <  1:1199136/151‑1 (MQ=255)
      actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                                  >  2:1042276/1‑150 (MQ=255)
       ctcctcCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGACtt                                                                                                                                                  >  1:1219818/1‑149 (MQ=255)
        tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTa                                                                                                                                                 >  1:959384/1‑149 (MQ=255)
             cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCtt                                                                                                                                          <  2:893121/151‑1 (MQ=255)
               cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTCTATCttt                                                                                                                                         >  2:1114298/1‑149 (MQ=255)
                      tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGc                                                                                                                                  <  1:163731/150‑1 (MQ=255)
                      tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGc                                                                                                                                  >  1:2154336/1‑150 (MQ=255)
                      tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGCTTTATTTCTTGaaaa                                                                                                                                     >  2:1749933/1‑145 (MQ=255)
                        cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTaaaaag                                                                                                                                   >  2:703750/1‑143 (MQ=255)
                         ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTaa                                                                                                                                       >  1:1129730/1‑142 (MQ=255)
                         ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTATTTTCTTTaaa                                                                                                                                      >  2:1883875/1‑142 (MQ=255)
                          gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGcc                                                                                                                             <  2:2040781/151‑1 (MQ=255)
                          gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTCTATCTTTAAAAAGCACGcc                                                                                                                             >  2:1645725/1‑151 (MQ=255)
                           aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCa                                                                                                                            <  2:365014/151‑1 (MQ=255)
                               tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTGAAAAGGCACGcccaaa                                                                                                                         >  1:2104413/1‑146 (MQ=255)
                               tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTTAAAAAGGACGcc                                                                                                                             >  2:2118452/1‑146 (MQ=255)
                                ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaaa                                                                                                                        <  1:572469/150‑1 (MQ=255)
                                 tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAAt                                                                                                                       <  2:310133/150‑1 (MQ=255)
                                 tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAAAAAGCACGCCCGaaa                                                                                                                        >  1:663775/1‑149 (MQ=255)
                                    aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTaa                                                                                                                   <  2:1184602/151‑1 (MQ=255)
                                    aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTaa                                                                                                                   <  1:1090377/151‑1 (MQ=255)
                                     cTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCAGGAATCTaa                                                                                                                   >  1:861153/1‑150 (MQ=255)
                                      tataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCAGAAATCTaa                                                                                                                   >  1:913061/1‑149 (MQ=255)
                                       ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTTAAAAAGCACGCCCGAAAACTaa                                                                                                                   >  1:280661/1‑148 (MQ=255)
                                             aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGc                                                                                                            <  2:1470595/149‑1 (MQ=255)
                                             aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCa                                                                                                           <  2:861142/150‑1 (MQ=255)
                                             aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTTTTTTTAACAAGCACCCCCGAAAACTAAAGccccc                                                                                                            >  2:112659/1‑147 (MQ=255)
                                              cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGc                                                                                                            <  2:1777916/148‑1 (MQ=255)
                                              cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATg                                                                                                         <  1:1180784/151‑1 (MQ=255)
                                                  ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTTTCTTTAACAGGCACGCCAGAAAACTAATGCCCGCATGa                                                                                                        >  2:1301894/1‑148 (MQ=255)
                                                  ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCATGatc                                                                                                      >  1:872473/1‑150 (MQ=255)
                                                      gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcat                                                                                                 <  2:72672/151‑1 (MQ=255)
                                                       ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAACAGGCACCCCCGAAAACTAAAGCCCGCATGattatc                                                                                                   >  1:1329458/1‑148 (MQ=255)
                                                           gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc                                                                                            <  2:1960815/151‑1 (MQ=255)
                                                           gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCCCATGATCATCATCGAt                                                                                             >  1:1700013/1‑150 (MQ=255)
                                                              ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcg                                                                                          <  1:653616/150‑1 (MQ=255)
                                                              ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcg                                                                                          <  1:1515039/150‑1 (MQ=255)
                                                              ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTTAAAAGGCACGCCAGAAATCTAATGCCCGCATGATCATCAACGATCcgc                                                                                         >  2:1239658/1‑151 (MQ=255)
                                                               gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgcg                                                                                        <  2:1388302/151‑1 (MQ=255)
                                                                         gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACCTACTATCTTTAAAAGGCACGCCAGAAATCTAAAGCCCGCATGATCATCATCGATCCGCGCTATACCGa                                                                               >  2:1713819/1‑150 (MQ=255)
                                                                         gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGGACGCCAGAAAACTAATGCCCGCATGATAATCATCGATCCGCGCTATAcc                                                                                 >  2:2106074/1‑148 (MQ=255)
                                                                         gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTGAAAAGGCACGCCCGAAAACTAATGCCCGCCTGATCCACATCGATCCgcgcctt                                                                                    >  2:425389/1‑142 (MQ=255)
                                                                             gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCg                                                                          <  1:1945950/151‑1 (MQ=255)
                                                                              atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGAGTTAATATCCTTAACAAGAACGCCAGAAATCTAATGCCCGCATAAACACCATCCATCCCCCCTATACCGacac                                                                            >  2:169041/1‑148 (MQ=255)
                                                                               tatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATCTTGAACAAGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCCCGCTATACCGACAcc                                                                           >  2:1795762/1‑148 (MQ=255)
                                                                                  cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTAAAATCTTTAAAAGGGAAGCCAGAAAACTAAAGCCCGCATGATCATCATTGATCCGCGCTATACCGACACCGGTg                                                                       >  2:1659319/1‑149 (MQ=255)
                                                                                   gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTACTATCTTGAACAGGCACGCCAGAAATCTAATTCCCGCATGATCATCATCGATCCGCGCTATACCCACACCCGGGcc                                                                     >  1:2016979/1‑150 (MQ=255)
                                                                                     aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCggg                                                                  <  2:809407/151‑1 (MQ=255)
                                                                                     aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAAAGGCACGCCAGAAAACTAAAGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCggg                                                                  >  1:72293/1‑151 (MQ=255)
                                                                                     aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCggg                                                                  >  1:1504405/1‑151 (MQ=255)
                                                                                      aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgc                                                                 <  2:609733/151‑1 (MQ=255)
                                                                                           tAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATTTTTAAAAGGCACGCCAGAAATCTAATGCCCGCCTGATTATCATTGATCCGCGCTATACCCACACCCGTGCCggg                                                                  >  2:267668/1‑145 (MQ=255)
                                                                                            aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGaa                                                            <  1:1867130/150‑1 (MQ=255)
                                                                                             aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTGAACAGGCACGCCCGAAAACTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGCCggg                                                                  >  2:1225163/1‑143 (MQ=255)
                                                                                               cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGa                                                          >  1:1249710/1‑149 (MQ=255)
                                                                                                   tAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGGAGAagaggg                                                    >  1:1906605/1‑149 (MQ=255)
                                                                                                    aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATCTTGAAAAGGCACGCCAGAAAACTAATACCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAgggg                                                   >  2:1702647/1‑151 (MQ=255)
                                                                                                     gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAAAAGGCACGCCCGAAAACTAAAGCCCGCCTGATAATCATTCATCCGCGCCATACCCACACCGGGGCCGGGGGCGAAGATGAGTgg                                                   >  2:363436/1‑150 (MQ=255)
                                                                                                      tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTTAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAgggg                                                   >  1:940938/1‑149 (MQ=255)
                                                                                                      tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGACCCGCGCTATACCCACACCGGGGCCGGGCGCGAAGATGAg                                                      >  1:281029/1‑146 (MQ=255)
                                                                                                       gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATc                                                >  1:877357/1‑151 (MQ=255)
                                                                                                       gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTTCTATCTTTAACAGGCACGCCAGAAAAATAATGCCCCCCTGATCCTCATCGATCCGCGCCTTACCGCCACCCGGGCCGGGGGCGAAAATGa                                                       >  2:717356/1‑144 (MQ=255)
                                                                                                       gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTACTATCTTGAACAAGCACGCCAGAAAAATAATGCCCGCATGATCATCATCGAACCGCGCTATACCCACACCGGGGCCGGGGGCGAAAAAGAgggg                                                   >  2:58419/1‑148 (MQ=255)
                                                                                                         tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGAAGATGAGTGGATccc                                              >  1:1753115/1‑151 (MQ=255)
                                                                                                          gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTGAAAAAGCCCGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATccc                                              >  1:527944/1‑150 (MQ=255)
                                                                                                            ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCGTACGATCTTTAACACGCGCGCCAAAAAACTAATGCCCGCAAGATCATCATCGATCCGCGCTATACCGACACCCGGGGCGGGCGCGAAGATGAGTGGATCCCTa                                            >  1:1974931/1‑150 (MQ=255)
                                                                                                            ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATTTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCCTCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGGTccc                                              >  1:1314133/1‑148 (MQ=255)
                                                                                                             tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt                                           <  1:1716340/150‑1 (MQ=255)
                                                                                                             tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCAGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt                                           <  1:790626/150‑1 (MQ=255)
                                                                                                              ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                         <  1:1771084/151‑1 (MQ=255)
                                                                                                               gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCg                                        <  2:1911504/151‑1 (MQ=255)
                                                                                                               gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCg                                        <  1:573380/151‑1 (MQ=255)
                                                                                                                taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAATCTAAAGACCGCATGATCCTCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGGAGATGAGTGGGTCCCTAtt                                          >  2:1990484/1‑148 (MQ=255)
                                                                                                                 aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGt                                       <  2:240815/150‑1 (MQ=255)
                                                                                                                 aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGt                                       <  1:1983331/150‑1 (MQ=255)
                                                                                                                 aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                                      <  2:447548/151‑1 (MQ=255)
                                                                                                                  ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGGACGCCCGAAATCTAAAGCCCGCATGATCATCATCGATCCGCGCTATAACCACACCGGTGCCGGGGGCGAAGATGAGAGGGTCCCTCTTCGTcc                                     >  2:1715883/1‑151 (MQ=255)
                                                                                                                   taaCCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATATAAAGCCCCCATGATCATCATCGATCCGCGCCATACCGACACCGGGGCCGGGGGCGAAGATGAGTGGATCCCTATTCGTcc                                     >  1:186118/1‑150 (MQ=255)
                                                                                                                    aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCCGAAATCTAATGCCCCCATGATCATCATCGGTCCGCGCTATACCGACACCGGTGCCGGGGGCGAAGATGAGGGGATCCCTATTCGTCCgg                                   >  1:863529/1‑151 (MQ=255)
                                                                                                                     aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCGGGGGCGGGCGCGAAGATGAGGGGATCCCTATTCGTCCggg                                  >  2:906296/1‑151 (MQ=255)
                                                                                                                          ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc                               <  1:1108849/149‑1 (MQ=255)
                                                                                                                          ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACa                              <  2:959370/150‑1 (MQ=255)
                                                                                                                          ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACa                              <  2:1331937/150‑1 (MQ=255)
                                                                                                                          ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAg                             <  1:351877/151‑1 (MQ=255)
                                                                                                                              aaaCGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGCACGCCAGAAATCTAATGCCCCCATGATCATCAACGATCCGCGCCATACCCACACCGGGGCCGGGGGCGAAGATGAGTGGGACCCTATTCGTCCgggga                                >  2:265358/1‑142 (MQ=255)
                                                                                                                                aCGCGAATGAGTGGCGGTGGGGGGACTTACTTTCTTGAACAAGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAACGACACCGGGGCCGGGCGCGAAGATGAGTGGGTTCCTATTCGTCCCGGGAAAGATGcc                        >  1:1256015/1‑150 (MQ=255)
                                                                                                                                  gcgAATGAGTGGCGGTGGGGTGGCTTTCTATCTTTAAAAAGCACGCCCGAAATCTAATGGCCCCATGATCATCATCGATCCCCGCTATACCGACACCGGTGGCGGGCGCGAAGATGAGTGGATCCCCATTCTTCCCGGGACAGATGcc                        >  2:2104026/1‑148 (MQ=255)
                                                                                                                                  gcgAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAAAGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGGTCCCTATTCGTCCGGGAACAGATGCCg                       >  2:1456855/1‑149 (MQ=255)
                                                                                                                                    gAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACt                   <  1:1042291/151‑1 (MQ=255)
                                                                                                                                    gAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACt                   <  1:2156890/151‑1 (MQ=255)
                                                                                                                                        gAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGt                <  2:797491/150‑1 (MQ=255)
                                                                                                                                            ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg           >  1:1315181/1‑151 (MQ=255)
                                                                                                                                                  ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg      <  1:1411734/150‑1 (MQ=255)
                                                                                                                                                    ggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCg   >  1:1389657/1‑151 (MQ=255)
                                                                                                                                                     gTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  <  1:812865/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941563‑941862

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: