Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,452,022 G→A P89S (CCG→TCG)  rplD ← 50S ribosomal subunit protein L4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,452,0220GA94.7% 49.3 / ‑3.9 19coding (265/606 nt)rplD50S ribosomal subunit protein L4
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/9);  minor base . (1/0);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CGTACGTCAACCTTGTGCAGGTTGCGCGCAGCCAGGAACAGGTTTTCGTCCAGCTCACCGGTGATGATCAGCACATCTTCCAGAGCCATGTCTTTCAGTTTCTGTGCCAGCAGCTTAGTTTTCGGCGCTTCTACAGAGAACTTCTCGACAACGATCAGACGATCCTGACGTACCAGTTCGGACAGGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGAC‑GAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCT  >  NC_000913/3451778‑3452259
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
cgTACGTCAACCTTGTGCAGGTTGCGCGCAGCCAGGAACAGGTTTTCGTCCAGCTCACCGGTGATGATCAGCACATCTTCCAGAGCCATGTCTTTCAGTTTCTGTGCCAGCAGCTTAGTTTTCGGCGCTTCTACAGAGAACTTCTCGACAACTATCAGACGATCCTGACGTACCAGTTCGGACAGGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGAAC‑Gagc                                                                                                                                                                                                                                         <  1:415235/251‑1 (MQ=255)
    cGTCAACCTTGTGCAGGTTGCGCGCAGCCAGGAACAGGTTTTCGTCCAGCTCACAGGTGATGATCAGCACATCTTCCAGAGCCATGTCTTTCAGTTTCTGTGCCAGCAGCTTAGTTTTCGGCGCTTCTACAGAGAACTTCTCGACAACGATCAGACGATCCTGACGTACCAGTTCGGACAGGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGAAC‑Gag                                                                                                                                                                                                                                          >  1:202428/1‑246 (MQ=255)
            tggtgcagttggggggcagcaAGGAACAGGTTT‑CGTCCAGCTCACCGGTGATGATCAGCACATCTTCCAGAGCCATGTCTTTCAGTTTCTGTGCCAGCAGCTTAGTTTTCGGCGCTTCTACAGAGAACTTCTCGACAACGATCAGACGATCCTGACGTACCAGTTCGGACAGGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGAAC‑GAGCAGCAAAGGTCan                                                                                                                                                                                                                             <  1:266248/235‑2 (MQ=255)
                                                                                                                gcTTAGTTTTCTGCGCTTCTACAGAGAACTTCTCGACAACGATCAGACGATCCTGACGTACCAGTTCGGACAGGATGCTTTTCCGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGAAC‑Gcgcagc                                                                                                                                                                                                                                      >  1:77564/1‑142 (MQ=255)
                                                                                                                     gTTTTCGGCGCTTCTACAGAGAACTTCTCGACAACGATCAGACGATCCTGACGTACCAGTTCGGACAGGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGAAC‑GAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGt                                                                                                                    <  1:175795/251‑1 (MQ=255)
                                                                                                                               cTTCTACAGAGAACTTCTCGACAACGATCAGACGATCCTGAGGTACGAGTTCGGACCGGCTGCTTGTCAGAGCGCCGCGGGACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCG‑ACAGAGCAGCACAGGACACGCCACCAGGACGTCAGATCGGGCTCTTGATAGAAACAGAACGCGCACGGCCGGTGCCTTTCTGGAGCCAAGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGc                                                                                                           >  1:432085/1‑250 (MQ=255)
                                                                                                                                        cnttnncgatntcAACGATCAGACGATCCTGACGTNCCAGTTCGGACAGGATGCTTTTCAGCGCGCCNCNGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGAAC‑GAGCAGCNAAGGTCACGCNACCNNAACGNCNGNTCGGGCTCTTGATAGAACCAGAACGCGCANGGCCGGTGCCTTTCTGGCGCCACNNTTNNNTANNGGAAc                                                                                                                                       <  2:333074/201‑1 (MQ=255)
                                                                                                                                            cTTCTCGACAACGATCAGACGATCCTGACGTACCCGTTCGGACAGGATGCTTTTCCGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGAAC‑GAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCtt                                                                                                                                                                                                  >  1:158257/1‑150 (MQ=255)
                                                                                                                                                    cAACGATCAGACGATCCTGACGTACCAGTTCGGACAGGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGAAC‑GAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAAc                                                                                                                                       >  1:333074/1‑201 (MQ=255)
                                                                                                                                                                                          ntGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGAAC‑GCGCAGCAAAGGTCACGCAACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCAGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCa                                                                                                                            >  1:85490/2‑174 (MQ=255)
                                                                                                                                                                                             cTTTTCAGCGCGCCGCGGTACATCTTNTTGTTANNTTNTTGACTGTGGTCCTGCGAAN‑GANNNGNNNANGTCACNCCNNCAGAACGCCNGATCGGGCTCTTGATAGAACCAGAACGCGCACGGTCGGTGCCTTTCTGGCGcc                                                                                                                                                         >  2:106688/1‑142 (MQ=255)
                                                                                                                                                                                             cTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACGTTTTGACTGTGGTCCTGCGAAC‑GAGCAGTAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGTCGGTGCCTTTCTTGCGcc                                                                                                                                                         <  1:106688/142‑1 (MQ=255)
                                                                                                                                                                                                              gTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGAAC‑GAGCAGCAAAGGTCACGCAACCAGCAAGCCAGATCGGGCTCTTGATAGAACCAGAAAGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTAACGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCATAAGcaacaacaa                                                                  >  1:296133/1‑212 (MQ=255)
                                                                                                                                                                                                                  aTCTTCTTGTTAACTTTCTGACGGTGGTCCTGCGAAC‑GAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAAc                                                                                                                                                                                          <  1:205812/88‑1 (MQ=255)
                                                                                                                                                                                                                                tttttGACTGTGGTCCTGCGAAC‑GAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGgcgc                                                                                                                                                          <  1:98132/106‑1 (MQ=255)
                                                                                                                                                                                                                                       ntgtgGTCCTGCGAAC‑GAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTAAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGct  >  1:273732/2‑251 (MQ=255)
                                                                                                                                                                                                                                               cTGCGAAC‑GAGCNGNAAAGGGCACGCCACCAGAACGCCAGATCGTGCTCTTGATAGAACCAGANCGCGCACGGCNGGTNNCTTTNTNGNGCCACAGTTTTGTGCCGGAACCAGTTACTNCAGCACGAGTCTTCTGAGCACGANNGCNNNGANNAGCAc                                                                                       <  2:337950/158‑1 (MQ=255)
                                                                                                                                                                                                                                               cTGCGAAC‑GAGCAGCAAAGGTCACGCAACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGCACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCAc                                                                                       >  1:337950/1‑158 (MQ=255)
                                                                                                                                                                                                                                                tGCGAAC‑GAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCg                                                                                                                                           <  1:29183/105‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
CGTACGTCAACCTTGTGCAGGTTGCGCGCAGCCAGGAACAGGTTTTCGTCCAGCTCACCGGTGATGATCAGCACATCTTCCAGAGCCATGTCTTTCAGTTTCTGTGCCAGCAGCTTAGTTTTCGGCGCTTCTACAGAGAACTTCTCGACAACGATCAGACGATCCTGACGTACCAGTTCGGACAGGATGCTTTTCAGCGCGCCGCGGTACATCTTCTTGTTAACTTTTTGACTGTGGTCCTGCGGAC‑GAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCT  >  NC_000913/3451778‑3452259

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