Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,875 A→G V513V (GTA→GTG)  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,8750AG98.3% 204.3 / ‑4.1 58V513V (GTA→GTG) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (26/31);  minor base T (1/0);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

TGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTT  >  NC_000913/1739‑2002
                                                                                                                                        |                                                                                                                               
tgtTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCAt                                                                                                                              >  2:473081/1‑140 (MQ=255)
    tgtgATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGcc                                                                                                                          >  2:687574/1‑140 (MQ=255)
     gtgATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCt                                                                                                                         <  2:1093016/140‑1 (MQ=255)
     gtgATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCt                                                                                                                         >  2:118223/1‑140 (MQ=255)
        aTTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTaa                                                                                                                      >  1:1043401/1‑140 (MQ=255)
          tGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATc                                                                                                                    >  1:304359/1‑140 (MQ=255)
            gCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTg                                                                                                                  >  1:999210/1‑140 (MQ=255)
            gCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTg                                                                                                                  >  1:1189762/1‑140 (MQ=255)
            gCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTg                                                                                                                  >  1:186592/1‑140 (MQ=255)
            gCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTg                                                                                                                  >  2:259356/1‑140 (MQ=255)
                  gTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAAc                                                                                                            >  1:672695/1‑140 (MQ=255)
                   tGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACt                                                                                                           <  2:428742/140‑1 (MQ=255)
                    ggCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTg                                                                                                          <  2:142600/140‑1 (MQ=255)
                     gCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTgg                                                                                                         <  1:82678/140‑1 (MQ=255)
                          ggCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGga                                                                                                    <  1:1084201/140‑1 (MQ=255)
                                   ctgctgGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGgc                                                                                           <  2:738543/140‑1 (MQ=255)
                                      ctgGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCa                                                                                        <  2:985339/140‑1 (MQ=255)
                                      ctgGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCa                                                                                        <  1:118223/140‑1 (MQ=255)
                                             aaCTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCAGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCaaa                                                                                 >  1:310055/1‑140 (MQ=255)
                                              aCTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAag                                                                                >  1:921070/1‑140 (MQ=255)
                                                tGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAagag                                                                              <  2:1024774/140‑1 (MQ=255)
                                                tGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAagag                                                                              <  2:742018/140‑1 (MQ=255)
                                                tGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAagag                                                                              >  1:327322/1‑140 (MQ=255)
                                                tGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAagag                                                                              <  2:1090672/140‑1 (MQ=255)
                                                tGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAagag                                                                              <  2:912140/140‑1 (MQ=255)
                                                tGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAagag                                                                              <  2:327322/140‑1 (MQ=255)
                                                                ctggctgAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGgc                                                              >  2:489739/1‑140 (MQ=255)
                                                                           aTAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGc                                                   <  1:987334/140‑1 (MQ=255)
                                                                             aaaCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCt                                                 <  2:666493/140‑1 (MQ=255)
                                                                               aCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCg                                               <  2:304359/140‑1 (MQ=255)
                                                                                  tatCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGa                                            <  1:473081/140‑1 (MQ=255)
                                                                                  tatCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGa                                            <  2:1189762/140‑1 (MQ=255)
                                                                                           aCGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGa                                        >  1:95391/1‑135 (MQ=255)
                                                                                           aCGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGa                                        <  2:95391/135‑1 (MQ=255)
                                                                                           aCGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATatc                                   >  1:617445/1‑140 (MQ=255)
                                                                                            cGTGTCAGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATatca                                  <  2:310055/140‑1 (MQ=255)
                                                                                              tgtCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTTCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATatcatc                                >  2:126706/1‑140 (MQ=255)
                                                                                              tgtCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATatcatc                                >  2:28940/1‑140 (MQ=255)
                                                                                                  tGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGaa                                           <  2:418784/125‑1 (MQ=255)
                                                                                                  tGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGaa                                           <  2:577635/125‑1 (MQ=255)
                                                                                                  tGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGaa                                           >  1:577635/1‑125 (MQ=255)
                                                                                                  tGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGaa                                           >  1:418784/1‑125 (MQ=255)
                                                                                                  tGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATctgct                            >  1:1171083/1‑140 (MQ=255)
                                                                                                  tGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATctgct                            >  1:88333/1‑140 (MQ=255)
                                                                                                     ggTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGaa                         <  2:102969/140‑1 (MQ=255)
                                                                                                      gtgtTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAAc                        >  2:1038704/1‑140 (MQ=255)
                                                                                                      gtgtTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAAc                        <  1:489739/140‑1 (MQ=255)
                                                                                                         ttGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCg                     >  2:893146/1‑140 (MQ=255)
                                                                                                          tGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCgg                    <  2:186592/140‑1 (MQ=255)
                                                                                                           gCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGt                   >  1:1140803/1‑140 (MQ=255)
                                                                                                             cAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCa                 <  1:259356/140‑1 (MQ=255)
                                                                                                                 tCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCAttgt             >  1:212843/1‑140 (MQ=255)
                                                                                                                      ggCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACt        >  1:270194/1‑140 (MQ=255)
                                                                                                                       gCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTg       <  2:709906/140‑1 (MQ=255)
                                                                                                                       gCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTg       <  2:1140803/140‑1 (MQ=255)
                                                                                                                       gCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTg       <  2:921070/140‑1 (MQ=255)
                                                                                                                       gCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTg       <  1:687574/140‑1 (MQ=255)
                                                                                                                           tGCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACt   >  1:204420/1‑140 (MQ=255)
                                                                                                                            gCTCACCAATGTGCATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACtt  <  2:1043401/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                               
TGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTT  >  NC_000913/1739‑2002

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑