Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7130TG38.2% 120.8 / 20.5 77V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (10/37);  new base G (29/0);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 8.18e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941566‑941862
                                                                                                                                                   |                                                                                                                                                     
actactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGAc                                                                                                                                                    <  1:1047621/151‑1 (MQ=255)
  tactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGGCtt                                                                                                                                                  >  2:295091/1‑151 (MQ=255)
   actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt                                                                                                                                                   <  1:1047411/149‑1 (MQ=255)
    ctcctcCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt                                                                                                                                                   <  1:668467/148‑1 (MQ=255)
         ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAt                                                                                                                                             <  2:990940/149‑1 (MQ=255)
          cgcgCAAATTGTGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCtt                                                                                                                                          <  2:498839/151‑1 (MQ=255)
          cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCtt                                                                                                                                          >  2:887024/1‑151 (MQ=255)
          cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCtt                                                                                                                                          >  2:1219151/1‑151 (MQ=255)
             gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTCCTAttttta                                                                                                                                        >  2:1176633/1‑148 (MQ=255)
              cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTGaa                                                                                                                                       >  2:1239145/1‑150 (MQ=255)
              cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTaa                                                                                                                                       >  2:415969/1‑150 (MQ=255)
                aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAg                                                                                                                                    <  1:1036818/151‑1 (MQ=255)
                     cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACg                                                                                                                               <  2:35232/151‑1 (MQ=255)
                        aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCa                                                                                                                            <  2:586139/151‑1 (MQ=255)
                             ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaaa                                                                                                                        <  1:1373650/150‑1 (MQ=255)
                             ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGCGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGCTTAATATCTTGCACAAGCagcccccaa                                                                                                                         >  2:990341/1‑141 (MQ=255)
                              tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAAt                                                                                                                       <  1:621270/150‑1 (MQ=255)
                              tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTATTATTTTTAACAGGCCAGCCAGcaa                                                                                                                        >  1:804390/1‑149 (MQ=255)
                               gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCt                                                                                                                     <  2:1325612/151‑1 (MQ=255)
                               gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCt                                                                                                                     <  2:1004704/151‑1 (MQ=255)
                                 aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTaa                                                                                                                   <  1:344614/151‑1 (MQ=255)
                                    ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGc                                                                                                                <  2:833187/151‑1 (MQ=255)
                                        cTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCTAGAAATCTAATGccc                                                                                                              <  1:1153228/149‑1 (MQ=255)
                                        cTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATATCTTGAACAGGCACGCCAGAAAACTAATGccccc                                                                                                            >  2:203061/1‑149 (MQ=255)
                                            ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGGACGCCAGAAAACTAATGgcc                                                                                                              >  2:20909/1‑145 (MQ=255)
                                            ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAAGCACGCCAGGAATCTAATGCCCGc                                                                                                            >  1:1355108/1‑147 (MQ=255)
                                               ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatca                                                                                                     <  1:1043050/151‑1 (MQ=255)
                                               ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGccccc                                                                                                            >  2:127482/1‑142 (MQ=255)
                                                 cTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTTAACAAGCACGCCCGAAATCTAATGCCCGCACGatc                                                                                                      >  1:1169212/1‑148 (MQ=255)
                                                     gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCCGAAATCTAATGCCCGCATGatcatcatc                                                                                                >  2:1315353/1‑150 (MQ=255)
                                                       aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCCTTATATCTTGAACAGGCACGCCCGAAATCTAATGGCCGCAt                                                                                                          >  1:827438/1‑138 (MQ=255)
                                                         aTGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATcc                                                                                           <  1:904152/151‑1 (MQ=255)
                                                         aTGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATcc                                                                                           <  1:970428/151‑1 (MQ=255)
                                                            gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgcg                                                                                        <  1:273387/151‑1 (MQ=255)
                                                            gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTCCCATCTTTAACAGGCACCCCCGCAATCTAATGCCCCCCTGATCCTCATCGCTcccc                                                                                         >  1:145227/1‑148 (MQ=255)
                                                                   cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAcc                                                                                 <  2:20655/151‑1 (MQ=255)
                                                                      gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGACTTAATATCTTTAACAAGCACGCCCGAAATCTAATTCCCCCATGATCATCATCCATCCGCGCTATAcc                                                                                 >  1:506243/1‑148 (MQ=255)
                                                                      gTCAGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAGTGAGTGGCGGTGGGGTGACATACTAGCCTGACCAGGCACGCCAGAAATCTGATACCCGCATGATCATCAACGATCCGCGTTATAcca                                                                                <  1:964128/149‑2 (MQ=255)
                                                                        cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTTAACAGGCACGCCCGAAAACTAATGCCCGCATGATCATCATCGATCCCCGCTATACCcaca                                                                             >  1:1092480/1‑150 (MQ=255)
                                                                           atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTTAACAAGCACGCCCGAAAAATAATGGCCGCATGATAATCTTTCATCCGCGCGATACCGACAccc                                                                          >  1:1270504/1‑149 (MQ=255)
                                                                            tatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCCATAACCACAcc                                                                           >  2:731429/1‑148 (MQ=255)
                                                                             atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGt                                                                        <  1:201369/150‑1 (MQ=255)
                                                                             atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTg                                                                       <  2:1073003/151‑1 (MQ=255)
                                                                               cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGcc                                                                     <  2:819515/151‑1 (MQ=255)
                                                                                gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCg                                                                    <  1:120184/151‑1 (MQ=255)
                                                                                gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGcc                                                                     >  2:1342356/1‑150 (MQ=255)
                                                                                   aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCCATAACCACACCGGGGCCggg                                                                  >  2:1122774/1‑150 (MQ=255)
                                                                                   aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCGGGGCCggg                                                                  >  2:873876/1‑150 (MQ=255)
                                                                                               gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGt                                                     <  1:1368334/151‑1 (MQ=255)
                                                                                                 aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTgg                                                   >  2:163088/1‑151 (MQ=255)
                                                                                                    gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAGGCACGCCAGAAAACTAATGCCCGCCTGATCATCATCGATCCCCGCTATACCCACACCGGGGCCGGGCGCGGAGATGAGTGGagc                                                >  2:119509/1‑149 (MQ=255)
                                                                                                      tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGAAGATGAGTGGATccc                                              >  2:645996/1‑151 (MQ=255)
                                                                                                        tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTAATATATTTGAAAGGCACGCCAGAAAACTAAAGCCCGCCTGATCCTCATCTATCCGCGCTATACCGACACCCGGGCGGGGGGCGAAGATGAGTGGGTCCCt                                             >  1:85871/1‑150 (MQ=255)
                                                                                                          tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt                                           <  1:713365/150‑1 (MQ=255)
                                                                                                            gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCg                                        <  1:844145/151‑1 (MQ=255)
                                                                                                              aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                                      <  2:655747/151‑1 (MQ=255)
                                                                                                              aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGGACGCCAGAAACCTAATGGCCGCATGATCATCATCGATCCGCGCTATACCGACAACGGTGCCCGGGGCGAAGATGAGTGGGTCCCCATTCGTc                                      >  1:1240192/1‑151 (MQ=255)
                                                                                                              aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCCGGGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                          >  2:868975/1‑147 (MQ=255)
                                                                                                              aataaCCCTGGCGAAACGCGAATGAGTGGCGGGGGGGCGACCTAATCTCTTTAACACGCGCGCCCGAAATCTTATGGCCGCATGATCATGATCGAACCGCGCTCTACCGACCCCCGTGCCGGGCGCGCAGAGGAGAGGGTCCCTATTCGTc                                      >  1:318773/1‑151 (MQ=255)
                                                                                                               ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTcc                                     <  1:254163/151‑1 (MQ=255)
                                                                                                                  aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCCTCATCGATCCGCGCTATAACCACACCGGGGCCGGGCGCGAAGATGAGTGGGTCCCTATTCGTCCggg                                  >  2:340186/1‑151 (MQ=255)
                                                                                                                   cccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGGCCGAAAACTAATGGCCGCATGATCATTATCGATCCGCGCTATACCGACACCGGGGGCGGGCGCGGAGATGAGGGGATCCCTATTCGTCCggg                                  >  1:918845/1‑150 (MQ=255)
                                                                                                                     cTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCgggga                                >  2:182603/1‑148 (MQ=255)
                                                                                                                      tGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTTCTATCTTCGACAGGCACGCCAGAAAAACAATGACCGGCTGGTAATCCACGATACGCGCCCAATCGCCACCGCTGCCGGGGGGGGAGAAGTGGGGACCCCCGTTCCGCCGGGAAc                               >  1:628405/1‑150 (MQ=255)
                                                                                                                       ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACa                              <  1:1019309/150‑1 (MQ=255)
                                                                                                                       ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAg                             <  2:1355988/151‑1 (MQ=255)
                                                                                                                         cGAAACGCGAATGAGTGGCGGTGGGGTGACTTCCTCTCTTTCACACGCACCCCAGAAATCTTATGCCCGCCTGCTCATCCTCGATCCTCGCCTTACCCCCCCCGCTGCCCGGGGCGAAGATGTGGGGGTCCCTCTTCGTCCGGGGAAAGa                            >  1:748163/1‑150 (MQ=255)
                                                                                                                         cGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACCCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCCATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGACa                              >  2:29556/1‑148 (MQ=255)
                                                                                                                           aaaCGCGAATGAGTGGCGGTGGGGTTACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGTTCATCATCGATCCGCGCTATACCAACTCCGGTGCCGCGCGCGAAGATGAGTGGATCCCTATTCATCCGTGTACAGAta                          <  1:1295728/150‑2 (MQ=255)
                                                                                                                            aaCGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGACAGATGcc                        >  2:1065906/1‑151 (MQ=255)
                                                                                                                                cgAATGAGTGGCGGGGGGGTGACTTTCTATCTTGAACCGGCACGCCAGAAATCTAATGCCCGCATGATCATAATCGACCCGCGCCATAACGACCACGGTGCCCGGCGCGAAGATGAGTGGAACCCTCTTCGTCCGGGAACAGATGTCGCAc                    >  1:622179/1‑151 (MQ=255)
                                                                                                                                cgAATGAGTGGCGGGGGGGTGACTTACTATCTTTAACAGGCACGCCCGAAATCTAATGCCCCCCTGATCACCATCGATCCGCGCCATACCCACACCGGTGCCGGGCGCGAAGAAGAGGGGGTCCCTATTCGGCCGGGAACAGATGCCGc                      >  1:657273/1‑149 (MQ=255)
                                                                                                                                  aaTGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTg                  <  2:174593/151‑1 (MQ=255)
                                                                                                                                          gCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACgg          <  2:1376711/151‑1 (MQ=255)
                                                                                                                                              tGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg      >  2:900187/1‑151 (MQ=255)
                                                                                                                                               ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg      <  2:433255/150‑1 (MQ=255)
                                                                                                                                               ggggTGACTTACTATCTTGAACAGGCACGCCAGAAAGCTAATGCCCGCATGATCAGCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg      <  1:338452/150‑1 (MQ=255)
                                                                                                                                                   tGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  1:847483/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941566‑941862

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: