Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,067,7450AG47.2% 33.8 / 10.5 36L511P (CTG→CCG) yihQputative sulpholipid alpha‑glucosidase; alpha‑glucosyl fluoride glucosidase
Reads supporting (aligned to +/- strand):  ref base A (10/9);  new base G (0/17);  total (10/26)
Fisher's exact test for biased strand distribution p-value = 4.40e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTT  >  NC_000913/4067600‑4067890
                                                                                                                                                    |                                                                                                                                                 
ggCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGggt                                                                                                                                               <  2:25225/150‑1 (MQ=255)
    aaaTGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGgtgt                                                                                                                                          <  2:997610/151‑1 (MQ=255)
    aaaTGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGgtgt                                                                                                                                          <  2:202415/151‑1 (MQ=255)
    aaaTGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGgtgt                                                                                                                                          <  1:425301/151‑1 (MQ=255)
    aaaTGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTC‑TTGCTGCGCTTCATCTCAAACAGGGTGgtgt                                                                                                                                          <  1:100975/150‑1 (MQ=255)
           cGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGc                                                                                                                                    >  2:776086/1‑150 (MQ=255)
                    cTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCg                                                                                                                            >  2:978421/1‑149 (MQ=255)
                            ccgaaaaccTGCTAGTTGTCGCCAGGCGGTTTCCCTTGGGGGGGCCGCATCATGGGCGGGAG‑GGCGCTAAAATCG‑CCCCA‑GCGCAGCAGCAGCTCTTTGCGCCCTTTCATCTAAACCAGGGTGGTGTACCCGCCAATGTCGCTGTGGTGc                                                                                                                   <  2:801576/142‑1 (MQ=255)
                             cgcaaaccgccgtgTTGGCGCCGGGACGTTTCCCTTCGGGGGTGCCCCTCCGCGGCGAGAA‑GGCGCTAAAATCG‑CCCCG‑GCGCAGCCGCACCTTTTGGCTGGGTTTCATCAAAAACGGGGTGGTGAAACCGCCAATGTCGCTGTGGTGCAg                                                                                                                 <  2:877497/138‑1 (MQ=255)
                                aaaccTGCCTGTTGGCGCCAGGACGGTTCCCTTGGTGGGTGCGCATCATGGGCGTGAA‑GGCGCTGAATTCG‑CACCA‑GCGCACCAGCCGCTCTTTGCTCCGTTTCCCCCCAAAAGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGcc                                                                                                              <  1:851147/147‑1 (MQ=255)
                                aaaaCGGCCGGTTGCCGCCGGGACGGTTCCTTGGGGGGGGGCGCCTCAGCGGGGGAAG‑GGCGCTAAAATCC‑CACCA‑GCGCAGCAGCAGCTTTTTGCTGCGCTTCATCTAAAACGGGGTGGTGTACCCGCCAATGTCGCTGTGGTGCAGGcc                                                                                                              <  2:396370/150‑1 (MQ=255)
                                     tgccgtttGCCGCCGGGCCGGTTCCCTCGGGGGGTGCGCATCATCGGGGTAAG‑GGCGCGAAAACCG‑CACCG‑GCGCAGCAGCAGCTTTTTGCTGCGCTTCATCCCAAACGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc                                                                                                         <  2:81527/145‑1 (MQ=255)
                                     cgcagttggccGCCGGGACGGTTACCTTCGTGGGTCCGCATCTTCGGCGGAAG‑GGCGCTAAAACCG‑CACCG‑GCGCAGCAGCAGTTTTTGGCGGCGTTTCATCTAAAACGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc                                                                                                         <  1:369837/141‑1 (MQ=255)
                                      gccgttggccGCAAGGACGGTTCCCTGGGGGGGTGCGCATCATCGGGGTGAG‑GGCGCTAAAGCCG‑CCCCA‑GCGGAGCAGCGGTTTTTTCCTGCGCTTACTCTCAAACGGGGGGGTATAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc                                                                                                         <  1:210698/141‑1 (MQ=255)
                                      gCCGGTTGTCGCGGGGAGGGTTCCCTGGGTGGGTGCGCACCAGCGGCGTGAG‑GGCGCTAAAATCG‑CCCGG‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACGGGGGGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc                                                                                                         <  1:408215/150‑1 (MQ=255)
                                      gCCATTTGTCGCAAGGGCGGTTACCTTGGGGGGGGGGCATCAGCGG‑GGGGA‑GGCCGTGAAAAGGACAACA‑GCGCAGCAGCACCTCTTGGCTCCGTTTCACCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCAGGTc                                                                                                         <  2:913329/150‑1 (MQ=255)
                                      cccgttgggcGCCGGGACGTTTCCCTTCGGGGGTGCGCATCACCGGGGTAAG‑GGCGTTGAAACCG‑CACCA‑GGGCGGCAGCGGCTTTTTGCGGGGCTCCATCCAAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc                                                                                                         <  1:648686/141‑1 (MQ=255)
                                           ttGTCGCCGGGGCGGTTCCCTGCGGGGGTGCGCATCCTCGGCGGGAA‑GGCGCTAAAATCG‑CCCCG‑GCGCAGCAGCAGCTTTTGGCTGCGCTACATCCAAAACGGGGTGGGTAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTc                                                                                                    <  1:547182/150‑1 (MQ=255)
                                                cccAGGCCGTTTCCCTTGGGGGGTGCGCATCTTCGGCGTGAA‑GGCGCTGAAACCG‑CACCG‑GCGCAGCAGCAGCTTTTTGCTGCGTTTCATCAAAAACGGGGGGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGcc                                                                                              <  1:474115/150‑1 (MQ=255)
                                                 ccAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCa                                                                                             >  1:1033525/1‑151 (MQ=255)
                                                       cgtttccctcggTGGGTGCGCATCACGGGCGGGAA‑GGCGCGAAAACCG‑CACCG‑GCGCAGCAGCACTTTTTTGCTGCGCTTCATCTCAAAGGGGGGGGTGAACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGAc                                                                                        <  1:439280/139‑1 (MQ=255)
                                                                                   gCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGggt                                                                                                                                               <  2:118995/67‑1 (MQ=255)
                                                                                   gCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGggt                                                                                                                                               >  1:118977/1‑67 (MQ=255)
                                                                                    ggTGAAGGGCGCTAAAAT‑G‑CACCAGGCCCAGCGGCCGCTTTTTGCTGCGTTTCATCTAAAACGGGGTGGTGAACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGc                                                            <  1:766105/149‑1 (MQ=255)
                                                                                    cGTGAG‑GGCGCTGAAATCG‑CACGG‑GCGCGGCAGCGGCTTTTTGCTGCGCTTCATCCAAAACGGGGTGGTGTACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGcc                                                           <  1:195716/150‑1 (MQ=255)
                                                                                             cgttaaaaCCG‑CCCCG‑GCGCACCAGCCGCTTTTTGCGGCGCTTCTCCCAAAACGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCgag                                                  <  2:174940/146‑1 (MQ=255)
                                                                                                   aaaaccgcccaggcgcAGAGGCGGCTTTTTGCTGCTTTTCATCTAAAAGGGGGTGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACt                                               <  2:796662/138‑1 (MQ=255)
                                                                                                                cggagccgcacttttttGCTGCCCTTCCCCCAAAACGGGGGGGGTTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTc                                  <  2:26559/136‑1 (MQ=255)
                                                                                                                     cagcagcTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcg                            >  2:859950/1‑151 (MQ=255)
                                                                                                                          gcTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgcccg                        >  1:857304/1‑150 (MQ=255)
                                                                                                                                     cgcTTCATCACAAACGGGGTGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCAtgg            <  2:158992/151‑1 (MQ=255)
                                                                                                                                      ctttaaTCCCAAAGGGGGTGGTGTACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCAtggt           <  2:256136/146‑1 (MQ=255)
                                                                                                                                           aTCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGt      >  2:29895/1‑151 (MQ=255)
                                                                                                                                           aTCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGt      >  1:778701/1‑151 (MQ=255)
                                                                                                                                           aTCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCCGGTCGCCCGCCCACATCATGGGGGAGt      >  2:1052554/1‑151 (MQ=255)
                                                                                                                                               cAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGCCGAGACTCCAGTCGACGTTCTGGGCCCCCGCCCACATCATGGCGGAGTAttt  >  2:510083/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTT  >  NC_000913/4067600‑4067890

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