Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,206 A→T A290A (GCA→GCT thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,2060AT85.3% 81.2 / 2.7 34A290A (GCA→GCTthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base T (17/12);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACG  >  NC_000913/961‑1451
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
tactCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAGTACCGGAAATCCTCAAGCACCAgg                                                                                                                                                                                                                                                  >  1:110783/1‑251 (MQ=255)
    ccgctgcGGTGGTGGCTGCCGGTTTACGCGCCGATTGTTGCGAGATGTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGGCAGGTTCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACg                                                                                                                                                                                                                                              <  2:44685/249‑1 (MQ=255)
    ccgcggcGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTAc                                                                                                                                                                                                                                               <  1:150776/245‑1 (MQ=255)
              gctggctgCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACTTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAACGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCATTGCCTGATTAAAACTACCGGAAATCATCAAGCTCCAGGTACGCTCATTGGTg                                                                                                                                                                                                                                    >  1:200139/1‑251 (MQ=255)
               ctggctgCATGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACATGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGCAGACCTACCAGGAAGCGGTGGAGCTTTCCTACTTCGGCGCAAAAGTTATTCACCCCCGCACCATTACCCCCCTCGCCCAGTTCCAGCTCCATTGCCTGATTAAAAATACCGGAAATCCTCAAGCACAAGTTACGCTCATTGGTGc                                                                                                                                                                                                                                   >  2:67663/1‑251 (MQ=255)
                                   cGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAAt                                                                                                                                                                                                               >  1:131743/1‑251 (MQ=255)
                                    ggttttttcGAGATTGTGACGGACGTTGACGGGGTCTATGCCTGCGATCTGCGTCAGGTGCCCGATGCGAGGTTGTTGAATGCGATTGCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGTAAAGCCTCAAGCATCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGaagt                                                                                                                                                                                                              <  2:77444/243‑3 (MQ=255)
                                                     gacggacgTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCAttt                                                                                                                                                                                             >  1:152520/1‑251 (MQ=255)
                                                     gacggacgTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCAttt                                                                                                                                                                                             >  1:311135/1‑251 (MQ=255)
                                                               gACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGa                                                                                                                                                                                    <  1:118926/250‑1 (MQ=255)
                                                                                   cccGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCg                                                                                                                                                               >  1:358183/1‑251 (MQ=255)
                                                                                                cccGATGCGAGGTTGGTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCggg                                                                                                                                                  <  1:247745/251‑1 (MQ=255)
                                                                                                 ccGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCgggg                                                                                                                                                 <  1:373678/251‑1 (MQ=255)
                                                                                                     tGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATg                                                                                                                                              >  1:292430/1‑250 (MQ=255)
                                                                                                              tttgAAGTCGATGTCCTACCAGGCAGCGATGGAGCTTTCATACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTGCCAGATCCCTTGCCTGATTAAAAATACGGGAAATCCTAAAGCTCCAGGTACGCTCATTGGTGCCAGCCGAGATGAAGACGAATTACCGGTTAAGGAAATTTCAACTCTGACTACTATGGCGATGTTCAGCGTTCCCGGCGCGGGGATGAAGGGGTTgg                                                                                                                                    >  2:141237/2‑251 (MQ=255)
                                                                                                                   aGTCGCTGTCCTACCAGGAAGCGATGGCGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGCTggctggc                                                                                                                               >  1:294027/1‑246 (MQ=255)
                                                                                                                      cGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATaat                                                                                                                                                                               >  1:145069/1‑199 (MQ=255)
                                                                                                                      cGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATaat                                                                                                                                                                               <  2:145069/200‑2 (MQ=255)
                                                                                                                       gATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTACCAGATCCCTTGCCTGCTTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGACGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGACTAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATgg                                                                                                                           >  2:147601/1‑251 (MQ=255)
                                                                                                                          gtccttacAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCCCCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTGCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATggcg                                                                                                                         <  1:144083/243‑1 (MQ=255)
                                                                                                                               tccAGGAAGCGATGGAGCTTTCCTGCTTCGGCGCTCAAGTTCTTAACCACCGCACCATTACCCCCCTCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACAGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGGGATGAAGACGAATTAACGGTCAAGAGCCTTTCCTATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATgtgtggcatggaggtagcg                                                                                                                   >  2:71942/2‑233 (MQ=255)
                                                                                                                                 cAGGAAGCGATTGAGCCTTCCTTCCTCGTCGCCAAAGTTCTTCACCCCCGCACCATTCCCCCCATCGCCCATGTCCAGATCCCTGGCCTGATTAAAAGTACCGGAAAACCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATATCATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGTGATGGTCGGCATGGCGGTGCGCGt                                                                                                                 <  2:110783/251‑1 (MQ=255)
                                                                                                                                   ggAAGCGATGGACCTTTCCTACTTCGGCGCTAAATTTCTTCCCCCCCGCATCATTTCCCCCATCGCCCAGTTCCAGATCCCGTGCCATATTAAAAATACCGGAAATCCTCAATCACCAGGTACGCTCATGGGGGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTGCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGtgt                                                                                                               <  2:292430/251‑3 (MQ=255)
                                                                                                                                       gCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGc                                                                                                           <  1:90780/251‑1 (MQ=255)
                                                                                                                                                  ttCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCAc                                                                                                >  1:326697/1‑251 (MQ=255)
                                                                                                                                                                            ccccccGCACCATTACCCCCACCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGAtt                                                                      <  1:106871/250‑1 (MQ=255)
                                                                                                                                                                                           cccccATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTcc                                                       <  1:226888/251‑1 (MQ=255)
                                                                                                                                                                                                       aGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATc                                           >  1:288428/1‑251 (MQ=255)
                                                                                                                                                                                                             aGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCCTCTTCCGAATTCAGTATCAGTTTc                                     >  1:75027/1‑251 (MQ=255)
                                                                                                                                                                                                              gATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCt                                    <  1:272852/251‑1 (MQ=255)
                                                                                                                                                                                                                              agaaagtccGGAATTCCTCAACCTCCAGTTACGCTCATTTGTTCCAGCCGTGATGATGACGATTTACCGGTCAAGGGCAGTTTCAATCTGAATAATATGGCAATGGTCAGCGTTTCCGGCCCGGGGATGAAAGGTATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCg                     <  2:105616/243‑1 (MQ=255)
                                                                                                                                                                                                                                         aaaTCCTCAAGCTCCAGGTACGCTCATTGGAGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAGTGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGAGTGTTTGCAGCGATTTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCCTCTTCCGAATACAGTTTCAGTTTCTGCGTTTCGCAAAGCGACTGTGTGctgg         >  1:316994/1‑248 (MQ=255)
                                                                                                                                                                                                                                             ccTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTgg     >  1:286869/1‑250 (MQ=255)
                                                                                                                                                                                                                                                cAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACg  >  1:50123/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
TACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACG  >  NC_000913/961‑1451

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