Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,077,129 G→A E39K (GAA→AAA)  dtd → D‑tyr‑tRNA(Tyr) deacylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,077,1290GA95.7% 69.1 / ‑3.4 23E39K (GAA→AAA) dtdD‑tyr‑tRNA(Tyr) deacylase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (14/8);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

ACGAACCATGATTGCATTAATTCAACGCGTAACCCGTGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTC  >  NC_000913/4077008‑4077257
                                                                                                                         |                                                                                                                                
acgaACCATGATTGCATTAATTCAACGCGTAACCCGTGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGt                                                                                                                 <  1:161047/139‑1 (MQ=255)
 cgaACCATGATTGCATTAATTCAACGCGTAACCCGTGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTc                                                                                                                >  1:391588/1‑139 (MQ=255)
           ttGCATTAATTCAACGCGTAACCCGTGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCg                                                                                                      >  2:430191/1‑139 (MQ=255)
              cATTAATTCAACGCGTAACCCGTGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCgtgt                                                                                                   >  1:409077/1‑139 (MQ=255)
                       aaCGCGTAACCCGTGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACAAAc                                                                                                                               >  2:548004/1‑102 (MQ=255)
                       aaCGCGTAACCCGTGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACAAAc                                                                                                                               <  1:548004/102‑1 (MQ=255)
                              aaCCCGTGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATc                                                                                   >  1:107723/1‑139 (MQ=255)
                                                        gagaAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTAcc                                                                                        >  1:556490/1‑108 (MQ=255)
                                                        gagaAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTAcc                                                                                        <  2:556490/108‑1 (MQ=255)
                                                          gaAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGc                                                                                                       >  1:277686/1‑91 (MQ=255)
                                                          gaAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGc                                                                                                       <  2:277686/91‑1 (MQ=255)
                                                            aGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAAtctc                                                     >  2:200816/1‑139 (MQ=255)
                                                                 cGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCa                                                                <  1:467893/123‑1 (MQ=255)
                                                                 cGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCa                                                                >  2:467893/1‑123 (MQ=255)
                                                                 cGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGt                                                >  1:223097/1‑139 (MQ=255)
                                                                                                aTTGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGtt                 >  1:65241/1‑139 (MQ=255)
                                                                                                aTTGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGtt                 >  2:103966/1‑139 (MQ=255)
                                                                                                 ttGGGTGTCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGttt                <  1:430191/139‑1 (MQ=255)
                                                                                                        tCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCaa                                                   <  1:146042/97‑1 (MQ=255)
                                                                                                        tCGAAAAGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCaa                                                   >  2:146042/1‑97 (MQ=255)
                                                                                                           aaaaGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCt                                                                                                <  2:219870/49‑1 (MQ=255)
                                                                                                           aaaaGGATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCt                                                                                                >  1:219870/1‑49 (MQ=255)
                                                                                                               ggATGACGACAAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTc  <  2:65241/139‑1 (MQ=255)
                                                                                                                         |                                                                                                                                
ACGAACCATGATTGCATTAATTCAACGCGTAACCCGTGCCAGCGTCACCGTGGAGGGAGAAGTGACGGGCGAAATTGGCGCGGGACTTTTGGTGTTATTGGGTGTCGAAAAGGATGACGACGAACAGAAAGCAAACCGTCTGTGCGAGCGTGTGCTCGGCTACCGCATCTTTAGCGATGCCGAAGGCAAGATGAATCTCAACGTGCAACAGGCGGGCGGCAGTGTGCTGGTGGTTTCCCAGTTTACCCTC  >  NC_000913/4077008‑4077257

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