Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 969 T→G A211A (GCT→GCG)  thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009139690TG96.0% 65.4 / ‑4.1 25A211A (GCT→GCG) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (16/8);  total (17/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

CGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCT‑AAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACG  >  NC_000913/725‑1215
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
cGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACACCGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGAAATGAAAAAGGCGAACTGGTGGTACTTGGGCGCAACGGTTCCGACTCCTCCGCGGCGGTg                                                                                                                                                                                                                                                   >  2:150776/1‑251 (MQ=255)
     gCCATTATGGCCGGCGTATTAGAAGCGCGCGGGCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTACGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGCCTACTCCGCGGCGGAgctggg                                                                                                                                                                                                                                              >  1:318659/1‑250 (MQ=255)
      ccATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAACAACTGCTGGCAGAGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGATGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTgctggct                                                                                                                                                                                                                                             >  1:44685/1‑251 (MQ=255)
              gCCGGCGTATTAGAAGAGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTggtt                                                                                                                                                                                                                                     >  1:329575/1‑248 (MQ=255)
                      aTTAGAAGCGAGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAACATGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGTGTCCACCCGCCGTATTGAGGAAAGTCGCATTCCGGCTGATCCCATGGTGCTGATGGCAGGTTTAACCGCCGGTAATGAAAAAGGCGTACTGGTGGTCCTGGGACGCAACGGTTCAGACTACGCCGCGGCGGTGAGGGCTACCTGTTTACGCGCCg                                                                                                                                                                                                                             >  2:226888/1‑251 (MQ=255)
                                           cGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCAGACTACTCCGCGGCGGTGCTGGCTGCCTGATTGAGCGGCGATTTTTGCGAGATTTggacgg                                                                                                                                                                                                        >  2:272852/1‑251 (MQ=255)
                                                                           tgGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGATTCCACCAGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGATTCACCGCCGGTAATGAAACAGGCGAACTGGTGGTACTTGGAAGAAACGGTTCCGACTAATCCGCGGCGGTGCTCGCTGGCTGGTTACTCGACGATTGTTACGGGATGTGGACGTACGTTGAGGGGGGCTTAACCTGCGAccccgtg                                                                                                                                                                        >  2:90780/1‑247 (MQ=255)
                                                                                       aTTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGAt                                                                                                                                                            >  1:105616/1‑251 (MQ=255)
                                                                                                   cTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTTTTGCGGGATTTGGGCGGACGTTGACGGGGTCTATACCTGCGACCCACGGCCGGTGCCCGATGCGAGGttgttg                                                                                                                                                >  1:338988/1‑251 (MQ=255)
                                                                                                           gATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCCCCGATTGTTGCGCGATTTTGCCGCACGGTGCCGGGGTCTATACATGAGCCCCCCGTCAGGGGCCCGATTCGAAGTTGATGGAGTCggt                                                                                                                                        >  1:247810/1‑249 (MQ=255)
                                                                                                                                   aTTGCGGCAAGTCGCATCCCGGCTGATCACAGGGTGCTGATGGCAGTTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGc                                                                                                                 <  1:93164/250‑1 (MQ=255)
                                                                                                                                         gCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTACTGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGCTGGAGCTTTCCt                                                                                                           >  1:129529/1‑250 (MQ=255)
                                                                                                                                             gtcGCATTCCGGTTGATCACATGGTGCGGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTc                                                                                                      <  1:208118/249‑1 (MQ=255)
                                                                                                                                                          tGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACCCGCCGATTGTTGCGAGGTTTGGAGGGGCGTTGGCGGGGTCTATACCTTCGACCCGCGGCAGGGGCCCGATGCGGGGGTTTTGAAGCCGATTCCCTACCAGGAAGGGATGGAGCTTTCCCATTTCAGCGTT‑AAAGttt                                                                                        >  1:178431/1‑250 (MQ=255)
                                                                                                                                                              cacaTGGTGCTGAGGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCT‑AAAGTTCTTCa                                                                                    <  1:54952/251‑1 (MQ=255)
                                                                                                                                                                        tGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGGTTGTTGCGGGATTTGGACGGACGTTGACGGGGGCTTTTCCTGCGACCCGCGTCAGGTGCCCGATGCggggg                                                                                                                                                     >  1:197309/1‑173 (MQ=255)
                                                                                                                                                                              cAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCT‑AAAGTTCTTCACCCCCGCACCATTaac                                                                    >  1:10682/1‑249 (MQ=255)
                                                                                                                                                                                    tCACCGCCGGTAATGAAAAAGGCGAACTGGTGTTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCT‑AAAGTTCTTCACCCCCGCACCATTACCCCCAt                                                               <  1:204621/250‑1 (MQ=255)
                                                                                                                                                                                       ccgccgGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGCCTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCT‑AAAGTTCTTCACCCCCGCACCATTACCCCCCTCGcc                                                           >  1:373718/1‑251 (MQ=255)
                                                                                                                                                                                        cgcaggTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTTCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCT‑AAAGTTCTTCACCCCCGCACCATTACCCCCATCGccc                                                          <  2:96321/247‑1 (MQ=255)
                                                                                                                                                                                         gccgGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCT‑AAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCa                                                         <  1:166970/251‑1 (MQ=255)
                                                                                                                                                                                                                ggtggtACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCT‑AAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGAt                                    <  1:44618/249‑1 (MQ=255)
                                                                                                                                                                                                                           gACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCT‑AAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGa                       <  1:245119/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ggTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTTACGGGGTCTATACCTGCGAACCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGACCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTTAAAGT‑CTTCACCCCCGCACCATTACCCCCCTAGCCCACTTCCAGAACCCCTGCCTGATTTAAAATACCGGaaaacattc               >  1:374959/1‑245 (MQ=255)
                                                                                                                                                                                                                                            tactCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCT‑AAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAGTACCGGAAATCCTCAAGCACCAgg      >  1:110783/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ccgctgcGGTGGTGGCTGCCGGTTTACGCGCCGATTGTTGCGAGATGTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGGCAGGTTCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCT‑AAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACg  <  2:44685/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
CGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCT‑AAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACG  >  NC_000913/725‑1215

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