Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,578,3010AC40.6% 177.5 / 16.9 96G142G (GGT→GGG) gntRd‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (49/8);  new base C (1/38);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 1.09e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGA  >  NC_000913/3578154‑3578445
                                                                                                                                                        |                                                                                                                                                
gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑Accc                                                                                                                                               <  1:814753/151‑1 (MQ=255)
     gaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCAc                                                                                                                                          >  1:1099778/1‑151 (MQ=255)
     gaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCAc                                                                                                                                          >  2:1355432/1‑151 (MQ=255)
        gTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAC‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCCCCACGGGAAT‑ACCCGCCACtt                                                                                                                                        <  1:6510/150‑1 (MQ=255)
        gTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTc                                                                                                                                       >  2:369735/1‑151 (MQ=255)
             ccGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAatc                                                                                                                                   >  2:230383/1‑150 (MQ=255)
             ccGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAatc                                                                                                                                   >  1:1263089/1‑150 (MQ=255)
             cagagaTAGGCAGTGTGGCGATCCCCGCGAGAAATAATGGCAGTGGTCATCTGGCGTGCTGCTTAAAAGTTATCAAAC‑CCGA‑CGGCAATATCAAGGCGTGGCGTCTTGCTGCCCAC‑CATTC‑CCCCCAGGGAAAT‑CCCCGCCAc                                                                                                                                          <  1:1376985/141‑1 (MQ=255)
                agaTAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCt                                                                                                                               >  1:750489/1‑151 (MQ=255)
                  aTAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCttt                                                                                                                             >  1:1027285/1‑151 (MQ=255)
                    aGGCAATGTGGCGATG‑CCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTa                                                                                                                            >  2:953413/1‑149 (MQ=255)
                     ggCAATGGGGCGATGCCCGCAAGAAAAAGGGGCGTTGTTCTTCGGGCGGGCTGCTTAAAATTTAAAAAAC‑CCGC‑GGGCAATATAAGGGCAGGGCGTTTTGCTGTCCAT‑CATTC‑CCCCCAGGGGAAT‑ACCCGCCACTTCAATCATCTTTaa                                                                                                                           <  1:198830/150‑1 (MQ=255)
                     gcaaagggggCGATGCCGGCGAAAAAAAAGGGCAGGGGTCTTCGGGCGTGCTGCTTAAAATTTAAAAAAC‑CCGA‑CGGCGATATCAGGGCGTGGCGATTTGCTCTCCAT‑CTGTC‑CCCCCGCGGAAAT‑CCCCGCCACTTCAATCATCTTTAAg                                                                                                                          <  1:262141/143‑1 (MQ=255)
                       cAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAgg                                                                                                                         >  1:494132/1‑150 (MQ=255)
                       cAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAgg                                                                                                                         >  1:386576/1‑150 (MQ=255)
                        aaTGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTg                                                                                                                       >  1:723912/1‑151 (MQ=255)
                         aTGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTgc                                                                                                                      >  2:847727/1‑151 (MQ=255)
                           ggggCGAGGCCCGCGAGAAATAAGGGCAGTGGTCATCTGGCGGGCTGCTTAAAATTAAAAAAAC‑CCGA‑CGGCGATATCAAGGCAGGCCGATTTGCTTTCCAT‑CGTTC‑CCCCCGGGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTgc                                                                                                                      <  1:1183892/147‑1 (MQ=255)
                             ggCGATCCCCGCGAGCAATAAGGGCAGTGGTCATCTGGCGTGCTCCTTAAAAGTTATAAAAA‑CCGC‑CGGCAATATCAAGGCATGGCGACTTGCTGCCCAT‑CAGTC‑CCACCACGGAAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                                                   <  2:1899062/150‑1 (MQ=255)
                                 aTGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgtg                                                                                                              >  2:115417/1‑151 (MQ=255)
                                 aTGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCAACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgtg                                                                                                              >  1:1446197/1‑151 (MQ=255)
                                  tccccGCGCGCAATAAGGGATGTGGTCATTGGGCGTAGTGCTTAAAATTTACCAAAC‑CCGG‑CGGGAATATCAACCCCTGCCTTTTTGCTGTCATT‑CCTTC‑CCCCCCCGGCATT‑CCCCGCCACTTCAATCATCTTTAAGGTGCgcgg                                                                                                                  <  1:422391/144‑1 (MQ=255)
                                       gcgAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGt                                                                                                         >  2:692479/1‑150 (MQ=255)
                                            cAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTc                                                                                                   >  1:1217217/1‑151 (MQ=255)
                                             aatagtGGCAGGGTTCATCGGGCGTGCTTCTAAAAAGTTAAAAAAC‑CCGC‑GGGCGATATCAAGGCAGGGCGTTTTGTTGTCCAT‑CATTC‑CCCCCAGGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTAc                                                                                                       <  1:771084/146‑1 (MQ=255)
                                             aataatGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCg                                                                                                  >  2:242272/1‑151 (MQ=255)
                                               aaagggCAGTGGTCATCGGGCTTGCTGCTAAAAATTAAAAAAAC‑CCGA‑CGGCAATAAAAAGGCAGGGCGATTTGCTGCCCAT‑CATTC‑CCCCCACGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGt                                                                                                <  2:820539/147‑1 (MQ=255)
                                                aatGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTg                                                                                               >  1:1184120/1‑151 (MQ=255)
                                                   ggCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAgg                                                                                            >  2:1833274/1‑151 (MQ=255)
                                                    gcgggggTCTTTTGGCGTCCTCCTAAAAAGTAAAAAAAC‑CCGA‑CGGCAATATCAAGGCAGGGCGATTTGCTTTCCAT‑CATTC‑CCACCACGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAg                                                                                             <  1:1325626/144‑1 (MQ=255)
                                                     cAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGAt                                                                                          >  1:570045/1‑151 (MQ=255)
                                                      aGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                         >  2:246134/1‑151 (MQ=255)
                                                      aGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                         >  2:723221/1‑151 (MQ=255)
                                                      aGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                         >  1:839723/1‑151 (MQ=255)
                                                        tGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                         >  1:606937/1‑149 (MQ=255)
                                                         gttCATCGGGCTTGCTGCTTCAAAGTAAAAAAAC‑CCGA‑CGGCAATATCAAGGGATGGCGCCTTGCTGCCCAT‑CATTC‑CCCCCAGGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCa                                                                                        <  1:1247987/147‑1 (MQ=255)
                                                           tcttCTGGCGTGCTGCTAAAAATTTACAAAAC‑CCGA‑CGGCAATACAAAGGCAGGGCGCTTTTCTGTCCAA‑CAGTC‑CCCCCAGGGGAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAg                                                                                       <  2:910208/145‑1 (MQ=255)
                                                           tCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGcc                                                                                    >  2:1490693/1‑151 (MQ=255)
                                                            cATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCg                                                                                   >  2:272407/1‑151 (MQ=255)
                                                              tCTGGCGGGCTGCTTCAAA‑TTACCAAAACCCGC‑CGGAGAAAGCAAGCCAGGGGGATTTGGTATCCAT‑CACTC‑CCACCACGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGt                                                                                                <  2:632807/136‑1 (MQ=255)
                                                                tGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCg                                                                                >  2:1481265/1‑150 (MQ=255)
                                                                 ggtgTGCTCCTAAAAATTTATCAAAC‑CCGC‑CGAAGATATCAAGCCATGGCTATTTGCTTCCCAT‑CAGTC‑CCCCCGCGGAATT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCg                                                                                <  1:1492706/146‑1 (MQ=255)
                                                                 ggCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGat                                                                              >  2:1595363/1‑151 (MQ=255)
                                                                 ggCGTGCTGCTAAAAATTTATAAAAC‑CCGA‑CGGCGATAAAAGGGCAGGGCGTTTTGCTCTCCAT‑CAGTC‑CCCCCAGGGGAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCg                                                                                <  1:260151/149‑1 (MQ=255)
                                                                 ggCGTGCTGCTAAAAATTAAAAAAAC‑CCGC‑GGGCGATAAAAAGGCAGGGCGCTTTGCTGCCCAT‑CAGTC‑CCCCCGGGGGAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCg                                                                                <  1:1242145/149‑1 (MQ=255)
                                                                        tgctTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTcc                                                                         >  2:125201/1‑149 (MQ=255)
                                                                         gctTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                                      >  1:692656/1‑151 (MQ=255)
                                                                         gctTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                                      >  1:1683858/1‑151 (MQ=255)
                                                                         gctTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                                      >  2:1948812/1‑151 (MQ=255)
                                                                            tCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                                      >  2:444966/1‑148 (MQ=255)
                                                                             caaaatttaaaaaa‑CCGCCCGGCAATATCAAGCAATGGCGAATTGC‑GTTCATACATTT‑CCACCAGGGAAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAg                                                                       <  2:1926946/137‑1 (MQ=255)
                                                                              aaaGTTAAAAAAC‑CCGG‑CGGGGATAACAAGGCGTGGCGTTTTGCTTCCCAT‑CATTC‑CCCCCGGGGGAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAg                                                                       <  1:470917/145‑1 (MQ=255)
                                                                                 gTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATg                                                              >  1:848042/1‑151 (MQ=255)
                                                                                   taccAAAA‑CCGA‑CGGTAATATAAGGGTAACGGGACTCGCTCCACAC‑AAGTC‑CCACCGCGGGCTT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATgg                                                             <  2:496712/147‑1 (MQ=255)
                                                                                   aataaaac‑CCGA‑GGGAGATACCAGGGCAGGGCGACTTCTTGTCCGT‑CTTTC‑CCCCCACGGGAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGa                                                            <  1:488980/147‑1 (MQ=255)
                                                                                    aacAAAC‑CCGG‑GGGCGATACCAGGGCTGGGCGTTTTGTTGCCCAT‑‑AGTTCCCCCCAGGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATg                                                              <  1:222653/146‑1 (MQ=255)
                                                                                    aTCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAt                                                           >  1:1591550/1‑151 (MQ=255)
                                                                                       aaac‑CCGC‑CGGCGATATAAGGGCAGGGCGCTTTGCTGCCAAT‑CATTC‑CCCCCACGGGAAC‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAtt                                                          <  1:950681/149‑1 (MQ=255)
                                                                                       aaac‑CCGA‑CGGCAATATCAGGGCGGGGCGCCTTGCTCCCCAT‑CATTC‑CCCCCAGGGGAAA‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAtt                                                          <  2:398395/149‑1 (MQ=255)
                                                                                       aaac‑CCGA‑CGGCAATATAAAGGCAGGGCGTTTTGCTGTCCAC‑CAGTC‑CCCCCACGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGa                                                            <  2:1179482/147‑1 (MQ=255)
                                                                                       aaaa‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCg                                                        >  1:905852/1‑151 (MQ=255)
                                                                                          aacccgacgtaaatatCAAGGCTTGGCGACTTGCTTCCCAT‑CATTC‑CCCCCAGGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCg                                                        <  1:192183/138‑1 (MQ=255)
                                                                                             ccgccggaaatacCAAGGCGGGGCGTCTTGTTCCCAAT‑CAGTC‑CCCCCGCGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAgg                                                     <  1:530058/141‑1 (MQ=255)
                                                                                              gA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGt                                                  >  1:1961261/1‑151 (MQ=255)
                                                                                                  ggCGATATCAAGGCATGGCGATTTGCTGCCCAT‑CAGTC‑CCCCCACGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTc                                                <  1:925176/150‑1 (MQ=255)
                                                                                                   gCGATAGCAAGGCATGGCGACTTGCGGCCCAT‑CAGTC‑CCCCCGGGGGAAG‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCtt                                              <  1:1493336/151‑1 (MQ=255)
                                                                                                       tatCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTc                                          >  1:248035/1‑151 (MQ=255)
                                                                                                          cAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAt                                       >  2:1775601/1‑151 (MQ=255)
                                                                                                             ggcggggCGCCTTGCTGCCCAC‑CTGTC‑CCCCCGGGGGAAC‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAt                                       <  1:1909419/143‑1 (MQ=255)
                                                                                                             ggCAGGGCGATTTGCTTTCCAT‑CAGTC‑CCCCCAGGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTc                                    <  2:1369377/151‑1 (MQ=255)
                                                                                                               cATGGCGTCTTGCTGCCCAT‑CAGTC‑CCCCCAGGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCg                                  <  2:52397/151‑1 (MQ=255)
                                                                                                                 tGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGt                                >  2:1852470/1‑151 (MQ=255)
                                                                                                                  ggCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGtt                               >  1:301271/1‑151 (MQ=255)
                                                                                                                     gACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtat                            >  1:1759763/1‑151 (MQ=255)
                                                                                                                      aCTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtat                            >  1:973604/1‑150 (MQ=255)
                                                                                                                       cTTGCTGTCCAT‑CAACT‑CCACCACCGGGAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGTGGCGTGTGAGTACGTTCGGTGAGGATCAGTCCATCGATATTCCATGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtata                           >  1:1939524/1‑150 (MQ=255)
                                                                                                                        ttGCTGTCCAC‑CATTC‑CCCCCACGGGAAC‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTta                             <  2:272376/147‑1 (MQ=255)
                                                                                                                         tgctgcCCAT‑CATTT‑CCCCCGGGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATaa                          <  2:686648/149‑1 (MQ=255)
                                                                                                                          ttctttcccaccgTT‑CCCCCCCGGGAATCCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtat                            <  2:59519/136‑1 (MQ=255)
                                                                                                                          gCTGCCCAT‑CAGTC‑CCCCCAGGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATaa                          <  2:905868/148‑1 (MQ=255)
                                                                                                                           cTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGt                      >  2:773154/1‑151 (MQ=255)
                                                                                                                           cTGCCCAT‑CTGTC‑CCACCAGGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAcc                        <  2:1102626/149‑1 (MQ=255)
                                                                                                                            tccccAT‑CATTT‑CCCCCAGGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAc                         <  2:494141/144‑1 (MQ=255)
                                                                                                                                cAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGACGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGCTTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTgc                 >  1:986168/1‑151 (MQ=255)
                                                                                                                                    cAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGcc              >  1:1311438/1‑151 (MQ=255)
                                                                                                                                     cgttcccaccgCGGGAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTgcgc               <  2:255938/145‑1 (MQ=255)
                                                                                                                                         ccccccagGGAAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGt                      <  1:1694509/135‑1 (MQ=255)
                                                                                                                                          ccaccaCGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCAt         >  2:1186137/1‑151 (MQ=255)
                                                                                                                                           ccccagGGGATT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATg        <  2:570054/149‑1 (MQ=255)
                                                                                                                                           caccacGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATg        >  1:1647603/1‑151 (MQ=255)
                                                                                                                                             ccacGAGAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCa             <  1:1067583/144‑1 (MQ=255)
                                                                                                                                              cacGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTc     >  1:1695380/1‑151 (MQ=255)
                                                                                                                                               agggaaaT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTc     <  1:1313626/145‑1 (MQ=255)
                                                                                                                                                 gggAAT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCt    <  1:869647/149‑1 (MQ=255)
                                                                                                                                                 ggaaaT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATg        <  2:1194368/142‑1 (MQ=255)
                                                                                                                                                 ggaaaT‑CCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGa  <  2:1728313/148‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                
GTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGA‑CGGCGATATCAAGGCATGGCGACTTGCTGTCCAT‑CAGTT‑CCACCACGGGAAT‑ACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGA  >  NC_000913/3578154‑3578445

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