Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 867 C→T S177S (AGC→AGT thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009138670CT96.5% 77.4 / ‑4.1 29S177S (AGC→AGTthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (28/0);  total (29/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

GAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGG  >  NC_000913/625‑1102
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            
gAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGCTATTGCTGAGTCCACCCGCCGTATTGCGTCAAGTCGCATTcc                                                                                                                                                                                                                                      >  1:437118/1‑251 (MQ=255)
             tAAAACATGTCCTGCATGGCATTAGTTTGTTGGGACAGTGCCCAGATAGCCTCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTCGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCCAGTCCCATTCCGGCTGATCACATg                                                                                                                                                                                                                         >  2:71102/1‑251 (MQ=255)
                      tCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGTTCCCATGGTGCTggtg                                                                                                                                                                                                                >  2:495253/1‑248 (MQ=255)
                      tCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTCCCGTGGCGAGAAACTATCGATCGCCATTATGACCGGCGTCTTAGCAGCGCGCGGTCACACCGTTACCGTTATCGATCCGGTCGACAAAATGCTGGCAGTGGGGCATTACCTCGAATCTACCATCGCTATTGCTGAATCCACCCGCCGAATTGCGGCCAGTCGCATTACGGCTGATCACATggccctgctg                                                                                                                                                                                                                >  1:228718/1‑243 (MQ=255)
                         tGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGc                                                                                                                                                                                                              >  1:428158/1‑250 (MQ=255)
                           cATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAgg                                                                                                                                                                                                           >  1:196053/1‑251 (MQ=255)
                                cATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGCTATTGCTGAGGCCACCCGCCGTATTGCGGCCAGTCGCATTCCGGCTGATACCATGGTGCTTATGGCCGGTTTc                                                                                                                                                                                                       >  2:339316/1‑250 (MQ=255)
                                      tttGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCAccgccg                                                                                                                                                                                                >  1:396160/1‑251 (MQ=255)
                                                  gTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGCTATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGATGATGGAAGGTTTCACCGCCGGTAATGacaacg                                                                                                                                                                                    >  2:408436/1‑246 (MQ=255)
                                                       cGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGa                                                                                                                                                                                >  1:121641/1‑250 (MQ=255)
                                                               aTCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGAGTATTAGAAGCGCGCGGTCACAACGATACCGTTATAGATCCGGTCGAAAAACTGCTGGCAGTGGGGAATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCACATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTCATGAAAAAGGCGAACtgggggg                                                                                                                                                                       >  1:506650/1‑247 (MQ=255)
                                                                  aaCGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTggtggtact                                                                                                                                                                    >  1:297945/1‑248 (MQ=255)
                                                                      ctgcgctgATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCgg                                                                                                                                                                                               >  1:100203/1‑220 (MQ=255)
                                                                       tgcgctgATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGTCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGAc                                                                                                                                                               >  1:43447/1‑251 (MQ=255)
                                                                                       ggCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGCTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGCGCTGATGGCAGGGTTCACCGCCGGTAATGAAAAAGGCAGACTGGTGGGACTTGGACGCACCGGTTCCGacc                                                                                                                                                >  2:482177/1‑249 (MQ=255)
                                                                                              aaaTGACGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTAACACCGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCCAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGACAAAGGCGAACTGGTGGGACTTGGACGCAACGGTTCCGACTActccgcg                                                                                                                                        >  2:17594/1‑247 (MQ=255)
                                                                                                aTGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGTTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTActccgcggc                                                                                                                                      >  1:195314/1‑245 (MQ=255)
                                                                                                    cGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGCTCACATGGTGCTGATGGCAGGTTTCCCCGCCGGGAATGAAAAAGGCGAACTGGTGAGACTTGGTCGCAACGGTTCCGACTACTCCGcggcgggg                                                                                                                                  >  2:273019/1‑249 (MQ=255)
                                                                                                      aTCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAAATGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGGgct                                                                                                                                >  2:74656/1‑251 (MQ=255)
                                                                                                                               aaGCGCGCAGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCTGGTTTCACCGCCGGTAATGAAAAAGGCGCACTGGGGGTCATTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGGCGGCTGCCTGTTTTAGCGCCGGttgt                                                                                                       >  2:283972/1‑251 (MQ=255)
                                                                                                                                aGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCgg                                                                                                                                                                                                                                    >  1:211983/1‑125 (MQ=255)
                                                                                                                                                gTTACCGTTATCGATCCGGTCGAAAAACAGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCAGACTACTCCGCGGCGGTGCTGGCTGACTGTTTTCGGGACGATTGTTGCGAGGTTTggacggg                                                                                      >  1:107071/1‑250 (MQ=255)
                                                                                                                                                  taccgTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGCTGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGCGCCGGCTGCCCGGTTACCCGCCGATTGGTGCGGGATTTgggcgggcg                                                                                    >  2:87270/5‑248 (MQ=255)
                                                                                                                                                    ccgTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGGGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCATGTTTACGCGCCGATTGTTGCGAGGTTTGGGCGGACGtt                                                                                  >  1:419674/3‑251 (MQ=255)
                                                                                                                                                                                    gTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTTCGCGCCGGTTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTccggg                                                 >  1:438957/1‑247 (MQ=255)
                                                                                                                                                                                                         aCCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGGGCCCGATGCGGGGTTGTTGa                             >  1:255160/1‑250 (MQ=255)
                                                                                                                                                                                                          ccGTCGGTATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTAATTGGACGCAACGGTTCCGACTACTCAGCGGCGGTGGTGGATGCCTGTTTACGCGCCGTTTGTTGCGGGAGTTGGGCGGGCGGTGACTGGGGTTATACCTGCGACCCGCGTCacgggaccgctgggagggtgtggggg                           >  2:38423/1‑226 (MQ=255)
                                                                                                                                                                                                                            ccaccaGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGGGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGGTTACGCGCCGATTGTTTCGAGCTTTGGACGGACGTTGCCGGGGTCTATACCAGGGACCCGCGGCAGGGGCCCCATTCGAGGGTGTTGAGGGTGATGTCCTCCCAggc          >  2:413144/1‑249 (MQ=255)
                                                                                                                                                                                                                                     gTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGGCGGACGTTGACGGGGGCTATACCTGCGACCCGCGGCAGGGGCCCGATGCGAGGTTGTTGAAGTTGATGTCCTACCAGGAAGcattgg  >  1:363392/1‑244 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            
GAATTTGCCCAAATAAAACATGTCCTGCATGGCATTAGTTTGTTGGGGCAGTGCCCGGATAGCATCAACGCTGCGCTGATTTGCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGG  >  NC_000913/625‑1102

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