Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,433,900 C→T L71L (CTA→TTA)  def → peptide deformylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,433,9000CT92.9% 33.4 / ‑6.6 14L71L (CTA→TTA) defpeptide deformylase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (9/4);  total (10/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGC  >  NC_000913/3433662‑3434109
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                 
tAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCAGCTTCGCCAAGATGCTAAACCGGTAGAAGAAGCGAAAGCAGACCTTCAGCGTATCGACGATGAAATGTTCGAGACGATGTACGCCGAAGAAGGTATTGGCCTGGCGGCAAACCAGGATGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGAGACGGACGGCTCGTGTTaatc                                                                                                                                                                                                        >  1:331457/1‑250 (MQ=255)
                                                                  gCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGaac                                                                                                                                     >  1:132106/1‑251 (MQ=255)
                                                                      cGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTTGGAAAACCGTGACGAACGGTTAgtgt                                                                                                                                                                                                             >  1:354911/1‑175 (MQ=255)
                                                                                cTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATACGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGc                                                                                                                                                                    >  1:107318/1‑206 (MQ=255)
                                                                                                                                                 tACGCAGAAGAAGGTATTAGCCTGGCGGCAACACCGGGTGATATCCATCACCGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTCGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGGGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTACACCATTTg                                                                >  1:117541/1‑241 (MQ=255)
                                                                                                                                                                   ggCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGNTTGNNNNNCgnnnnnnnnnn                                                                                                                                     >  2:213078/1‑144 (MQ=255)
                                                                                                                                                                        ggcggcAACCCAGGTTGATATCCATCAACGTATCCTTGTAACTGATGTTTCGGAAAACCGTGACGACCGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGCACACCGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGCCGGTAAACCATTTGCACTGGAAGCAGACGGTCTGTTAGCCATCTGt                                >  1:311409/1‑250 (MQ=255)
                                                                                                                                                                          cggcAACCCAGGTTGATATCCATCAACGTCTCATTGTTATTGATGTTTCGGAAACCCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTACACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTAtt                             >  1:138580/1‑251 (MQ=255)
                                                                                                                                                                                        gATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTCAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATgg                >  1:120618/1‑250 (MQ=255)
                                                                                                                                                                                        gATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAAACCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATggn               <  1:213078/251‑2 (MQ=255)
                                                                                                                                                                                                   aCGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCATAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTc                                              <  1:203221/209‑1 (MQ=255)
                                                                                                                                                                                                                   gATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGAc                                                  >  1:140882/1‑189 (MQ=255)
                                                                                                                                                                                                                           ggAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTGTTAGAAAAAAGCGGCGAAACAGGCAGTGAAGAAGGTTGCTTGTCGAGCCCTGAACAACGTTCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCACTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGc  <  1:268768/229‑1 (MQ=255)
                                                                                                                                                                                                                             aaaaCCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGt                                                                                                                                                        <  1:118087/77‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                 
TAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGC  >  NC_000913/3433662‑3434109

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: