Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,769,928 C→G 15.4% S285R (AGC→AGG)  ydiK → UPF0118 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,769,9280CG15.4% 54.2 / 5.0 26S285R (AGC→AGG) ydiKUPF0118 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (13/9);  new base G (2/2);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

CGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGT  >  NC_000913/1769781‑1770073
                                                                                                                                                   |                                                                                                                                                 
cGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCg                                                                                                                                                  >  2:293870/1‑149 (MQ=255)
cGGCGTACCTTATGCAACGTTGCTAACGGTGGTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCg                                                                                                                                                  <  8:273109/149‑1 (MQ=255)
       ccTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGGGGTGTGGt                                                                                                                                           >  3:272814/1‑149 (MQ=255)
                 cGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTg                                                                                                                                 >  2:52485/1‑149 (MQ=255)
                              gtTAATGATCCTCTCCTGCCTTGTCCAGCTTGGACCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATcc                                                                                                                    >  6:351498/1‑149 (MQ=255)
                                  aTGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGcc                                                                                                                 >  2:181620/1‑148 (MQ=255)
                                  aTGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGcc                                                                                                                 <  1:52485/148‑1 (MQ=255)
                                   tGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGccc                                                                                                                >  7:51234/1‑148 (MQ=255)
                                                ccTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCg                                                                                                    <  1:293870/147‑1 (MQ=255)
                                                     tCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGt                                                                                             >  1:5153/1‑149 (MQ=255)
                                                               cccGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTAc                                                                                   >  4:190650/1‑149 (MQ=255)
                                                                          tACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGGGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCt                                                                         >  5:272221/1‑148 (MQ=255)
                                                                            cTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGGGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGa                                                                       <  6:272221/148‑1 (MQ=255)
                                                                              gATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCTTGGGTGCCGATTTACCGCTGATCCTGATtc                                                                    >  5:325065/1‑149 (MQ=255)
                                                                                                    tACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAAATCGCATGGGTGCCGATTTACCGCTCATCCTGATTCTCACTGGCGTTATTGGTGGttt                                              >  4:305995/1‑149 (MQ=255)
                                                                                                    tACTGGACTGGCGATACCACCTGGGGAACGGTATTGGTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACTTCATCAGCCCAATGTTAATTCGCATGGGTGCCGATATCACGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGttt                                              >  4:303658/1‑149 (MQ=255)
                                                                                                      cttgactggCGATACCACCTGGGGAACGGTTTTTTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGGGCCTATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGa                                            <  7:358647/146‑1 (MQ=255)
                                                                                                         gactggCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGGGTGGTTGGCACGCTGGATAACGGCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGAtt                                          <  4:317353/148‑1 (MQ=255)
                                                                                                                gATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCgg                                  >  4:223132/1‑149 (MQ=255)
                                                                                                                  tACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGga                                <  3:190650/149‑1 (MQ=255)
                                                                                                                      accTGGGGAACGGTATTGTTAGTGTGGAGGGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGgatgat                            <  4:272814/149‑1 (MQ=255)
                                                                                                                           gggAACGGTATTGTTAGTGTAGAGCGCTGTGGTTGGCAAACTAGATAACGTAATCCGCCCAATGTTAATTCGAATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTc                       >  2:136542/1‑149 (MQ=255)
                                                                                                                              aaCGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGGCATCCGCCCAATGGTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGt                    <  6:323340/149‑1 (MQ=255)
                                                                                                                                       gTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTc                                    >  1:18799/1‑124 (MQ=255)
                                                                                                                                       gTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTc                                    <  2:18799/124‑1 (MQ=255)
                                                                                                                                                gAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGt  <  8:51234/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGATTATCTGGCTCTACTGGACTGGCGATACCACCTGGGGAACGGTATTGTTAGTGTGGAGCGGTGTGGTTGGCACGCTGGATAACGTCATCCGCCCAATGTTAATTCGCATGGGTGCCGATTTACCGCTGATCCTGATTCTCTCTGGCGTTATTGGTGGTTTGATTGCTTTCGGGATGATCGGTCTGTTTATTGGTCCGGTTCTGT  >  NC_000913/1769781‑1770073

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