Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,797 G→A Q487Q (CAG→CAA 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,7970GA97.1% 91.9 / ‑4.1 34Q487Q (CAG→CAAthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (13/20);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-01

GCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTG  >  NC_000913/1554‑2040
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
gcgcACCTTGCTTGGGATCTTGGCGAAAGTCTTTGCCGCGCTGGCCCTCGCCCATATCAACATTGGCGCCATTTCTCAGGGATCTTCTGAACGCGCAATCTCTTTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTGCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACaaatc                                                                                                                                                                                                                                               <  2:339509/250‑3 (MQ=255)
           gTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTAATCATCAGATGCTGTTCAATACCGATCAGGATATCGAGGTGTTTGTGATTGGCGTCGGGGGCGTTGGCGGTGCGCTGCTGGGGCAAATGAAGCGTCAACAAAGCGGGCTgaaga                                                                                                                                                                                                                                     >  2:137301/1‑249 (MQ=255)
                    cGGCGAAATTCTTTCCCGCGCTGGCCCGCGCCACTATCAACCTTGACGCCATTGCTCAGGGCTCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGAACACTGGCGTGCGCGTTCCTCATCAGATGCTGTTCAATACCGCTCAGGTTATCGCAGTGTTTGTGCTTGGCGTCGGTGGAGTTGTCGGTGCGCTGCTGGAGCAACTGAACCGTCAACCAAGCTGGCTGATGAATAAACCTATc                                                                                                                                                                                                                          >  2:171287/1‑251 (MQ=255)
                        gAAATTCCTTGCCGCGCTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACt                                                                                                                                                                                                                      <  1:71942/251‑1 (MQ=255)
                                ttgtcgcgcTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGAGCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTTGCTGAAGAATAAACATATCGACTTACgtgt                                                                                                                                                                                                               <  1:166565/242‑1 (MQ=255)
                                ttGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACgtgt                                                                                                                                                                                                               <  1:250346/250‑1 (MQ=255)
                                 tgccacgctgtcccccgcCAATATCAACATTGTCGCCATTGCTCAGGGTCCTTCTGAACGCTTAATCCCTGTCGTGGTACATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTTTTCAATACCGATCAGGTTATCGAAGTGTTGGTGATTTGCGTCGGTGGTGTTGGCTGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCTACTTACgtgta                                                                                                                                                                                                              <  2:358183/235‑2 (MQ=255)
                                                                gTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCAc                                                                                                                                                                              <  1:218887/251‑1 (MQ=255)
                                                                    ccATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAAt                                                                                                                                                                          >  1:244003/1‑251 (MQ=255)
                                                                     cATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGAGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGAGTTGGCGGTGCGCTGCTGGGGCAACTGTAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAAATAGAAGGCACTGCGCACCAATg                                                                                                                                                                         >  1:102613/1‑251 (MQ=255)
                                                                             aGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGACGAATAAACATATCGACTTACGTGTATGCGGTGTTGCCAACTCGAATGCACAGCTCACCAATGTGCATGGc                                                                                                                                                                 >  2:51367/1‑251 (MQ=255)
                                                                              gggATCGTCTGAACTCTCAAGCCCGGGCGTGGTAAATAACGAGGATGCGACCCCTGGCGTGCGCGTTACGCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGAGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGGGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGTCTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTTACAAATGTGCATGGcc                                                                                                                                                                <  2:360943/251‑1 (MQ=255)
                                                                                gctcttctGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCAGCAGATGCTGTTCAATAACGATCAGGTTATCGCTGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGGGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCGGGCTGAAGAATAAACAAATCGACTTACTGGTCTGCGGTGTTGCCAAATCGAAGGCACTGCTCACCAATGtgacgggcct                                                                                                                                                               >  2:239662/3‑241 (MQ=255)
                                                                                                     ctGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGga                                                                                                                                          >  1:163939/1‑250 (MQ=255)
                                                                                                       gTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGgaag                                                                                                                                        <  1:333226/250‑1 (MQ=255)
                                                                                                               aaaTAACGATGATGCGACCACTGGCGTGCGCGTTAATCATCAGATGCTGTTCAATACCGCTCAGGTTATCGAAGTGATTGTGATTGGCGTCGGTGGCGTAGGCGGTGCGCTGCTGGGGCAAATGAAGCGTCAAAAAAGCTGGCGGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGACAGAGCTGCTCACCCATGTGCATGGCCTAAATCTGGAAAAATGGCAGGAACAACTGtcgc                                                                                                                               >  2:328530/1‑251 (MQ=255)
                                                                                                                    aCGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGcc                                                                                                                          <  1:108430/251‑1 (MQ=255)
                                                                                                                            gCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCTCCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGcc                                                                                                                  <  1:276576/251‑1 (MQ=255)
                                                                                                                                   cTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGTAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAATGCACTGCTCACCACTGTGCATGGCCTATATCTGGCAAACTGGAAGGAAGAACTGGCGCAAGCCAACGAGCCGTTTAAt                                                                                                           >  1:14359/1‑251 (MQ=255)
                                                                                                                                     ggCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAAtct                                                                                                         <  1:350159/251‑1 (MQ=255)
                                                                                                                                       cGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCg                                                                                                       <  1:141237/251‑1 (MQ=255)
                                                                                                                                               ttACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGAAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCtt                                                                                                >  1:126605/1‑250 (MQ=255)
                                                                                                                                               ttACTCATCAGATGCTGCTCCATACCAATCAGGTTATCGAAGTGCTTATCATTGGCGTCGGTGGCGTTGACCGTGCGCTCCTCGAGCAGCTGAACCGTCAACAAAGATGGCTGAAGAATAAACATGTCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGAACTGCTCACAAACGTGCCTGGCCTAAATCTAGAAAACTGGCAGGGAAAACTGGCGCAAGCCAAGGAGACGATTAATCTCgagagcct                                                                                                >  2:252619/1‑243 (MQ=255)
                                                                                                                                                    catcaGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                           <  1:280262/250‑1 (MQ=255)
                                                                                                                                                               gtccAATCCCGATCAGGTTATCGAAGTGTTTGTGATTGGTGGCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCCGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCGCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAATAACGGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGaa                                                                                 <  1:215293/246‑1 (MQ=255)
                                                                                                                                                                  cAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAat                                                                            <  1:137301/251‑1 (MQ=255)
                                                                                                                                                                            cAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATctgc                                                                   <  2:244003/250‑1 (MQ=255)
                                                                                                                                                                                    cGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCg                                                           <  1:308242/250‑1 (MQ=255)
                                                                                                                                                                                                     gCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACt                                         <  1:239662/251‑1 (MQ=255)
                                                                                                                                                                                                       gTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAAAAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTc                                       >  1:319574/1‑251 (MQ=255)
                                                                                                                                                                                                             ggCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTccagcc                                  <  1:352052/250‑1 (MQ=255)
                                                                                                                                                                                                                   ggCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGt                           >  1:203298/1‑251 (MQ=255)
                                                                                                                                                                                                                           gcgtctgaaGCAACTTGAGCGTCAACAAAGCTGGCTGAAGAAGAAACATATCGACTTACGTGTCTGCTGTGTTTCCAACGCGAAGGCACGGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGATCCGTTTAAGCTCGGGCGCTTAATTCGCCTCTTGAAAGAATATCATCTGCTGAACCCGGTCAGTGTTGACTGCACTTCCAGCCAGGCATTGGCGGATc                   <  2:371715/243‑1 (MQ=255)
                                                                                                                                                                                                                                            aGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTcttg  <  1:203420/251‑4 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
GCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTG  >  NC_000913/1554‑2040

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