Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,154 C→T 75.0% D606D (GAC→GAT thrA → fused aspartokinase I and homoserine dehydrogenase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,1540CT75.0% 31.6 / 9.6 20D606D (GAC→GATthrAfused aspartokinase I and homoserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (6/9);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

GATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATT  >  W3110S.gb/2095‑2222
                                                           |                                                                    
gATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTgg                                                           >  1:665104/1‑71 (MQ=255)
gATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTgg                                                           >  1:1512551/1‑71 (MQ=255)
     ctacCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTg                                                      <  1:1435546/71‑1 (MQ=255)
     ctacCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTg                                                      <  1:1834774/71‑1 (MQ=255)
            cAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTAcc                                               >  1:189678/1‑71 (MQ=255)
             aGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCg                                              >  1:1172617/1‑71 (MQ=255)
             aGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCg                                              >  1:521977/1‑71 (MQ=255)
                     tATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTAtt                                        >  1:1375774/1‑69 (MQ=255)
                                aaaaTCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCaaa                           <  1:1518702/71‑1 (MQ=255)
                                aaaaTCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCaaa                           <  1:660349/71‑1 (MQ=255)
                                   aTCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATc                        >  1:1650143/1‑71 (MQ=255)
                                      gcggcgTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTg                      <  1:909005/70‑1 (MQ=255)
                                             aaaTTCCTCTATGATACCAACGTTGGGGCTGGATTAcc                                               <  1:1683971/38‑1 (MQ=255)
                                             aaaTTCCTCTATGATACCAACGTTGGGGCTGGATTAcc                                               <  1:351168/38‑1 (MQ=255)
                                             aaaTTCCTCTATGATACCAACGTTGGGGCTGGATTAcc                                               <  1:543038/38‑1 (MQ=255)
                                              aaTTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCa             <  1:354368/71‑1 (MQ=255)
                                              aaTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATCGAGAACCTGCAAAATCTg                      >  1:294662/1‑62 (MQ=255)
                                              aaTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATCGAGAACCTGCAAAATCTg                      >  1:820216/1‑62 (MQ=255)
                                                  ccTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGt          <  1:597958/70‑1 (MQ=255)
                                                  ccTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGt          <  1:1832822/70‑1 (MQ=255)
                                                          acacCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAAtt  <  1:940051/70‑1 (MQ=255)
                                                           |                                                                    
GATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATT  >  W3110S.gb/2095‑2222

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