Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,323 C→T V329V (GTC→GTT 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,3230CT93.8% 82.9 / ‑3.0 32V329V (GTC→GTTthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (15/15);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

GATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGC  >  NC_000913/1080‑1555
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                        
gATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTACCAGATCCCTTGCCTGCTTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGACGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGACTAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATgg                                                                                                                                                                                                                                   >  2:147601/1‑251 (MQ=255)
   gtccttacAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCCCCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTGCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATggcg                                                                                                                                                                                                                                 <  1:144083/243‑1 (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)
                                                                                                                           agctccAGGTACGCGCATTGGTGCCAGCCTTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTTCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACggg                                                                                                        <  1:71846/247‑1 (MQ=255)
                                                                                                                            gctccAGGTACGCTCATTGGTGCCAGCCGTGATGACGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGACAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCCGCGCTGTCACGCGCCCGTATTTCCGTGGTGATGATTACGCAATCATCTTCCGGATCCAGTATCAGTTTCTGCGTTCCGCAAAGCGCCTGTGTGCGAGCTGAACGGGc                                                                                                       >  1:59534/4‑251 (MQ=255)
                                                                                                                               ccatggaCGCTCCTTGGTGCCCGGCGTGATGAAGCCGAAGTACCGGTCATGGGCATTGCCAATCTAAATAAGATGGCAATGTACAGCGTTTCCGGCCCGGGGATGAAAGGGATGGGTGGCATGGCGGCGCGCGTGTTTGCAGCTATGGCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCa                                                                                                      <  2:131743/243‑1 (MQ=255)
                                                                                                                               ccAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAAt                                                                                                    <  1:321434/251‑1 (MQ=255)
                                                                                                                                           aTTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGtt                                                                                        >  1:840/1‑251 (MQ=255)
                                                                                                                                                       cGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGTCGCGCGTGTTTTCAGCGATGTCCCGCGCCCGTATTTCCGTGGTGCTGCTTACGCAATCCTCTTCCGAATACAGAATCAGTTTCTGCGTTCCTCAAAGCGACTGTGTGCGCTCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACt                                                                            >  1:307711/1‑251 (MQ=255)
                                                                                                                                                          gatgaAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCTGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGaa                                                                         <  1:15845/251‑1 (MQ=255)
                                                                                                                                                                    aaTTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGCAGGGGATGGTCGGCATGGCGGAGCGCGTCTTTGCAGCGATGTCACGCGCCAGTATTTCCGGGGTGCTGATTACGCAATCATCTTCCGACTACAGCATAAGTTTCTGCTGTCCCCCAAGCGACTGAGGGCGAGCTGAACGGGCAATTCAGGAACAGTTCTACCTGGAAATGAaagaagacgt                                                                >  2:250346/1‑246 (MQ=255)
                                                                                                                                                                      ttACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCAGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGATGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGACGGCTTACt                                                             >  1:161515/1‑251 (MQ=255)
                                                                                                                                                                             tCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCg                                                      >  1:292225/1‑251 (MQ=255)
                                                                                                                                                                                    cATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGcgg                                               <  1:177631/251‑3 (MQ=255)
                                                                                                                                                                                         ccAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGCTGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGAGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCCGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGGGCCGCTGGCGGTGACg                                          >  1:357986/1‑251 (MQ=255)
                                                                                                                                                                                         ccAATCTGAATAATATGGCAATGGTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGAc                                           <  1:147601/250‑1 (MQ=255)
                                                                                                                                                                                                gAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAatcatc                                                                                                                                                                   <  1:177701/123‑1 (MQ=255)
                                                                                                                                                                                                    aataTGGCAATGGTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTAGCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTgg                                <  2:840/247‑1 (MQ=255)
                                                                                                                                                                                                                       tttccggcccGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGACGGAACGGCGGGCCCTTATCTCggggggcggg            >  1:360943/9‑243 (MQ=255)
                                                                                                                                                                                                                                 ggggATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCCTTTTCTCGGTGGGAGGGGATGGTATgg  >  1:11471/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                        
GATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGC  >  NC_000913/1080‑1555

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: