Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,840,9460C.22.2% 42.4 / 10.1 18coding (458/621 nt)ynjFCDP‑alcohol phosphatidyltransferase family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (6/8);  new base . (2/2);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTTGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAAACAACAACCAGCCGCCTGCCAGCGCATTTTG  >  NC_000913/1840818‑1841085
                                                                                                                                |                                                                                                                                           
cACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGAc                                                                                                                         <  6:164089/149‑1 (MQ=255)
        acCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCAGATAGTAAAAc                                                                                                                            <  5:1579/137‑1 (MQ=255)
        acCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCAGATAGTAAAAc                                                                                                                            >  6:1579/1‑137 (MQ=255)
                          aTCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTATGTGCATAACCg                                                                                                          <  5:114470/138‑1 (MQ=255)
                          aTCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTATGTGCATAACCg                                                                                                          >  6:114470/1‑138 (MQ=255)
                               gCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGa                                                                                                                                     <  3:35145/106‑1 (MQ=255)
                               gCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGa                                                                                                                                     >  4:35145/1‑106 (MQ=255)
                                        aCCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTtgctg                                                                                 >  5:92652/1‑149 (MQ=255)
                                           gAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCAGATAGTAAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTTGCTGCCAg                                                                             >  5:164297/1‑149 (MQ=255)
                                               atccacgcaatccacgcaAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCa                                                                                                                                       <  2:30876/88‑1 (MQ=255)
                                                                                      cccAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTTGCTGCCAGTGCAGCAAAGGCGAGAAAActg                                                       <  6:82452/129‑1 (MQ=255)
                                                                                      cccAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTTGCTGCCAGTGCAGCAAAGGCGAGAAAActg                                                       >  5:82452/1‑129 (MQ=255)
                                                                                             caAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTTGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAaacaa                             >  2:46848/1‑148 (MQ=255)
                                                                                                           gTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTTGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAAACAACAACCAGCCGCCTGc              <  6:92652/149‑1 (MQ=255)
                                                                                                            tttCTGTGCCTTCGGTTAATCCCCCAGATAGTAAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTTGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAAACAACAACCAGCCGCCTGCCa            <  6:164297/149‑1 (MQ=255)
                                                                                                                        cGGTTAATCCCCCCAGATAGTAAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTTGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAAACAACAACCAGCCGCCTGCCAGCGCATTTTg  <  1:46848/148‑1 (MQ=255)
                                                                                                                              aTCCCCCCAGATAGTAAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTtgct                                                                                  >  2:30395/1‑62 (MQ=255)
                                                                                                                              aTCCCCCCAGATAGTAAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTtgct                                                                                  <  1:30395/62‑1 (MQ=255)
                                                                                                                                |                                                                                                                                           
CACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCCCCAGATAGTAAAACGACTTATGTGCATAACCGGGGTTATCAATCTGATGTTTTGCTGCCAGTGCAGCAAAGGCGAGAAAACTGCTGCCCGTGCCGATAAACGCAAACAACAACCAGCCGCCTGCCAGCGCATTTTG  >  NC_000913/1840818‑1841085

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