Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,795,562 C→T D795N (GAT→AAT)  pheT ← phenylalanine‑‑tRNA ligase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,795,5620CT92.0% 65.7 / ‑0.3 25D795N (GAT→AAT) pheTphenylalanine‑‑tRNA ligase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (17/6);  total (18/7)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

GGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTCGGATAAAATATCCGCTGCGGGAACGTTTTCTGCG  >  NC_000913/1795338‑1795803
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
ggAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCtt                                                                                                                                                                                                                         >  1:546636/1‑251 (MQ=255)
       ccTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCAAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGAGCCATAGGTTCAATTCCTCATTGATGCCTGGAATCGCTCTTTTAATGc                                                                                                                                                                                                                  >  2:419447/1‑251 (MQ=255)
                                      ttGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAAtc                                                                                                                                                                                   >  1:430352/1‑251 (MQ=255)
                                                     cAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAgt                                                                                                                                                                    >  2:410542/1‑251 (MQ=255)
                                                           gAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACAtt                                                                                                                                                                                                       <  2:69759/210‑1 (MQ=255)
                                                           gAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACAtt                                                                                                                                                                                                       >  1:69759/1‑210 (MQ=255)
                                                                tACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTATTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAACCAGATATTCTGCCATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACCCATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGt                                                                                                                                                         >  2:456827/1‑251 (MQ=255)
                                                                        ccAGAGAGTTTCACCTGGTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTATGTAAGCGCCATAGGTTCAATTCCTCAAGGATGCCTGGAATCGATCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGTTATCTTGc                                                                                                                                                  >  2:157256/1‑250 (MQ=255)
                                                                                                          cgcgACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTc                                                                                                               >  2:266431/1‑251 (MQ=255)
                                                                                                                 gATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAac                                                                                                         >  1:53481/1‑250 (MQ=255)
                                                                                                                    ctctTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTg                                                                                                                                                                                                     <  1:186584/155‑1 (MQ=255)
                                                                                                                    ctctTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTg                                                                                                                                                                                                     >  2:186584/1‑155 (MQ=255)
                                                                                                                          gAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACcg                                                                                               >  1:236267/1‑251 (MQ=255)
                                                                                                                             aaaCAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCgg                                                                                            >  1:299866/1‑251 (MQ=255)
                                                                                                                                                     aTCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCCTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCa                                                                    <  2:210967/251‑1 (MQ=255)
                                                                                                                                                        ccGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACt                                                                 >  1:343023/1‑251 (MQ=255)
                                                                                                                                                                 aaGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGAttt                                                        <  1:20846/251‑1 (MQ=255)
                                                                                                                                                                          cTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTAAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGGATCTTGCAGGATCAGGCTTATAGAGAGGCTCTTATACCCCTCCTCAACACCCTTACAGCGGTCAACGTCAAATAAGTTTAAGCCCAATAACTGATTTACGCCAACt                                               >  2:436320/1‑251 (MQ=255)
                                                                                                                                                                            taacAAACAGATATCCTGACATTTCATCGTTTGTAAGCGACATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTGTTAATGCCTCTACACACTTGGCGACGGTAGCGGCTATCTCCTCGTCTGCGAGTGTACGGCTTGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTTCGcccaca                                               <  2:133350/246‑5 (MQ=255)
                                                                                                                                                                                aaaCAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCtt                                          >  1:319417/1‑250 (MQ=255)
                                                                                                                                                                                   cAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTCCCTGATTTACGCCACCTTTCTTACAt                                      >  2:570944/1‑251 (MQ=255)
                                                                                                                                                                                             gACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTAAATTCGGATaaa                             >  1:541605/1‑250 (MQ=255)
                                                                                                                                                                                                        ttttGTAAGCGCCATAGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTCGGATAAAATATCCGCTGCg                 >  1:355170/1‑251 (MQ=255)
                                                                                                                                                                                                                 cgcCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTCGGATAAAATATCCGCTGCGGGAACGttt        <  1:266431/251‑1 (MQ=255)
                                                                                                                                                                                                                       aGGTTCAATTCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTCGGATAAAATATCCGCTGCGGGAACGTTTTCTGCg  <  1:410542/251‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
GGAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCTTCTTCGAGTGTACGGCTGGTATCTTGCAGGATCAGGCTTATGGCGAGGCTCTTATACCCCTCCGCAACACCCTTACCGCGGTACACGTCAAATAAGTTTACGCCAACTACCTGATTTACGCCAACTTTCTTACATTCGGATAAAATATCCGCTGCGGGAACGTTTTCTGCG  >  NC_000913/1795338‑1795803

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