Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,912,078 A→G S167P (TCC→CCC)  htpX ← zinc dependent endoprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,912,0780AG95.7% 73.5 / NA 23S167P (TCC→CCC) htpXzinc dependent endoprotease
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (12/10);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTC  >  NC_000913/1911943‑1912193
                                                                                                                                       |                                                                                                                   
tGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATa                                                                                                                 <  1:449456/140‑1 (MQ=255)
    cGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGa                                                                                                             >  1:485692/1‑140 (MQ=255)
               cAAACACCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAg                                                                                                     <  2:390601/137‑1 (MQ=255)
               cAAACACCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAg                                                                                                     >  1:390601/1‑137 (MQ=255)
                 aaCACCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGtt                                                                                                >  2:24835/1‑140 (MQ=255)
                    acCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTcac                                                                                             >  1:554225/1‑140 (MQ=255)
                          tCCAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCACGcc                                                                                       <  1:402295/140‑1 (MQ=255)
                          tCCAGAACCGTTGCAACAGCAATGTAGATCACCGGTTTGTCGTTGCTCCCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCACGcc                                                                                       <  2:178986/140‑1 (MQ=255)
                              gAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCACGCCCTGa                                                                                   <  1:24835/140‑1 (MQ=255)
                                    ttGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGc                                                                             >  1:365658/1‑140 (MQ=255)
                                             cAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCTATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTg                                                                    <  2:18816/140‑1 (MQ=255)
                                                            ggTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCATGCCCTGAGTCAGCGGCATGGTGACCATATCACCATTg                                                     >  1:168331/1‑140 (MQ=255)
                                                            ggTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCACGCCCTGACTCATCGTCATGGTCACCATATCACCATTg                                                     >  1:227121/1‑140 (MQ=255)
                                                            ggTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTg                                                     >  1:146719/1‑140 (MQ=255)
                                                               tGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCg                                                  >  1:121247/1‑140 (MQ=255)
                                                                                      tcaCGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGACATGGTGAGCATTTCAACATTGGCGGTGTGACTGATTTCGTGAGCAAt                           >  2:194400/1‑140 (MQ=255)
                                                                                           aTTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCg                      >  1:526737/1‑140 (MQ=255)
                                                                                                      tAAAACCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGAGGTGGCTGGTTTGGTGAGCAATGACCGCCTCGGCTTCa           >  1:182966/1‑140 (MQ=255)
                                                                                                         aaCCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCCCCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATcc        <  2:121247/140‑1 (MQ=255)
                                                                                                          aCCCGCGGCAAGCTGTGCCAGACTATGGGAAATAAAGCTAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCg       <  1:194400/140‑1 (MQ=255)
                                                                                                          aCCCGCGGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCCCCACGCCCTGAATCAGCGTCATGGTGACCATATCACCTTTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCg       <  2:485692/140‑1 (MQ=255)
                                                                                                             cgcgGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGc    <  2:661308/140‑1 (MQ=255)
                                                                                                               cgGCAAGCTGCGCCAGAATACGGGGAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTc  <  2:146719/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                   
TGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGCAAAGTAGATCAGCGGGTTGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTC  >  NC_000913/1911943‑1912193

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