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

AGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATG  >  NC_000913/1801‑2218
                                                                                                                                                                                                                                              |                                                                                                                                                                                    
aGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCCGCCAGGCAGTGGCGGATCCATAT‑GCCGACTTCTTGCGCGAAGGttt                                                                                                                                                                         >  2:232900/1‑251 (MQ=255)
    ggctgAAGAATAAACATATCGACTTACGTGTCTGCGGTGGTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATGCGCCTCGTGAAAGAATATCATCTGCTGATCCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCAc                                                                                                                                                                     <  1:23856/251‑1 (MQ=255)
      ctgAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGCTCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGt                                                                                                                                                                   >  1:217794/1‑251 (MQ=255)
       tgAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAAATGGCGCAAGCCAAAGAGCAGTTTAATCTCGGACGCTTAATTCGACTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGGTTCCACggt                                                                                                                                                                  >  1:159352/1‑249 (MQ=255)
            aaTAAACATATCTACTTACGTGTCTGCGGTGTTGCCAACTCGAATGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCTTTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCAc                                                                                                                                                             <  1:33506/251‑1 (MQ=255)
                 aCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGACGTTTAATCTCGGGCGCTTAATTCGCCTCGTGACAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCCATAT‑GCAGACTTCTTGAGCACAGGTTTCCACGTTGTCACGCCg                                                                                                                                                         >  2:148335/1‑250 (MQ=255)
                             aCGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGcc                                                                                                                                             <  1:44213/250‑1 (MQ=255)
                              cGTGTCGGCGGTGTTGGCAACGCGAAGGCACTGCTCACCAATGAGCATGGCCTAAATCTAGAAAACTGGCAGTAAGATCTTGCTCAAGCCAAAGATCCGTTTAATCTCGGGTGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCAGGTCATTGTTGGCTGCACTTCCAGCCAGGCATGGGCGGAGCAATAT‑GCCGATTGCGCGCGCGACGGTTACCCCGTGGTCACTCCGAAcgaaccggccaa                                                                                                                                           <  2:110444/251‑12 (MQ=255)
                                tgtCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCCTGGCCTAAATCTGGAAAACTGGCAGGAAGAACTAGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAACGAAGATCATCTGCTGAACCCGGTCATTGTTGAGTGCACTTCCAGCCAGGCAGTGGCGGCCCAATAT‑GCCGGCGTCTTGCGCGCAGGTTTCCACGTTGTCACGCCGCACCAAAAGGCCaaa                                                                                                                                          >  2:184166/1‑249 (MQ=255)
                                                           actGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATctgt                                                                                                              <  1:55153/250‑4 (MQ=255)
                                                             tGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTg                                                                                                            >  1:171968/1‑251 (MQ=255)
                                                                tCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCAAACACCTCGTCGATGGATTACTACAATCTGTTGCGt                                                                                                         >  1:186638/1‑251 (MQ=255)
                                                                   ccAATGTGCATGGCCTAAATCTGGAAAACTTGCAGGAAGAACTGGCGCAAGTCAAAGAGCCGGTTAAGCTCGGGCGCATAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCTCGTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGAGTTCGTGCGCGAAGGGTTCCGCGTTGTCACGCCTAACAAAAATGCCAACACCGCGGCGAGGGAttgcttccagctgtggtgtcat                                                                                                      <  2:139935/251‑21 (MQ=255)
                                                                               gCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGGTGTCACGCCGAACAAAAAGGCCAACAACTCGTCGATGGATTACTACCATCTGTTTCGTCATGCGGCGGaaaa                                                                                           >  1:161953/1‑250 (MQ=255)
                                                                                    aaTCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATcgcg                                                                                     <  1:1273/251‑1 (MQ=255)
                                                                                         ggAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCCTCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCTCGTTGTCACGCCGAACAAAAAGGCCAACACCGCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCgcgg                                                                                    <  1:8281/247‑1 (MQ=255)
                                                                                                 ggCAGGAAGAACTTGCTCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTTTTGACTTCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCGGACTTCTTGCGCGAAGGTTGCCACGTGGTCTCGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAAtt                                                                            <  1:43917/247‑1 (MQ=255)
                                                                                                                cgcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGGCTTCTTGCGCGCAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGtt                                                         >  1:124544/1‑251 (MQ=255)
                                                                                                                            agCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGAtt                                              >  1:182701/1‑250 (MQ=255)
                                                                                                                              ccGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTAcc                                           <  1:18652/251‑1 (MQ=255)
                                                                                                                              ccGTTTAATCGCGGTGGCTTTATTCGCCTCGTTAAAGAATATCAGCTGCTGAACCGCGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCAGTTCAATAG‑GCCGACTTCGTGCGCGAAGGTGGCCATGTTGTCACGCCGAACAAAAATGTCAACGCCTCGGCGATGGTTTACTTTTATCTGTTTCGTCATGCGGCGGATAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGAGGCTGGATTAcc                                           <  1:239620/251‑1 (MQ=255)
                                                                                                                                         ccggCGCTTAATTCGCCTCGTGAAAGAATACCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGGCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGGCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAAGCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGTCTGGATTACCGTTTATTgag                                 <  2:126537/248‑1 (MQ=255)
                                                                                                                                         cGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTAtt                                    <  1:97620/247‑1 (MQ=255)
                                                                                                                                         cGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGGTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTgaga                                <  2:171968/251‑1 (MQ=255)
                                                                                                                                               cTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGa                                                                                                                            >  1:213701/1‑153 (MQ=255)
                                                                                                                                               cTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGa                                                                                                                            <  2:213701/153‑1 (MQ=255)
                                                                                                                                                         cTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCCATAT‑GCCGTCTTCTTGCGCGAAGGTTTCCACGTTGTCCCGCCGAACAAAAAGGCCAACACCTCGTCGATGGCTTACTACCATCTGtt                                                                                                             >  2:123112/1‑158 (MQ=255)
                                                                                                                                                         cTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGtt                                                                                                             <  1:123112/158‑1 (MQ=255)
                                                                                                                                                            gTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTggg                                                      >  1:201534/1‑210 (MQ=255)
                                                                                                                                                            gTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTggg                                                      <  2:201534/210‑1 (MQ=255)
                                                                                                                                                             tGAAAGAATATCATCTGCTGACCCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATCGCGGA‑TTCTTGCGCGAAGGTTTCCACGGTGTCACGCAGAACAAACAGGCCAACACCTCGTCGAAGGATTACTACCATATGTTGCGTCCTGCGGCGGAAAAATCGAGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAACAACTGCTCacg            >  2:264836/1‑249 (MQ=255)
                                                                                                                                                                        catcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGtgatg  <  1:176866/250‑1 (MQ=255)
                                                                                                                                                                                          aTTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCCCGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGtt                                                                                                             >  2:234208/1‑125 (MQ=255)
                                                                                                                                                                                          aTTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGtt                                                                                                             <  1:234208/125‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                    
AGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATAT‑GCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATG  >  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: ‑