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

GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941563‑941862
                                                                                                                                                      |                                                                                                                                                     
gCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                                                                                       <  2:637165/151‑1 (MQ=255)
gCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                                                                                       <  1:351196/151‑1 (MQ=255)
 cGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTg                                                                                                                                                      <  2:1820677/151‑1 (MQ=255)
         cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTAc                                                                                                                                                <  2:1609887/149‑1 (MQ=255)
            ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCt                                                                                                                                           <  1:1969952/151‑1 (MQ=255)
             cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTTCTATCtt                                                                                                                                          >  2:674501/1‑151 (MQ=255)
               cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTAtttt                                                                                                                                          >  2:1552920/1‑149 (MQ=255)
                gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGAAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGTGGCctactaaat                                                                                                                                           >  1:676575/1‑144 (MQ=255)
                  aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTCTATCTTTaaa                                                                                                                                      >  1:982843/1‑149 (MQ=255)
                    aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAgg                                                                                                                                   <  2:2031889/151‑1 (MQ=255)
                    aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAgg                                                                                                                                   <  2:447046/151‑1 (MQ=255)
                      tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGc                                                                                                                                  <  2:939389/150‑1 (MQ=255)
                        cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTACTATTTTTaaaaaggac                                                                                                                                >  1:1434939/1‑143 (MQ=255)
                         ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGc                                                                                                                              <  1:1429747/151‑1 (MQ=255)
                         ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTTCTATTTTTAAAAGGCa                                                                                                                                 >  2:1327420/1‑148 (MQ=255)
                         ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGAACTACAATTTTTAACAAGCAc                                                                                                                                >  2:1476957/1‑149 (MQ=255)
                          gAAGGTTTGAACTATACCTACGGCGGCTGGACAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGcc                                                                                                                             <  2:1149604/151‑1 (MQ=255)
                           aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATTTTTAAAAAGCCCGcc                                                                                                                             >  2:2016433/1‑150 (MQ=255)
                           aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGAATTACTATCttttaaaaggaccccc                                                                                                                            >  2:493644/1‑137 (MQ=255)
                            aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAACAGGCCCGcc                                                                                                                             >  2:210815/1‑149 (MQ=255)
                             gttttGAACTGTATCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCTAAAACAGTAAGCTGGTAGTGCTGTTTGGTGATAACCCTGGCGAAACGCGAATGAGTGTCCGTGTGGTGACTTACTATCTTGAACAGGCACTCCat                                                                                                                           <  1:2033417/148‑2 (MQ=255)
                                ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaaa                                                                                                                        <  1:2050282/150‑1 (MQ=255)
                                ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaaa                                                                                                                        <  2:243345/150‑1 (MQ=255)
                                 tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGCTTACTATCTTTAACAAGCACGCCAGaaa                                                                                                                        >  1:79636/1‑149 (MQ=255)
                                    aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTaa                                                                                                                   <  1:881029/151‑1 (MQ=255)
                                    aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTaa                                                                                                                   <  1:1260113/151‑1 (MQ=255)
                                    aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGTGGATTACTATCTTGAAAAAGAACGCCCGAAAACTa                                                                                                                    >  2:1379341/1‑150 (MQ=255)
                                          ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCCTACTATCTTTAACAAGGACGCCAGAAATCTAATGccc                                                                                                              >  2:1023223/1‑150 (MQ=255)
                                          ccTACGGCGGCTGGGCAGAAGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCAGGCGAAACGCGAATGAGGGGCGGAGGGGGGGCGTATCAAATTGAACAGGGACGCCCGaaaggtgcatcac                                                                                                              >  2:22770/1‑140 (MQ=255)
                                            tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTTTTATCTTTAACAAGGACGCCCGAAAAATAATAccccc                                                                                                            >  2:945467/1‑148 (MQ=255)
                                                 cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAAGCACGCCAGAAATCTAAAGGCCGccaga                                                                                                        >  2:1427071/1‑145 (MQ=255)
                                                 cggcTGGGCAGATGGCAACAGCCCGTCGGAGATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatc                                                                                                      <  2:1153040/151‑1 (MQ=255)
                                                  ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatca                                                                                                     <  2:1995311/151‑1 (MQ=255)
                                                  ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatca                                                                                                     <  2:324943/151‑1 (MQ=255)
                                                  ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatca                                                                                                     <  1:378163/151‑1 (MQ=255)
                                                   gcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcat                                                                                                    <  2:358265/151‑1 (MQ=255)
                                                    cTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAGGCACGCCCAAAAAATAATGCCCGc                                                                                                            >  1:1899869/1‑142 (MQ=255)
                                                         cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAAGCACGCCAGAAAACTAATGCCCCCa                                                                                                           >  2:1364294/1‑138 (MQ=255)
                                                          aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGGTTACTATCTTTAAAAGGGACCCCCGAAAAATAATTACCGCAAGattattat                                                                                                 >  2:730234/1‑147 (MQ=255)
                                                          aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatc                                                                                                   >  2:366402/1‑145 (MQ=255)
                                                              ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgc                                                                                         <  2:1217653/151‑1 (MQ=255)
                                                               gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTGAAAAGGCACGCCCGAAATCTAATGCCCGCCTGATCATCATCGCACcgc                                                                                         >  1:805653/1‑150 (MQ=255)
                                                                  aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCCTTCTATTTTGAACAGGCCCCCCCCACATTTAAAGCCCCCACGatccac                                                                                                   >  2:638268/1‑134 (MQ=255)
                                                                     gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAc                                                                                  <  1:1542167/151‑1 (MQ=255)
                                                                       ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGGACGCCAGAAATCTAAAGCCCGCCTGATTATCATCGATCCCCGCTATAcc                                                                                 >  1:393387/1‑150 (MQ=255)
                                                                        cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAcc                                                                                 <  2:1753034/149‑1 (MQ=255)
                                                                         gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGac                                                                              <  1:1889839/151‑1 (MQ=255)
                                                                         gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTCTATCTTTAACAGGGACGCCAGAAATCTAATTGCCCCATTATCATCATCGTTCCGCCCTATAcc                                                                                 >  1:1459618/1‑148 (MQ=255)
                                                                           cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGaca                                                                             <  1:443971/150‑1 (MQ=255)
                                                                           cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTATTATCTTTAAAAGGGACGCCAGAAAACTAATGGCCGCCTGATCATCATCCATCCGCGCtataac                                                                                 >  1:2060335/1‑144 (MQ=255)
                                                                                atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTg                                                                       <  2:408456/151‑1 (MQ=255)
                                                                                   gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTGAACAGGCACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCTATACCCACACCGGGGcc                                                                     >  2:844394/1‑150 (MQ=255)
                                                                                    aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCgg                                                                   <  1:1646328/151‑1 (MQ=255)
                                                                                      aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTTAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCggg                                                                  >  2:821036/1‑150 (MQ=255)
                                                                                       aCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcg                                                                <  1:1014765/151‑1 (MQ=255)
                                                                                         aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGcgcg                                                              <  2:1367897/151‑1 (MQ=255)
                                                                                         aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAAACTAATGCCCGCATGATTATCATCGATCCGCGCTATACCGACACCCGGGCCGGGggcg                                                              >  2:785804/1‑151 (MQ=255)
                                                                                            aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATTTTTAACAGGCACGCCACAAATCTAATTGCCGCATGATCATCCTCGGTCCGCGCTTTACCCACACCGGTGCCGGGggcg                                                              >  1:1321581/1‑148 (MQ=255)
                                                                                            aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGCTTACCATCTTTAAAAGGCACGCCCGAAATCTAATGCCCCCATGATCATCAACCAACCGCGCCAAACCGACACccgggccggg                                                                  >  2:1040168/1‑144 (MQ=255)
                                                                                             aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGATGGGGGAGATATTAAATTTAAAAAGCACGCCCGAAATCTAAAAGACTCATGATTATCAATGAGCCGCGCTATACAGACACAGGGGGCGGGcgcg                                                              >  1:1576189/1‑147 (MQ=255)
                                                                                              gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                         <  1:1368655/151‑1 (MQ=255)
                                                                                              gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                         <  2:571752/151‑1 (MQ=255)
                                                                                              gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                         <  1:950394/151‑1 (MQ=255)
                                                                                               cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCCCGCCCGAAAAATAAAGCCCGCATGATCATCATCGATACGCGCCATAACGACACCGGGGGCGGGGGCGgaga                                                          >  1:975503/1‑149 (MQ=255)
                                                                                               cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAAACTAATTCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGGCGGGCGCGGAGAt                                                         >  2:1449907/1‑150 (MQ=255)
                                                                                                       gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATc                                                <  1:1021871/151‑1 (MQ=255)
                                                                                                       gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATc                                                <  1:355482/151‑1 (MQ=255)
                                                                                                            ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTa                                            <  2:1442996/150‑1 (MQ=255)
                                                                                                              ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACCTACTATCCTTAAAAAGCACCCCCGAAAACTAATGCCCCCATGAATATCCTCCATCCGCGCCATACCCACACCGGGGCCCGGGGCGAAGATGaaggggtcc                                               >  1:1404712/1‑137 (MQ=255)
                                                                                                               ggaataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGGGGagcaagaggag                                        <  2:829914/149‑11 (MQ=255)
                                                                                                                taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCCTACTATTTTTAACAGGCACGCCAGAAATCTAATGCCCCCATGATTATCATCGCTCCGCGCTATACCGACACCCGGGCCGGGCGCGAAGATGAGTGGATCCCCAtttct                                       >  2:1231864/1‑148 (MQ=255)
                                                                                                                 aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTATCTTTAACAGGCACGCCCGAAAACCAATGGCCGCATGGGCAGCATCGCTCCGCGCTATACCCACACCGGTGCCGGGCGCGCAGATGAGGGGGGCTCTATTCGt                                       >  2:1358106/1‑150 (MQ=255)
                                                                                                                   taaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTAATTTTAAAAAGCACGCCAGAAAACTAAATCCCGCATGATCATCCTCGGTCCCCGCCATTCCGACACCCGGGGCGGGGGCGGAGAAGAGTGGAGCCCCTTTCg                                        >  1:1179366/1‑147 (MQ=255)
                                                                                                                    aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTTTATCTTTAACAGGCACGCCCGAAAAATAATGCCCGCCCGATCATCATCGATCCGCGCTATACCCCCACCGGTGCCGGGGGCGAAGATGAGGGGGTCCCCATTCGTCCgg                                   >  1:1680648/1‑151 (MQ=255)
                                                                                                                    aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTACTATCCTGAAAAAGCACGCCAGAAAAGTAATGCCCGCCTGAGCATCATCCATCCGCGCTATACCCACACCCGTGGCGGGGGCGAAGATGAGTGGATCTCCATTCGTc                                      >  1:1064379/1‑148 (MQ=255)
                                                                                                                     aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                                  <  1:1250853/151‑1 (MQ=255)
                                                                                                                      cccTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTTAACAGGGACGCCAGAAAACTAAAAACCCCCTGATCATCATCGGTCCGCGCTATACCGCCACCCGGGGCCGGCGGGGAGAAGAGGGGGTCCCCATTCGTCCcggg                                 >  1:139808/1‑150 (MQ=255)
                                                                                                                       ccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTGAACAGGCACGCCAGAAAAACAATGCCCCCATGATCATCATCGCTCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGGGGATCCCTCTTCGTCCggg                                  >  1:388377/1‑149 (MQ=255)
                                                                                                                         tGGCGAAACGCGAATGAGTGGCGGAGGGGTGGGTTTCTATCTTTAAAAGGGCAGCCAGAAAACCAAAGCCCGCCCGTAGATCATATTTCCGCGCCTTAACAACACCGGGGCCGGGCGCCAAAAAGTGTGGGCTCCTATTTGCCCggg                                  >  1:1589349/1‑147 (MQ=255)
                                                                                                                          ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc                               <  2:438345/149‑1 (MQ=255)
                                                                                                                             gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATg                          <  2:1434987/151‑1 (MQ=255)
                                                                                                                             gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATg                          <  2:1900371/151‑1 (MQ=255)
                                                                                                                             gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATg                          <  1:504702/151‑1 (MQ=255)
                                                                                                                             gAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGACAGa                            >  2:286011/1‑149 (MQ=255)
                                                                                                                                 cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCg                       <  2:1231199/150‑1 (MQ=255)
                                                                                                                                  gcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCa                     <  2:1741423/151‑1 (MQ=255)
                                                                                                                                     aaTGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACt                   <  1:546987/150‑1 (MQ=255)
                                                                                                                                     aaTGAGTGGCGGTGGGGGGACTTACTATTTTGAAAAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGGAGATGGGGGGGTTCCTTCTCTTCCCGGGACAGaaggccccc                    >  1:1917072/1‑141 (MQ=255)
                                                                                                                                      aTGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGTATCCCTATTCGTCCGGGAACAGATGCCGCACTgg                 <  1:1693294/151‑1 (MQ=255)
                                                                                                                                            ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc            <  1:1907641/150‑1 (MQ=255)
                                                                                                                                            ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg           >  1:1200479/1‑151 (MQ=255)
                                                                                                                                             gCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACgg          >  2:1747464/1‑151 (MQ=255)
                                                                                                                                             gCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACgg          <  1:998222/151‑1 (MQ=255)
                                                                                                                                             gCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACgg          >  2:111844/1‑151 (MQ=255)
                                                                                                                                              cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt         <  1:1093644/151‑1 (MQ=255)
                                                                                                                                               ggTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTc        >  2:1884123/1‑151 (MQ=255)
                                                                                                                                                  ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg      <  2:730657/150‑1 (MQ=255)
                                                                                                                                                   gggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc    <  2:1449341/151‑1 (MQ=255)
                                                                                                                                                      tGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  2:998962/1‑150 (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: