Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,111,7470GT53.3% ‑1.3 / 23.5 15N156K (AAC→AAAyiiQDUF1454 domain‑containing protein YiiQ
Reads supporting (aligned to +/- strand):  ref base G (6/1);  new base T (0/8);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 1.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCC  >  NC_000913/4111612‑4111878
                                                                                                                                        |                                                                                                                                   
ttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGtttt                                                                                                                                 >  1:311463/1‑140 (MQ=255)
                                     cccAGCTTAAGGGTTCCAAGAGCGAAGGTCACCCCCTTTTCGCGTTTGCCCGCGACAACATACCGCAGTGCCCCTTCTGTCTCGGTTTAATATCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                            <  2:21623/139‑1 (MQ=255)
                                               cggctaaaaCAGCGAGGGGCACCCCTTTTTCCCCGTGGCCCGCGAAAACATAACGCGTTGCCCCTTCCGTCTGGGTGAAATAACTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                                                                                  <  2:279399/133‑1 (MQ=255)
                                                  gTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTc                                                                               >  1:530498/1‑140 (MQ=255)
                                                   taaaaCAGGGAAGGTCACCCCTTTTGCGCCTTTGTCCGCGACAACATAACGCGGTGCCCTTTCTGTCTGGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                              <  2:311463/137‑1 (MQ=255)
                                                    tCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt                                                                             >  1:304114/1‑140 (MQ=255)
                                                              aagctcgCCCCTTTTCCGCGTTTGCCCGGGACAAAATACCGCAGTCCCCCTCCTGCCTGGGTGAAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                    <  2:304114/134‑1 (MQ=255)
                                                                 ccacccctTTTTCCCCGTTCCCCGCAACAACAAAACGCATTGCCCCTTCTTTCCGGTTGTAATACTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgt                                                                <  2:441972/136‑1 (MQ=255)
                                                                         ttttCCCCTTTCCCCGCAACAACATACCGCGGTCCCCTTTCTGTTTGGGTTTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGAt                                                        <  2:169807/140‑1 (MQ=255)
                                                                                  tttGTCCGCAAAAACATACCGCGGTCCCCTTTCTGTCTGGGTGTAATACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                <  2:465044/139‑1 (MQ=255)
                                                                                  tttGCCCGCAAAAAAATACCGCAGTCCCCTTTCTGTCTGGTTGTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTAATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                <  2:647520/139‑1 (MQ=255)
                                                                                    tccccccgAAAACATACCGCGGCCCCCCTTCTGTTTGGGTTTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGcc                                              <  2:92572/134‑1 (MQ=255)
                                                                                      cccGCAAAAACATAACGCAGTCCCCTTTCTGTCTGGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCAt                                            <  2:530498/139‑1 (MQ=255)
                                                                                                          tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                        >  2:540355/1‑140 (MQ=255)
                                                                                                                    gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt              >  1:244259/1‑140 (MQ=255)
                                                                                                                           tgtAATAACGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct       <  1:540355/140‑1 (MQ=255)
                                                                                                                               ataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc   >  1:428381/1‑140 (MQ=255)
                                                                                                                                tacttttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGccc  <  2:264741/135‑1 (MQ=255)
                                                                                                                                tacctttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGccc  <  1:218780/135‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                   
TTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCC  >  NC_000913/4111612‑4111878

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: