Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,154 C→T 89.8% 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,1540CT89.8% 98.4 / 4.9 39D606D (GAC→GATthrAfused aspartokinase I and homoserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (12/23);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAAT  >  W3110S.gb/2086‑2221
                                                                    |                                                                   
tcgtcgATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGatac                                                                   >  1:932613/1‑71 (MQ=255)
       tGGATTACTACCATCTGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGtt                                                            <  1:2031130/71‑1 (MQ=255)
            tactacCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGc                                                       >  1:2553549/1‑71 (MQ=255)
            tactacCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGc                                                       >  1:728067/1‑71 (MQ=255)
            tactacCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGc                                                       >  1:2008039/1‑71 (MQ=255)
                   aTCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTa                                                >  1:2494946/1‑71 (MQ=255)
                      aGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTgg                                                          >  1:147933/1‑58 (MQ=255)
                      aGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTgg                                                          >  1:1489457/1‑58 (MQ=255)
                      aGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTgg                                                          >  1:1584193/1‑58 (MQ=255)
                               aTGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTAtt                                       >  1:202657/1‑68 (MQ=255)
                                tGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTgag                                    <  1:332350/70‑1 (MQ=255)
                                          aaaTCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTgaga                                   <  1:1628583/61‑1 (MQ=255)
                                          aaaTCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTgaga                                   <  1:1083176/61‑1 (MQ=255)
                                          aaaTCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTgaga                                   <  1:267187/61‑1 (MQ=255)
                                           aaTCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTgaga                                   <  1:1408259/60‑1 (MQ=255)
                                             tCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGa                                                   >  1:65494/1‑42 (MQ=255)
                                             tCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGa                                                   >  1:1890258/1‑42 (MQ=255)
                                                 ggcgTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTc                  <  1:1945072/71‑1 (MQ=255)
                                                   cgTAAATTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTgaga                                   <  1:571550/52‑1 (MQ=255)
                                                     tAAATTCCTCTATGATACCATCGTTGGGGCTGGATTAcc                                              <  1:703279/39‑1 (MQ=38)
                                                     tAAATTCCTCTATGATACCAACGTTGGGGCTGGATTAcc                                              <  1:712070/39‑1 (MQ=255)
                                                     tAAATTCCTCTATGATACCAACGTTGGGGCTGGATTAcc                                              <  1:1226238/39‑1 (MQ=255)
                                                     tAAATTCCTCTATGATACCAACGTTGGGGCTGGATTAcc                                              <  1:534148/39‑1 (MQ=255)
                                                       aaTTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAAt               >  1:1830145/1‑68 (MQ=255)
                                                       aaTTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATg              >  1:532831/1‑69 (MQ=255)
                                                       aaTTCCTCTATGATACAAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATg              >  1:421283/1‑69 (MQ=255)
                                                           ccTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGt         <  1:821094/70‑1 (MQ=255)
                                                           ccTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGt         <  1:2399871/70‑1 (MQ=255)
                                                           ccTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGt         <  1:1958602/70‑1 (MQ=255)
                                                            ctctATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGtg        <  1:1216489/70‑1 (MQ=255)
                                                                 tGATACCAACGTTGGGGCTGGATTACCGGTTATTTAGAACCTGCAAAATCTGCTCAATGc             <  1:83424/60‑1 (MQ=255)
                                                                 tGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTGAATGc             <  1:1125440/60‑1 (MQ=255)
                                                                 tGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGc             <  1:312187/60‑1 (MQ=255)
                                                                 tGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGc             <  1:173858/60‑1 (MQ=255)
                                                                 tGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGtgatg     <  1:293822/68‑1 (MQ=255)
                                                                 tGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGtgatg     <  1:2590468/68‑1 (MQ=255)
                                                                 tGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGtgatg     <  1:1315476/68‑1 (MQ=255)
                                                                 tGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGtgatg     <  1:842643/68‑1 (MQ=255)
                                                                 tGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAAt  >  1:1842235/1‑71 (MQ=255)
                                                                 tGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGtgatg     <  1:1656013/68‑1 (MQ=255)
                                                                    |                                                                   
TCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAAT  >  W3110S.gb/2086‑2221

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