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_0009139690TG93.3% 79.4 / ‑3.1 30A211A (GCT→GCG) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (23/5);  total (24/6)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

CCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATC  >  NC_000913/731‑1198
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
ccATTAAGGCCGGCGTATTAGACGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGTTATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGATTCACTGCCGGTAATGACAAAGGCGACCTGGTGTTACTTTGACGCAACGGTTCCGACTACTCCGCGGCGGTgccggct                                                                                                                                                                                                                            >  1:302491/1‑251 (MQ=255)
              gTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTCCTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGATGCCTGTTTACgcgc                                                                                                                                                                                                              >  1:21671/1‑251 (MQ=255)
                     aaGCGCGAGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGACATTACCTCGAATCTACCGTCGAAATTGCTGAGTCCCCCCGCCGTATTGCGCCAAGTCGCATTCCGGCTGCTCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGCCTACTACGCGGCGGAGCTGGCTGCCAGTTTATGCGACGActgt                                                                                                                                                                                                       >  1:313276/1‑251 (MQ=255)
                             ggTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGGTGGGTGCCCTTTTTCGCGCCCGTTGTTGCgaggt                                                                                                                                                                                               >  2:175749/1‑249 (MQ=255)
                                 acaACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTACACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTAATCCGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGctagcttttg                                                                                                                                                                                           >  2:142420/1‑242 (MQ=255)
                                 acaACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTgg                                                                                                                                                                                           >  1:318425/1‑251 (MQ=255)
                                      gTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTggacggg                                                                                                                                                                                      >  1:183809/1‑250 (MQ=255)
                                       ttACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCAACCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGCTGGCAGGTTTCACCGCCGGTCATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCAGCGGCGGGTCTGGCTGCCTGTTTTAGCGCCGGttgttggcggtattgtgaggac                                                                                                                                                                                     >  2:466144/1‑235 (MQ=255)
                                                     ccGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATCTTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCTCCGATTGTTGCGAGATTTGGATGG‑ACGTTGTCGGGGTCtt                                                                                                                                                                      >  1:46057/1‑250 (MQ=255)
                                                      cGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTTGTTACCCGCCGCCTTTTGCGTGATTTGGACGG‑GCGtttgcggagtgttt                                                                                                                                                                     >  2:184590/1‑239 (MQ=255)
                                                          cGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGGGGCGGCTGCCTGTTTACGCGCCGATTTTTGCGGGGTTGTGACGGGACGTTGAGGGGGtgtttacc                                                                                                                                                                  >  2:332005/1‑244 (MQ=255)
                                                              aaaCTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCACGTCGCATTCCGGCTGATCCCATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCCGTTTCCCCGCCGGTTGTTGCGGGGTTTGGACGG‑ACGTTGACGGGGGCTTTACCTGCgc                                                                                                                                                             >  1:163257/1‑250 (MQ=255)
                                                                   gctgGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCGGGGTGCCTGTTTTCGCGCCGATTGGTGCGGGGATTGGGCGG‑ACGTTGACGGGGGCTATACCTGCGGCCcgc                                                                                                                                                        >  2:126821/1‑251 (MQ=255)
                                                                                                  gTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTc                                                                                                                         >  1:98955/1‑251 (MQ=255)
                                                                                                                    aCCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGGGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGGTGACGGGGGCTATACCTGGGACCCGCGTCAGGTGCCCGATGCGGGGTTGGTGAAGTCGATGGCCTACCAGGAAGc                                                                                                       >  1:127743/1‑251 (MQ=255)
                                                                                                                       cgcgttgtggcgtaaagcctctttCCGGCGGATCACATGGTGCTGATGGCAGTTTTCCCCGTCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGAt                                                                                                    <  1:85033/229‑1 (MQ=255)
                                                                                                                          cgTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATggg                                                                                                 >  1:83974/1‑250 (MQ=255)
                                                                                                                              ttGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGGGCCCGATGCGAGGGTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGGGCtt                                                                                             >  1:304090/1‑251 (MQ=255)
                                                                                                                                                      aTCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTAAGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGGGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCCGCGCTTAAGttctt                                                                     >  1:239051/1‑251 (MQ=255)
                                                                                                                                                                            tttCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGGCTATACCTGCGACCCGCGTCAGGGGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTAcccccc                                               >  1:457735/1‑250 (MQ=255)
                                                                                                                                                                                   gccgGTAATGAAAAAGGCGAACTGGAGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATAACCCCCATCGccc                                         >  1:352514/1‑250 (MQ=255)
                                                                                                                                                                                      ggTAATGAAAAAGGCGAACTGGTGGTCCTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGtt                                     <  2:83974/251‑1 (MQ=255)
                                                                                                                                                                                      ggTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGCAGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGGTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCCAAAGTTCTTCACCCCCGCACCATTACCCCCCTCGCCCAGtt                                     >  1:25293/1‑251 (MQ=255)
                                                                                                                                                                                             aaaaaGGCGAACTGGTGGTACTTGGACGCACCGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCCGTTTACGCGCCGGTTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGGGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGGTGGGGCTTTCCTACTTCGGCGCTTAAGTTCTTCACCCCCGCACCCTTTCCCCCCTCGccccggtccacgtc                              >  2:230236/1‑240 (MQ=255)
                                                                                                                                                                                               aaaGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATccc                            <  1:391195/251‑1 (MQ=255)
                                                                                                                                                                                                   gCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCCGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGc                        <  1:142420/251‑1 (MQ=255)
                                                                                                                                                                                                           gtggtaCTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGCGATTTGGGCGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGGCCTACCAGGAAGCGGTGGAGCTTTCCTAATTCGGCGCTAAAGTTCTTCACCCCCGCACCCTTACCCCCATCGCCCAGTTCCCGATCCCTTGCCTGATTa                 >  2:12345/1‑250 (MQ=255)
                                                                                                                                                                                                                         cAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATc  <  2:189956/251‑1 (MQ=255)
                                                                                                                                                                                                                             ggTTCCGGCTACGCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCTCCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTaaaa              <  2:286452/235‑1 (MQ=255)
                                                                                                                                                                                                                             ggTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGCTTaaaa              >  1:286452/1‑235 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
CCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATC  >  NC_000913/731‑1198

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: ‑