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

GGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTT  >  NC_000913/4067600‑4067890
                                                                                                                                                    |                                                                                                                                                 
ggCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGggt                                                                                                                                               <  2:25225/150‑1 (MQ=255)
    aaaTGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGgtgt                                                                                                                                          <  2:997610/151‑1 (MQ=255)
    aaaTGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGgtgt                                                                                                                                          <  2:202415/151‑1 (MQ=255)
    aaaTGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGgtgt                                                                                                                                          <  1:425301/151‑1 (MQ=255)
    aaaTGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTGCTGCGCTTCATCTCAAACAGGGTGgtgt                                                                                                                                          <  1:100975/150‑1 (MQ=255)
           cGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGc                                                                                                                                    >  2:776086/1‑150 (MQ=255)
                    cTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCg                                                                                                                            >  2:978421/1‑149 (MQ=255)
                            ccgaaaaccTGCTAGTTGTCGCCAGGCGGTTTCCCTTGGGGGGGCCGCATCATGGGCGGGAGGGCGCTAAAATCGCCCCAGCGCAGCAGCAGCTCTTTGCGCCCTTTCATCTAAACCAGGGTGGTGTACCCGCCAATGTCGCTGTGGTGc                                                                                                                   <  2:801576/142‑1 (MQ=255)
                             cgcaaaccgccgtgTTGGCGCCGGGACGTTTCCCTTCGGGGGTGCCCCTCCGCGGCGAGAAGGCGCTAAAATCGCCCCGGCGCAGCCGCACCTTTTGGCTGGGTTTCATCAAAAACGGGGTGGTGAAACCGCCAATGTCGCTGTGGTGCAg                                                                                                                 <  2:877497/138‑1 (MQ=255)
                                aaaccTGCCTGTTGGCGCCAGGACGGTTCCCTTGGTGGGTGCGCATCATGGGCGTGAAGGCGCTGAATTCGCACCAGCGCACCAGCCGCTCTTTGCTCCGTTTCCCCCCAAAAGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGcc                                                                                                              <  1:851147/147‑1 (MQ=255)
                                aaaaCGGCCGGTTGCCGCCGGGACGGTTCCTTGGGGGGGGGCGCCTCAGCGGGGGAAGGGCGCTAAAATCCCACCAGCGCAGCAGCAGCTTTTTGCTGCGCTTCATCTAAAACGGGGTGGTGTACCCGCCAATGTCGCTGTGGTGCAGGcc                                                                                                              <  2:396370/150‑1 (MQ=255)
                                     tgccgtttGCCGCCGGGCCGGTTCCCTCGGGGGGTGCGCATCATCGGGGTAAGGGCGCGAAAACCGCACCGGCGCAGCAGCAGCTTTTTGCTGCGCTTCATCCCAAACGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc                                                                                                         <  2:81527/145‑1 (MQ=255)
                                     cgcagttggccGCCGGGACGGTTACCTTCGTGGGTCCGCATCTTCGGCGGAAGGGCGCTAAAACCGCACCGGCGCAGCAGCAGTTTTTGGCGGCGTTTCATCTAAAACGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc                                                                                                         <  1:369837/141‑1 (MQ=255)
                                      gccgttggccGCAAGGACGGTTCCCTGGGGGGGTGCGCATCATCGGGGTGAGGGCGCTAAAGCCGCCCCAGCGGAGCAGCGGTTTTTTCCTGCGCTTACTCTCAAACGGGGGGGTATAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc                                                                                                         <  1:210698/141‑1 (MQ=255)
                                      gCCGGTTGTCGCGGGGAGGGTTCCCTGGGTGGGTGCGCACCAGCGGCGTGAGGGCGCTAAAATCGCCCGGGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACGGGGGGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc                                                                                                         <  1:408215/150‑1 (MQ=255)
                                      gCCATTTGTCGCAAGGGCGGTTACCTTGGGGGGGGGGCATCAGCGGGGGGAGGCCGTGAAAAGGACAACAGCGCAGCAGCACCTCTTGGCTCCGTTTCACCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCAGGTc                                                                                                         <  2:913329/150‑1 (MQ=255)
                                      cccgttgggcGCCGGGACGTTTCCCTTCGGGGGTGCGCATCACCGGGGTAAGGGCGTTGAAACCGCACCAGGGCGGCAGCGGCTTTTTGCGGGGCTCCATCCAAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc                                                                                                         <  1:648686/141‑1 (MQ=255)
                                           ttGTCGCCGGGGCGGTTCCCTGCGGGGGTGCGCATCCTCGGCGGGAAGGCGCTAAAATCGCCCCGGCGCAGCAGCAGCTTTTGGCTGCGCTACATCCAAAACGGGGTGGGTAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTc                                                                                                    <  1:547182/150‑1 (MQ=255)
                                                cccAGGCCGTTTCCCTTGGGGGGTGCGCATCTTCGGCGTGAAGGCGCTGAAACCGCACCGGCGCAGCAGCAGCTTTTTGCTGCGTTTCATCAAAAACGGGGGGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGcc                                                                                              <  1:474115/150‑1 (MQ=255)
                                                 ccAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCa                                                                                             >  1:1033525/1‑151 (MQ=255)
                                                       cgtttccctcggTGGGTGCGCATCACGGGCGGGAAGGCGCGAAAACCGCACCGGCGCAGCAGCACTTTTTTGCTGCGCTTCATCTCAAAGGGGGGGGTGAACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGAc                                                                                        <  1:439280/139‑1 (MQ=255)
                                                                                   gCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGggt                                                                                                                                               <  2:118995/67‑1 (MQ=255)
                                                                                   gCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGggt                                                                                                                                               >  1:118977/1‑67 (MQ=255)
                                                                                    ggTGAAGGGCGCTAAAATGCACCAGGCCCAGCGGCCGCTTTTTGCTGCGTTTCATCTAAAACGGGGTGGTGAACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGc                                                            <  1:766105/149‑1 (MQ=255)
                                                                                    cGTGAGGGCGCTGAAATCGCACGGGCGCGGCAGCGGCTTTTTGCTGCGCTTCATCCAAAACGGGGTGGTGTACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGcc                                                           <  1:195716/150‑1 (MQ=255)
                                                                                             cgttaaaaCCGCCCCGGCGCACCAGCCGCTTTTTGCGGCGCTTCTCCCAAAACGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCgag                                                  <  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)
                                                                                                                                                    |                                                                                                                                                 
GGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTT  >  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: