Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,038 C→T L568L (CTG→TTG)  thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,0380CT91.2% 90.9 / NA 34L568L (CTG→TTG) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/16);  minor base G (0/3);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATG  >  NC_000913/1801‑2218
                                                                                                                                                                                                                                              |                                                                                                                                                                                    
aGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCCGCCAGGCAGTGGCGGATCCATATGCCGACTTCTTGCGCGAAGGttt                                                                                                                                                                         >  2:232900/1‑251 (MQ=255)
    ggctgAAGAATAAACATATCGACTTACGTGTCTGCGGTGGTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATGCGCCTCGTGAAAGAATATCATCTGCTGATCCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCAc                                                                                                                                                                     <  1:23856/251‑1 (MQ=255)
      ctgAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGCTCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGt                                                                                                                                                                   >  1:217794/1‑251 (MQ=255)
       tgAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAAATGGCGCAAGCCAAAGAGCAGTTTAATCTCGGACGCTTAATTCGACTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGGTTCCACggt                                                                                                                                                                  >  1:159352/1‑249 (MQ=255)
            aaTAAACATATCTACTTACGTGTCTGCGGTGTTGCCAACTCGAATGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCTTTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCAc                                                                                                                                                             <  1:33506/251‑1 (MQ=255)
                 aCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGACGTTTAATCTCGGGCGCTTAATTCGCCTCGTGACAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCCATATGCAGACTTCTTGAGCACAGGTTTCCACGTTGTCACGCCg                                                                                                                                                         >  2:148335/1‑250 (MQ=255)
                             aCGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGcc                                                                                                                                             <  1:44213/250‑1 (MQ=255)
                              cGTGTCGGCGGTGTTGGCAACGCGAAGGCACTGCTCACCAATGAGCATGGCCTAAATCTAGAAAACTGGCAGTAAGATCTTGCTCAAGCCAAAGATCCGTTTAATCTCGGGTGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCAGGTCATTGTTGGCTGCACTTCCAGCCAGGCATGGGCGGAGCAATATGCCGATTGCGCGCGCGACGGTTACCCCGTGGTCACTCCGAAcgaaccggccaa                                                                                                                                           <  2:110444/251‑12 (MQ=255)
                                tgtCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCCTGGCCTAAATCTGGAAAACTGGCAGGAAGAACTAGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAACGAAGATCATCTGCTGAACCCGGTCATTGTTGAGTGCACTTCCAGCCAGGCAGTGGCGGCCCAATATGCCGGCGTCTTGCGCGCAGGTTTCCACGTTGTCACGCCGCACCAAAAGGCCaaa                                                                                                                                          >  2:184166/1‑249 (MQ=255)
                                                           actGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATctgt                                                                                                              <  1:55153/250‑4 (MQ=255)
                                                             tGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTg                                                                                                            >  1:171968/1‑251 (MQ=255)
                                                                tCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCAAACACCTCGTCGATGGATTACTACAATCTGTTGCGt                                                                                                         >  1:186638/1‑251 (MQ=255)
                                                                   ccAATGTGCATGGCCTAAATCTGGAAAACTTGCAGGAAGAACTGGCGCAAGTCAAAGAGCCGGTTAAGCTCGGGCGCATAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCTCGTCCAGCCAGGCAGTGGCGGATCAATATGCCGAGTTCGTGCGCGAAGGGTTCCGCGTTGTCACGCCTAACAAAAATGCCAACACCGCGGCGAGGGAttgcttccagctgtggtgtcat                                                                                                      <  2:139935/251‑21 (MQ=255)
                                                                               gCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGGTGTCACGCCGAACAAAAAGGCCAACAACTCGTCGATGGATTACTACCATCTGTTTCGTCATGCGGCGGaaaa                                                                                           >  1:161953/1‑250 (MQ=255)
                                                                                    aaTCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATcgcg                                                                                     <  1:1273/251‑1 (MQ=255)
                                                                                         ggAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCCTCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCTCGTTGTCACGCCGAACAAAAAGGCCAACACCGCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCgcgg                                                                                    <  1:8281/247‑1 (MQ=255)
                                                                                                 ggCAGGAAGAACTTGCTCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTTTTGACTTCACTTCCAGCCAGGCAGTGGCGGATCAATATGCGGACTTCTTGCGCGAAGGTTGCCACGTGGTCTCGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAAtt                                                                            <  1:43917/247‑1 (MQ=255)
                                                                                                                cgcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGGCTTCTTGCGCGCAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGtt                                                         >  1:124544/1‑251 (MQ=255)
                                                                                                                            agCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGAtt                                              >  1:182701/1‑250 (MQ=255)
                                                                                                                              ccGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTAcc                                           <  1:18652/251‑1 (MQ=255)
                                                                                                                              ccGTTTAATCGCGGTGGCTTTATTCGCCTCGTTAAAGAATATCAGCTGCTGAACCGCGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCAGTTCAATAGGCCGACTTCGTGCGCGAAGGTGGCCATGTTGTCACGCCGAACAAAAATGTCAACGCCTCGGCGATGGTTTACTTTTATCTGTTTCGTCATGCGGCGGATAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGAGGCTGGATTAcc                                           <  1:239620/251‑1 (MQ=255)
                                                                                                                                         ccggCGCTTAATTCGCCTCGTGAAAGAATACCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGGCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGGCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAAGCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGTCTGGATTACCGTTTATTgag                                 <  2:126537/248‑1 (MQ=255)
                                                                                                                                         cGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTAtt                                    <  1:97620/247‑1 (MQ=255)
                                                                                                                                         cGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGGTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTgaga                                <  2:171968/251‑1 (MQ=255)
                                                                                                                                               cTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGa                                                                                                                            >  1:213701/1‑153 (MQ=255)
                                                                                                                                               cTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGa                                                                                                                            <  2:213701/153‑1 (MQ=255)
                                                                                                                                                         cTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCCATATGCCGTCTTCTTGCGCGAAGGTTTCCACGTTGTCCCGCCGAACAAAAAGGCCAACACCTCGTCGATGGCTTACTACCATCTGtt                                                                                                             >  2:123112/1‑158 (MQ=255)
                                                                                                                                                         cTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGtt                                                                                                             <  1:123112/158‑1 (MQ=255)
                                                                                                                                                            gTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTggg                                                      >  1:201534/1‑210 (MQ=255)
                                                                                                                                                            gTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTggg                                                      <  2:201534/210‑1 (MQ=255)
                                                                                                                                                             tGAAAGAATATCATCTGCTGACCCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATCGCGGATTCTTGCGCGAAGGTTTCCACGGTGTCACGCAGAACAAACAGGCCAACACCTCGTCGAAGGATTACTACCATATGTTGCGTCCTGCGGCGGAAAAATCGAGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAACAACTGCTCacg            >  2:264836/1‑249 (MQ=255)
                                                                                                                                                                        catcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGtgatg  <  1:176866/250‑1 (MQ=255)
                                                                                                                                                                                          aTTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCCCGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGtt                                                                                                             >  2:234208/1‑125 (MQ=255)
                                                                                                                                                                                          aTTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGtt                                                                                                             <  1:234208/125‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                    
AGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATG  >  NC_000913/1801‑2218

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: