Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 393 T→G R19R (CGT→CGG thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133930TG92.0% 50.9 / ‑0.8 25R19R (CGT→CGGthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (13/10);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

CGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGGTTCCCGCTGGCGCAATTGAAAACTTTCGTCGATC  >  NC_000913/152‑622
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                     
cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGat                                                                                                                                                                                                                              >  1:270520/1‑250 (MQ=255)
     cgacagaTACAAATTACAGAGTACCCAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAACAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGCAACGAAGTAACAACCATGCGAGTGTTGAAGTACGGAGGTACATCAGTGGCAAATGAAGAACGGTTTCTGGGGGTTGGCGATTTTCtt                                                                                                                                                                                                                        >  2:245490/2‑250 (MQ=255)
               aaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGAGCTTTTTTTTCGACCAAAGTTAACGGGGTAACAACCATGCGAGTTTTGAAGTTCTGCGGTCCATCAGTGGCTAATGCAAAACGTTTTCTGCGGGTTGCCGATCTTCTGGAAAGCAATg                                                                                                                                                                                                              >  1:137131/1‑251 (MQ=255)
                     cAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCcaggc                                                                                                                                                                                                        <  1:245490/251‑1 (MQ=255)
                         gTACACAACCTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCAAggcaggg                                                                                                                                                                                                    >  1:264046/1‑251 (MQ=255)
                                acaTCCAATAAACGCATTAGCCCCCCCATTCACACCACCATCACCATTACCACAGGTAACGGTGAGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg                                                                                                                                                                                             <  2:81191/251‑1 (MQ=255)
                                                       caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccccct                                                                                                                                                                       >  1:167112/1‑250 (MQ=255)
                                                            ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGCGGCCACCGTCCTCTCTGCCCCCGCCaa                                                                                                                                                                   >  1:95047/1‑250 (MQ=255)
                                                                  caaccatctccatTACCAAGGGAAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGACCGCCCATGACAGGGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTAACAACAATGCGAGTGTTTAAGTTCGGCGGTACATCAGTGGCAAAAGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGTAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAAGCacca                                                                                                                                                           <  2:108720/249‑1 (MQ=255)
                                                                   caacatcaccatTACCACAGGAAACGGTGCGAGCTAACGCGGAAAGGAAACACAGAAAAAAGCCCGCACCTGAAAGTGCGGGCTTTTAATTAGACCAAAGGTAAAACAGAAACAACGATGAGAGAGTTGAAGTTCGGAGGTACATCAATGGCACAAGACAGACGGTTTATGCGGGGATCCGATATTCTGACAATCACAGGCAGGCAGGGGCAGGAGGCCACCGTACATTCTGGCACCGGCAAAATCaccg                                                                                                                                                           >  2:375699/4‑249 (MQ=255)
                                                                            catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTAACGCGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGGCCTCTCTGCCCCCGCCCAAATCACCAACCACCtggt                                                                                                                                                 >  2:69040/1‑251 (MQ=255)
                                                                                   acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgatg                                                                                                                                          <  1:281161/251‑1 (MQ=255)
                                                                                            ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAAGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCTCCAAAATCACCAACCAACTGGTTGCGATGTATGaaaa                                                                                                                                  >  2:270049/1‑250 (MQ=255)
                                                                                                                           aaaaaGCCCCCACCTTACAGTGCGTGTTTGTTTTTCGACAAATGTTAATGGGGTAACAACAAGGCGAGGGTTGACGTTCGGCGTTACATCAGTGGCAAATGCAGAACGTTTTCGGGGGGTTGCCGATATTCTGAAAAGCAATGCCAGGCAGGGGGAGGTGGCCACCGTCCTCTCTTCCCCCGCCAAAATCACCAACCCCCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAata                                                                                                  <  2:326251/251‑1 (MQ=255)
                                                                                                                                        cTGACAGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCg                                                                                      >  1:122160/1‑250 (MQ=255)
                                                                                                                                          gACAGTGCGGGCTTTTTTTTCGACCAAAGGGAATGAGGTAACAACCATGCGAGTGTGGAAGTTCGGCGGTACACCAGTGGCAAATGCAGAACGTTGTCTGCGGGTTGCCGATATTCTGGAAAGCAAGGCCAGGCAGGGTCAGTTGGCACCCGTCCTCACGGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGa                                                                                     <  2:336001/249‑1 (MQ=255)
                                                                                                                                             aGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCttt                                                                                                          <  1:78381/225‑1 (MQ=255)
                                                                                                                                             aGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTagttcggaagggcgtcgtgtagggg                                                                                 >  2:78381/1‑226 (MQ=255)
                                                                                                                                                                          aaCGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTcgc                                                    >  1:182570/1‑251 (MQ=255)
                                                                                                                                                                          aaCGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTcg                                                     <  1:57010/250‑1 (MQ=255)
                                                                                                                                                                                    caacCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGccc                                            >  1:4730/1‑249 (MQ=255)
                                                                                                                                                                                              aGTGTTGAAGTTCGGCGGTACAGCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCGGAAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGGGGCGATGATTGAAAAAACCATTAGCGACCAGGATGCTTTACCCCATATCAGCGATGCCGAACGGATTTTTGCCGAACTTTTGCCGGGACCCGCCCCCGCCCAGCCgggaatac                                >  2:192985/1‑246 (MQ=255)
                                                                                                                                                                                                   tGAAGTTCGGCGGTACATAAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCAAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCt                            <  1:268501/250‑1 (MQ=255)
                                                                                                                                                                                                    gAAGTTCGGCGTTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTg                           <  1:114903/250‑1 (MQ=255)
                                                                                                                                                                                                                            aaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCTCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATc  >  1:349926/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                     
CGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGGTTCCCGCTGGCGCAATTGAAAACTTTCGTCGATC  >  NC_000913/152‑622

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