Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 52,852 G→A H189Y (CAT→TAT)  pdxA ← 4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb52,8520GA100.0% 113.8 / NA 37H189Y (CAT→TAT) pdxA4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/18);  total (19/18)

GGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCT  >  W3110S.gb/52788‑52898
                                                                |                                              
ggATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAAt                                          <  1:840977/71‑1 (MQ=255)
ggATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAAt                                          <  1:2270895/71‑1 (MQ=255)
ggATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAAt                                          <  1:2321924/71‑1 (MQ=255)
 gATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAACACCAAATTTGGTCCGCAAATCGTGATACAAAAt                                          <  1:3448854/70‑1 (MQ=255)
     cAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCa                                      <  1:3545968/70‑1 (MQ=255)
     cAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCa                                      <  1:1705660/70‑1 (MQ=255)
     cAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCa                                      <  1:2322136/70‑1 (MQ=255)
            cAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTc                              <  1:982575/71‑1 (MQ=255)
            cAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTc                              <  1:3564011/71‑1 (MQ=255)
            cAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTc                              <  1:3341144/71‑1 (MQ=255)
            cAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTc                              <  1:3178505/71‑1 (MQ=255)
            cAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTc                              <  1:3102013/71‑1 (MQ=255)
             aGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTc                              <  1:2538078/70‑1 (MQ=255)
                        cgcgGTTCGGCAATACCAAATTTGGTCCGCAAATCGTCATACAAAAt                                          <  1:3795903/47‑1 (MQ=39)
                               cGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGt            >  1:4076861/1‑70 (MQ=255)
                               cGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGt            >  1:2160510/1‑70 (MQ=255)
                               cGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGt            >  1:2133509/1‑70 (MQ=255)
                               cGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGt            >  1:1883512/1‑70 (MQ=255)
                               cGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGt            >  1:887695/1‑70 (MQ=255)
                               cGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGt            >  1:2980354/1‑70 (MQ=255)
                               cGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGt            >  1:614170/1‑70 (MQ=255)
                               cGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGt            >  1:3828858/1‑70 (MQ=255)
                               cGGCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGt            >  1:3928736/1‑70 (MQ=255)
                                ggCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTg           <  1:3400925/70‑1 (MQ=255)
                                ggCAATACCAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCCTGCAAAAGTGCAGGGGTg           <  1:651990/70‑1 (MQ=255)
                                       ccAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCa                 >  1:3850707/1‑57 (MQ=255)
                                       ccAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGt    >  1:3716575/1‑70 (MQ=255)
                                       ccAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTc   >  1:3421535/1‑71 (MQ=255)
                                       ccAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTc   >  1:420924/1‑71 (MQ=255)
                                       ccAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTc   >  1:818728/1‑71 (MQ=255)
                                       ccAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTc   >  1:3970455/1‑71 (MQ=255)
                                       ccAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTc   >  1:3710859/1‑71 (MQ=255)
                                       ccAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTc   >  1:107133/1‑71 (MQ=255)
                                       ccAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTc   >  1:3009308/1‑71 (MQ=255)
                                       ccAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTc   >  1:2111116/1‑71 (MQ=255)
                                        cAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCt  <  1:1566401/71‑1 (MQ=255)
                                        cAAATTTGGTCCGCAAATCGTGATACAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCt  <  1:1424689/71‑1 (MQ=255)
                                                                |                                              
GGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCT  >  W3110S.gb/52788‑52898

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