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→GGGgntRd‑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

GTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGA  >  NC_000913/3578154‑3578445
                                                                                                                                                        |                                                                                                                                                
gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATAccc                                                                                                                                               <  1:814753/151‑1 (MQ=255)
     gaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCAc                                                                                                                                          >  1:1099778/1‑151 (MQ=255)
     gaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCAc                                                                                                                                          >  2:1355432/1‑151 (MQ=255)
        gTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAACCCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCCCCACGGGAATACCCGCCACtt                                                                                                                                        <  1:6510/150‑1 (MQ=255)
        gTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTc                                                                                                                                       >  2:369735/1‑151 (MQ=255)
             ccGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAatc                                                                                                                                   >  2:230383/1‑150 (MQ=255)
             ccGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAatc                                                                                                                                   >  1:1263089/1‑150 (MQ=255)
             cagagaTAGGCAGTGTGGCGATCCCCGCGAGAAATAATGGCAGTGGTCATCTGGCGTGCTGCTTAAAAGTTATCAAACCCGACGGCAATATCAAGGCGTGGCGTCTTGCTGCCCACCATTCCCCCCAGGGAAATCCCCGCCAc                                                                                                                                          <  1:1376985/141‑1 (MQ=255)
                agaTAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCt                                                                                                                               >  1:750489/1‑151 (MQ=255)
                  aTAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCttt                                                                                                                             >  1:1027285/1‑151 (MQ=255)
                    aGGCAATGTGGCGATGCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTa                                                                                                                            >  2:953413/1‑149 (MQ=255)
                     ggCAATGGGGCGATGCCCGCAAGAAAAAGGGGCGTTGTTCTTCGGGCGGGCTGCTTAAAATTTAAAAAACCCGCGGGCAATATAAGGGCAGGGCGTTTTGCTGTCCATCATTCCCCCCAGGGGAATACCCGCCACTTCAATCATCTTTaa                                                                                                                           <  1:198830/150‑1 (MQ=255)
                     gcaaagggggCGATGCCGGCGAAAAAAAAGGGCAGGGGTCTTCGGGCGTGCTGCTTAAAATTTAAAAAACCCGACGGCGATATCAGGGCGTGGCGATTTGCTCTCCATCTGTCCCCCCGCGGAAATCCCCGCCACTTCAATCATCTTTAAg                                                                                                                          <  1:262141/143‑1 (MQ=255)
                       cAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg                                                                                                                         >  1:494132/1‑150 (MQ=255)
                       cAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg                                                                                                                         >  1:386576/1‑150 (MQ=255)
                        aaTGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTg                                                                                                                       >  1:723912/1‑151 (MQ=255)
                         aTGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTgc                                                                                                                      >  2:847727/1‑151 (MQ=255)
                           ggggCGAGGCCCGCGAGAAATAAGGGCAGTGGTCATCTGGCGGGCTGCTTAAAATTAAAAAAACCCGACGGCGATATCAAGGCAGGCCGATTTGCTTTCCATCGTTCCCCCCGGGGGAATACCCGCCACTTCAATCATCTTTAAGGTgc                                                                                                                      <  1:1183892/147‑1 (MQ=255)
                             ggCGATCCCCGCGAGCAATAAGGGCAGTGGTCATCTGGCGTGCTCCTTAAAAGTTATAAAAACCGCCGGCAATATCAAGGCATGGCGACTTGCTGCCCATCAGTCCCACCACGGAAATACCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                                                   <  2:1899062/150‑1 (MQ=255)
                                 aTGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgtg                                                                                                              >  2:115417/1‑151 (MQ=255)
                                 aTGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCAACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgtg                                                                                                              >  1:1446197/1‑151 (MQ=255)
                                  tccccGCGCGCAATAAGGGATGTGGTCATTGGGCGTAGTGCTTAAAATTTACCAAACCCGGCGGGAATATCAACCCCTGCCTTTTTGCTGTCATTCCTTCCCCCCCCGGCATTCCCCGCCACTTCAATCATCTTTAAGGTGCgcgg                                                                                                                  <  1:422391/144‑1 (MQ=255)
                                       gcgAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGt                                                                                                         >  2:692479/1‑150 (MQ=255)
                                            cAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTc                                                                                                   >  1:1217217/1‑151 (MQ=255)
                                             aatagtGGCAGGGTTCATCGGGCGTGCTTCTAAAAAGTTAAAAAACCCGCGGGCGATATCAAGGCAGGGCGTTTTGTTGTCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTAc                                                                                                       <  1:771084/146‑1 (MQ=255)
                                             aataatGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCg                                                                                                  >  2:242272/1‑151 (MQ=255)
                                               aaagggCAGTGGTCATCGGGCTTGCTGCTAAAAATTAAAAAAACCCGACGGCAATAAAAAGGCAGGGCGATTTGCTGCCCATCATTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGt                                                                                                <  2:820539/147‑1 (MQ=255)
                                                aatGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTg                                                                                               >  1:1184120/1‑151 (MQ=255)
                                                   ggCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAgg                                                                                            >  2:1833274/1‑151 (MQ=255)
                                                    gcgggggTCTTTTGGCGTCCTCCTAAAAAGTAAAAAAACCCGACGGCAATATCAAGGCAGGGCGATTTGCTTTCCATCATTCCCACCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAg                                                                                             <  1:1325626/144‑1 (MQ=255)
                                                     cAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGAt                                                                                          >  1:570045/1‑151 (MQ=255)
                                                      aGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                         >  2:246134/1‑151 (MQ=255)
                                                      aGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                         >  2:723221/1‑151 (MQ=255)
                                                      aGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                         >  1:839723/1‑151 (MQ=255)
                                                        tGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                         >  1:606937/1‑149 (MQ=255)
                                                         gttCATCGGGCTTGCTGCTTCAAAGTAAAAAAACCCGACGGCAATATCAAGGGATGGCGCCTTGCTGCCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCa                                                                                        <  1:1247987/147‑1 (MQ=255)
                                                           tcttCTGGCGTGCTGCTAAAAATTTACAAAACCCGACGGCAATACAAAGGCAGGGCGCTTTTCTGTCCAACAGTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAg                                                                                       <  2:910208/145‑1 (MQ=255)
                                                           tCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGcc                                                                                    >  2:1490693/1‑151 (MQ=255)
                                                            cATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCg                                                                                   >  2:272407/1‑151 (MQ=255)
                                                              tCTGGCGGGCTGCTTCAAATTACCAAAACCCGCCGGAGAAAGCAAGCCAGGGGGATTTGGTATCCATCACTCCCACCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGt                                                                                                <  2:632807/136‑1 (MQ=255)
                                                                tGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCg                                                                                >  2:1481265/1‑150 (MQ=255)
                                                                 ggtgTGCTCCTAAAAATTTATCAAACCCGCCGAAGATATCAAGCCATGGCTATTTGCTTCCCATCAGTCCCCCCGCGGAATTCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCg                                                                                <  1:1492706/146‑1 (MQ=255)
                                                                 ggCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGat                                                                              >  2:1595363/1‑151 (MQ=255)
                                                                 ggCGTGCTGCTAAAAATTTATAAAACCCGACGGCGATAAAAGGGCAGGGCGTTTTGCTCTCCATCAGTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCg                                                                                <  1:260151/149‑1 (MQ=255)
                                                                 ggCGTGCTGCTAAAAATTAAAAAAACCCGCGGGCGATAAAAAGGCAGGGCGCTTTGCTGCCCATCAGTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCg                                                                                <  1:1242145/149‑1 (MQ=255)
                                                                        tgctTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTcc                                                                         >  2:125201/1‑149 (MQ=255)
                                                                         gctTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                                      >  1:692656/1‑151 (MQ=255)
                                                                         gctTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                                      >  1:1683858/1‑151 (MQ=255)
                                                                         gctTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                                      >  2:1948812/1‑151 (MQ=255)
                                                                            tCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                                      >  2:444966/1‑148 (MQ=255)
                                                                             caaaatttaaaaaaCCGCCCGGCAATATCAAGCAATGGCGAATTGCGTTCATACATTTCCACCAGGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAg                                                                       <  2:1926946/137‑1 (MQ=255)
                                                                              aaaGTTAAAAAACCCGGCGGGGATAACAAGGCGTGGCGTTTTGCTTCCCATCATTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAg                                                                       <  1:470917/145‑1 (MQ=255)
                                                                                 gTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATg                                                              >  1:848042/1‑151 (MQ=255)
                                                                                   taccAAAACCGACGGTAATATAAGGGTAACGGGACTCGCTCCACACAAGTCCCACCGCGGGCTTCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATgg                                                             <  2:496712/147‑1 (MQ=255)
                                                                                   aataaaacCCGAGGGAGATACCAGGGCAGGGCGACTTCTTGTCCGTCTTTCCCCCCACGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGa                                                            <  1:488980/147‑1 (MQ=255)
                                                                                    aacAAACCCGGGGGCGATACCAGGGCTGGGCGTTTTGTTGCCCATAGTTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATg                                                              <  1:222653/146‑1 (MQ=255)
                                                                                    aTCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAt                                                           >  1:1591550/1‑151 (MQ=255)
                                                                                       aaacCCGCCGGCGATATAAGGGCAGGGCGCTTTGCTGCCAATCATTCCCCCCACGGGAACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAtt                                                          <  1:950681/149‑1 (MQ=255)
                                                                                       aaacCCGACGGCAATATCAGGGCGGGGCGCCTTGCTCCCCATCATTCCCCCCAGGGGAAAACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAtt                                                          <  2:398395/149‑1 (MQ=255)
                                                                                       aaacCCGACGGCAATATAAAGGCAGGGCGTTTTGCTGTCCACCAGTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGa                                                            <  2:1179482/147‑1 (MQ=255)
                                                                                       aaaaCCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCg                                                        >  1:905852/1‑151 (MQ=255)
                                                                                          aacccgacgtaaatatCAAGGCTTGGCGACTTGCTTCCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCg                                                        <  1:192183/138‑1 (MQ=255)
                                                                                             ccgccggaaatacCAAGGCGGGGCGTCTTGTTCCCAATCAGTCCCCCCGCGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAgg                                                     <  1:530058/141‑1 (MQ=255)
                                                                                              gACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGt                                                  >  1:1961261/1‑151 (MQ=255)
                                                                                                  ggCGATATCAAGGCATGGCGATTTGCTGCCCATCAGTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTc                                                <  1:925176/150‑1 (MQ=255)
                                                                                                   gCGATAGCAAGGCATGGCGACTTGCGGCCCATCAGTCCCCCCGGGGGAAGCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCtt                                              <  1:1493336/151‑1 (MQ=255)
                                                                                                       tatCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTc                                          >  1:248035/1‑151 (MQ=255)
                                                                                                          cAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAt                                       >  2:1775601/1‑151 (MQ=255)
                                                                                                             ggcggggCGCCTTGCTGCCCACCTGTCCCCCCGGGGGAACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAt                                       <  1:1909419/143‑1 (MQ=255)
                                                                                                             ggCAGGGCGATTTGCTTTCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTc                                    <  2:1369377/151‑1 (MQ=255)
                                                                                                               cATGGCGTCTTGCTGCCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCg                                  <  2:52397/151‑1 (MQ=255)
                                                                                                                 tGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGt                                >  2:1852470/1‑151 (MQ=255)
                                                                                                                  ggCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGtt                               >  1:301271/1‑151 (MQ=255)
                                                                                                                     gACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtat                            >  1:1759763/1‑151 (MQ=255)
                                                                                                                      aCTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtat                            >  1:973604/1‑150 (MQ=255)
                                                                                                                       cTTGCTGTCCATCAACTCCACCACCGGGATCCCCGCCACTTCAATCATCTTTAAGGTGCGTGGCGTGTGAGTACGTTCGGTGAGGATCAGTCCATCGATATTCCATGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtata                           >  1:1939524/1‑150 (MQ=255)
                                                                                                                        ttGCTGTCCACCATTCCCCCCACGGGAACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTta                             <  2:272376/147‑1 (MQ=255)
                                                                                                                         tgctgcCCATCATTTCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATaa                          <  2:686648/149‑1 (MQ=255)
                                                                                                                          ttctttcccaccgTTCCCCCCCGGGAATCCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtat                            <  2:59519/136‑1 (MQ=255)
                                                                                                                          gCTGCCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATaa                          <  2:905868/148‑1 (MQ=255)
                                                                                                                           cTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGt                      >  2:773154/1‑151 (MQ=255)
                                                                                                                           cTGCCCATCTGTCCCACCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAcc                        <  2:1102626/149‑1 (MQ=255)
                                                                                                                            tccccATCATTTCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAc                         <  2:494141/144‑1 (MQ=255)
                                                                                                                                cATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGACGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGCTTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTgc                 >  1:986168/1‑151 (MQ=255)
                                                                                                                                    cAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGcc              >  1:1311438/1‑151 (MQ=255)
                                                                                                                                     cgttcccaccgCGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTgcgc               <  2:255938/145‑1 (MQ=255)
                                                                                                                                         ccccccagGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGt                      <  1:1694509/135‑1 (MQ=255)
                                                                                                                                          ccaccaCGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCAt         >  2:1186137/1‑151 (MQ=255)
                                                                                                                                           ccccagGGGATTACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATg        <  2:570054/149‑1 (MQ=255)
                                                                                                                                           caccacGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATg        >  1:1647603/1‑151 (MQ=255)
                                                                                                                                             ccacGAGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCa             <  1:1067583/144‑1 (MQ=255)
                                                                                                                                              cacGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTc     >  1:1695380/1‑151 (MQ=255)
                                                                                                                                               agggaaaTCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTc     <  1:1313626/145‑1 (MQ=255)
                                                                                                                                                 gggAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCt    <  1:869647/149‑1 (MQ=255)
                                                                                                                                                 ggaaaTCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATg        <  2:1194368/142‑1 (MQ=255)
                                                                                                                                                 ggaaaTCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGa  <  2:1728313/148‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                
GTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGA  >  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: