Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,647,049 C→T Q102* (CAA→TAA)  bcp → thiol peroxidase, thioredoxin‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,647,0490CT100.0% 223.7 / NA 67Q102* (CAA→TAA) bcpthiol peroxidase, thioredoxin‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/38);  total (29/38)

AGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAG  >  minE/1646999‑1647113
                                                  |                                                                
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGAGCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:1807451/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:266487/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:3851158/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:3850453/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:3733714/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:3709716/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:370007/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:3689639/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:3655974/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:2964720/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:4309843/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:2684316/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:2234873/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:2120425/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:756449/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:115157/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:1404284/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:182079/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:1519543/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:1711092/65‑1 (MQ=255)
agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:1696761/65‑1 (MQ=255)
 gctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:4080540/64‑1 (MQ=255)
  ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:1968860/1‑65 (MQ=255)
  ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:3723766/1‑65 (MQ=255)
  ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:1566976/1‑65 (MQ=255)
  ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:4077779/1‑65 (MQ=255)
  ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:3476754/1‑65 (MQ=255)
  ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:3443179/1‑65 (MQ=255)
  ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:1654293/1‑65 (MQ=255)
  ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:4344681/1‑65 (MQ=255)
  ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:582602/1‑65 (MQ=255)
  ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTAGGCGTCTgggg                                                  >  1:1418398/1‑65 (MQ=255)
         aCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCtt                                      <  1:4088827/70‑1 (MQ=255)
         aCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCtt                                      <  1:4323257/70‑1 (MQ=255)
         aCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCtt                                      <  1:1619927/70‑1 (MQ=255)
         aCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCtt                                      <  1:2383450/70‑1 (MQ=255)
         aCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCtt                                      <  1:812331/70‑1 (MQ=255)
         aCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCtt                                      <  1:140293/70‑1 (MQ=255)
         aCTTTACGCTCCTGTCTGATGAGGACCACCAGGGGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCtt                                      <  1:2900827/70‑1 (MQ=255)
                    cTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAAc                          <  1:4234048/71‑1 (MQ=255)
                    cTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAAc                          <  1:3437399/71‑1 (MQ=255)
                                gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATgg                 >  1:812615/1‑68 (MQ=255)
                                gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATgg                 >  1:576513/1‑68 (MQ=255)
                                gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATgg                 >  1:3947111/1‑68 (MQ=255)
                                gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATgg                 >  1:2775967/1‑68 (MQ=255)
                                gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATgg                 >  1:1239491/1‑68 (MQ=255)
                                gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATgg                 >  1:1829055/1‑68 (MQ=255)
                                gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACCATgg                 >  1:2991157/1‑68 (MQ=255)
                                   caccaGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAAc                          >  1:2911917/1‑56 (MQ=255)
                                    accaGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCAt          <  1:3956866/71‑1 (MQ=255)
                                    accaGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCAt          <  1:1897528/71‑1 (MQ=255)
                                    accaGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCAt          <  1:2433377/71‑1 (MQ=255)
                                       aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:3999313/1‑65 (MQ=255)
                                       aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:653749/1‑65 (MQ=255)
                                       aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:641857/1‑65 (MQ=255)
                                       aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:4377201/1‑65 (MQ=255)
                                       aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:2548188/1‑65 (MQ=255)
                                       aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:3924825/1‑65 (MQ=255)
                                       aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:4283407/1‑65 (MQ=255)
                                       aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:3525863/1‑65 (MQ=255)
                                       aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:3985278/1‑65 (MQ=255)
                                            tgCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAg  <  1:455203/71‑1 (MQ=255)
                                            tgCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAg  <  1:204531/71‑1 (MQ=255)
                                            tgCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAg  >  1:3077129/1‑71 (MQ=255)
                                            tgCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAg  >  1:720182/1‑71 (MQ=255)
                                            tgCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAg  <  1:378254/71‑1 (MQ=255)
                                             gCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAg  <  1:3139324/70‑1 (MQ=255)
                                                  |                                                                
AGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAG  >  minE/1646999‑1647113

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: