Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,265,554 T→C 7.4% P110P (CCA→CCG tdcB ← catabolic threonine dehydratase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,265,5540TC 7.4% 141.5 / 3.4 54P110P (CCA→CCGtdcBcatabolic threonine dehydratase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (17/33);  new base C (2/2);  total (19/35)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATCACCATGCAGAACGACTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATCGCGCAGGAGAGGGAAACCCCTTG  >  W3110S.gb/3265495‑3265622
                                                           |                                                                    
ttATCACCATGCAGAACGACTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCttt                                                           <  1:913307/71‑1 (MQ=255)
 tATCACCATGCAGAACGACTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTCGGCGCACCtttt                                                          >  1:3524609/1‑71 (MQ=255)
  aTCACCATGCAGAACGACTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTg                                                                     >  1:2437877/1‑59 (MQ=255)
  aTCACCATGCAGAACGACTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTg                                                         >  1:3002016/1‑71 (MQ=255)
   tCACCATGCAGAACGACTTCTGCGGAGTAGTCGCAGGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTgg                                                        <  1:1472159/71‑1 (MQ=255)
         tGCAGAACGACTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATcac                                                  <  1:2956812/71‑1 (MQ=255)
          gCAGAACGACTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCAcc                                                              >  1:1486848/1‑58 (MQ=255)
           cAGAACGACTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATcacca                                                <  1:1631435/71‑1 (MQ=255)
             gAACGACTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCAc                                                               <  1:2239338/54‑1 (MQ=255)
             gAACGACTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATcac                                                  <  1:814540/67‑1 (MQ=255)
               aCGACTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACttt                                            <  1:1585581/71‑1 (MQ=255)
               aCGACTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACtt                                             >  1:2310464/1‑70 (MQ=255)
                  aCTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTa                                           >  1:2066432/1‑69 (MQ=255)
                   cTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCtt                                                            <  1:1470028/51‑1 (MQ=255)
                   cTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACtt                                             <  1:2208733/66‑1 (MQ=255)
                    ttCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATcaccac                                               <  1:294356/63‑1 (MQ=255)
                    ttCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGt                                       >  1:2398106/1‑71 (MQ=255)
                      cTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTc                                      >  1:41749/1‑70 (MQ=255)
                      cTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTc                                      >  1:1488741/1‑70 (MQ=255)
                            agtagtCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATAcc                                >  1:3382432/1‑70 (MQ=255)
                              tagtCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTAc                                          <  1:303071/58‑1 (MQ=255)
                              tagtCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAg                             <  1:971796/71‑1 (MQ=255)
                                gtCGCACGTTGCCGCTACTTTGGATTTTGGCGCAcc                                                              >  1:3056745/1‑36 (MQ=255)
                                gtCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCttt                                                           <  1:3711142/39‑1 (MQ=255)
                                gtCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGc                                                       <  1:1306718/43‑1 (MQ=255)
                                  cGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACt                                              >  1:2810834/1‑50 (MQ=255)
                                  cGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACGCAGCATc                         >  1:3271425/1‑71 (MQ=255)
                                  agCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACt                                              >  1:573733/2‑50 (MQ=255)
                                     aCGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATcgcg                      <  1:1212495/71‑1 (MQ=255)
                                       gTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTc                                      <  1:1689214/53‑1 (MQ=255)
                                        ttGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCa                           <  1:3051061/63‑1 (MQ=255)
                                        ttGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCAATACCCAGCATCGCGCa                    <  1:3202345/70‑1 (MQ=255)
                                           ccGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATCGCGCAgta                 <  1:980921/70‑3 (MQ=255)
                                           ccGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATCGCGCAGgag                <  1:2618763/71‑1 (MQ=255)
                                           ccGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATCGCGCAGga                 <  1:3704655/70‑1 (MQ=255)
                                            cGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATCGCGCAGga                 <  1:2873570/69‑1 (MQ=255)
                                            cGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATCGCGCAGga                 <  1:1352040/69‑1 (MQ=255)
                                             gCTACTTTGGATTTTGGCGCACCTTTTGGCATcac                                                  <  1:718183/35‑1 (MQ=255)
                                             gCTACTTTGGATTTTGGCGCACCTTTTGGCATcac                                                  <  1:220614/35‑1 (MQ=255)
                                             gCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATCGCGCAGgag                <  1:197230/69‑1 (MQ=255)
                                              cTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATa                                  >  1:2131527/1‑50 (MQ=255)
                                              cTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATcgcg                      <  1:2447126/62‑1 (MQ=255)
                                              cTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATCGCGCAg                   >  1:1938654/1‑65 (MQ=255)
                                               tACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATCGCGCAGGAGAgg             >  1:1330786/1‑70 (MQ=255)
                                                aCTTTGGATTTCGGCGCACCTTTTGGCATCACCACTTTACCGTCAATACCCAGCATcgc                       >  1:301626/1‑59 (MQ=255)
                                                  tttGGATTTTGGCGCACCTTTTGGCATCACCACttt                                            <  1:1872728/36‑1 (MQ=255)
                                                  tttGGATTTTGGCGCACCTTTTGGCATCACCACttt                                            <  1:2794106/36‑1 (MQ=255)
                                                  tttGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATCGCGCAGGagag              <  1:822109/66‑1 (MQ=255)
                                                  tttGGATTTCGGCGCACCTTTTGGCATCACCACTTTACCGTCAATACCCAGCATCgcgc                     <  1:2195699/59‑1 (MQ=255)
                                                  tttGGATTTCGGCGCACCTTTTGGCATCACCACTTTACCGTCAATACCCAGCATCgcgc                     <  1:748806/59‑1 (MQ=255)
                                                    tGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATAc                                 >  1:922032/1‑45 (MQ=255)
                                                      gATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATCGCGCAGGAGAGGGAAAcccc     <  1:3397013/71‑1 (MQ=255)
                                                       aTTTTGGCGCACCTTTTGGCATCACCACTTTACCGTc                                      <  1:3032357/37‑1 (MQ=255)
                                                       aTTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATCGCGCAGga                 <  1:1309050/58‑1 (MQ=255)
                                                         tttGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATCGCGCAGGAGAGGGAAACCCCTTg  >  1:1397739/1‑71 (MQ=255)
                                                           |                                                                    
TTATCACCATGCAGAACGACTTCTGCGGAGTAGTCGCACGTTGCCGCTACTTTGGATTTTGGCGCACCTTTTGGCATCACCACTTTACCGTCGATACCCAGCATCGCGCAGGAGAGGGAAACCCCTTG  >  W3110S.gb/3265495‑3265622

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: