Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,154 C→T 18.2% 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,1540CT18.2% 77.6 / 10.7 33D606D (GAC→GATthrAfused aspartokinase I and homoserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base C (9/18);  new base T (3/3);  total (12/21)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGA  >  W3110S.gb/2088‑2224
                                                                  |                                                                      
gtcgATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATAcc                                                                     <  1:3381117/70‑1 (MQ=255)
gtcgATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAcc                                                                     <  1:3728276/70‑1 (MQ=255)
   gATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACg                                                                 >  1:3345745/1‑71 (MQ=255)
       gATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGATACCAACGTTgg                                                             >  1:609236/1‑71 (MQ=255)
       gATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTg                                                              <  1:4506169/70‑1 (MQ=255)
            ctacCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTg                                                        <  1:1838536/71‑1 (MQ=255)
                cATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGAtt                                                    <  1:843735/71‑1 (MQ=255)
                cATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGAtt                                                    <  1:4481460/71‑1 (MQ=255)
                    aGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCg                                                >  1:954291/1‑71 (MQ=255)
                        gCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTa                                            >  1:1245457/1‑71 (MQ=255)
                        gCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTa                                            >  1:1368227/1‑71 (MQ=255)
                          gTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTAtt                                          <  1:17185/71‑1 (MQ=255)
                          gTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTAtt                                          <  1:147035/71‑1 (MQ=255)
                                cggcggAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAAc                                    >  1:3577094/1‑71 (MQ=255)
                                  gcggAAAAATCGCGGCGTAAATTCCTATAAGACACCAACg                                                                 <  1:483929/40‑1 (MQ=38)
                                   cggAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCg                                                >  1:730806/1‑56 (MQ=255)
                                     gAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTgag                                       <  1:2472555/63‑1 (MQ=255)
                                      aaaaaTCGCGGCGTAAATTCCTCTATGATACCAACGTTGGGGCTg                                                        <  1:886410/45‑1 (MQ=255)
                                       aaaaTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGt                                              <  1:3759386/54‑1 (MQ=255)
                                       aaaaTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTga                                        >  1:4074287/1‑60 (MQ=255)
                                       aaaaTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTga                                        >  1:1620676/1‑60 (MQ=255)
                                       aaaaTCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTga                                        >  1:2231810/1‑60 (MQ=255)
                                                 cgTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTa                                            <  1:3172134/46‑1 (MQ=255)
                                                 cgTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTa                                            <  1:1030457/46‑1 (MQ=255)
                                                 cgTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCaa                   <  1:3909705/71‑1 (MQ=255)
                                                  gTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTa                                            <  1:2209594/45‑1 (MQ=255)
                                                   tAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGtt                                             <  1:4330248/43‑1 (MQ=255)
                                                   tAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGt                                              <  1:3872416/42‑1 (MQ=255)
                                                     aaTTCCTCTATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCa               <  1:3130319/71‑1 (MQ=255)
                                                          ctctATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCaa                   >  1:2525160/1‑62 (MQ=255)
                                                           tctATGATACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCaa                   >  1:2413276/1‑61 (MQ=255)
                                                             tATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCaaa                             <  1:4664812/49‑1 (MQ=255)
                                                             tATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCaaa                             <  1:730463/49‑1 (MQ=255)
                                                                  cacCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATtga  <  1:520319/71‑1 (MQ=255)
                                                                  |                                                                      
GTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGA  >  W3110S.gb/2088‑2224

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: