Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,609,336 (C)7→8 intergenic (+13/‑78) yjjG → / → prfC dUMP phosphatase/peptide chain release factor RF‑3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,609,3291.C92.9% 41.3 / ‑2.8 14intergenic (+6/‑85)yjjG/prfCdUMP phosphatase/peptide chain release factor RF‑3
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (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 = 8.39e-01

TCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTG  >  NC_000913/4609242‑4609443
                                                                                         |                                                                                                                  
tCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcg                                                                  >  1:95667/1‑139 (MQ=255)
                    tCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTCCTc                                              >  1:801738/1‑139 (MQ=255)
                     cgcgCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCa                                             >  1:17023/1‑139 (MQ=255)
                     cgcgCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTc                                                                                                                   >  1:959282/1‑68 (MQ=255)
                     cgcgCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTc                                                                                                                   <  2:959282/68‑1 (MQ=255)
                           cACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga                                       <  2:17023/139‑1 (MQ=255)
                                      tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAATAATTATGa                            >  2:689379/1‑139 (MQ=255)
                                      tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGa                            >  1:1024654/1‑139 (MQ=255)
                                      tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGa                            >  1:479214/1‑139 (MQ=255)
                                      tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGa                            >  1:64707/1‑139 (MQ=255)
                                        tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                          <  2:64707/139‑1 (MQ=255)
                                               tGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCg                                                     >  1:920488/1‑105 (MQ=255)
                                               tGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCg                                                     <  2:920488/105‑1 (MQ=255)
                                                                cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATg                                                        <  1:243198/85‑1 (MQ=255)
                                                                cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATg                                                        >  2:243198/1‑85 (MQ=255)
                                                                cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTg  >  1:550608/1‑139 (MQ=255)
                                                                                         |                                                                                                                  
TCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTG  >  NC_000913/4609242‑4609443

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