Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 259,305 G→A A7V (GCA→GTA)  phoE ← outer membrane phosphoporin protein E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb259,3050GA100.0% 105.3 / NA 34A7V (GCA→GTA) phoEouter membrane phosphoporin protein E
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/19);  total (15/19)

CCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCGG  >  W3110S.gb/259239‑259364
                                                                  |                                                           
ccGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATACCa                                                          <  1:1903274/70‑1 (MQ=255)
ccGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATACCa                                                          <  1:2045310/70‑1 (MQ=255)
            tatatTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATACCa                                                          <  1:86932/58‑1 (MQ=255)
                  tCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATACCAGAGTGCTCTttttcatttt                                       >  1:3533899/1‑71 (MQ=255)
                      cAGCCTTTACAGATGCAGATGCCACAATGCCCATCACCACTAATACCAGAGTGCTCTTTttcattttcatt                                   <  1:2048101/71‑1 (MQ=255)
                      cAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATACCAGAGTGCTCTTTttcattttcatt                                   <  1:4114963/71‑1 (MQ=255)
                      cAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATACCAGAGTGCTCTTTttcattttcatt                                   <  1:3181042/71‑1 (MQ=255)
                      cAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATACCAGAGTGCTCTTTttcattttcatt                                   <  1:3118080/71‑1 (MQ=255)
                      cAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATACCAGAGTGCTCTTTttcattttcatt                                   <  1:3105089/71‑1 (MQ=255)
                      cAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATACCAGAGTGCTCTTTttcattttcatt                                   <  1:2829860/71‑1 (MQ=255)
                      cAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATACCAGAGTGCTCTTTttcattttcatt                                   <  1:2538374/71‑1 (MQ=255)
                       aGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATACCAGAGTGCTCTTTttcattttcatt                                   <  1:1031950/70‑1 (MQ=255)
                                agatgcagatgcCACAATGCCCATCACCACTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTtta                         <  1:3545468/71‑1 (MQ=255)
                                 gatgcagatgcCACAATGCCCATCACCACTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTtta                         <  1:22445/70‑1 (MQ=255)
                                       gatgcCACAATGCCCATCACCACTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAAcg                  >  1:2888434/1‑71 (MQ=255)
                                       gatgcCACAATGCCCATCACCACTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAAcg                  >  1:2425051/1‑71 (MQ=255)
                                       gatgcCACAATGCCCATCACCACTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAAcg                  >  1:1877406/1‑71 (MQ=255)
                                         tgcCACAATGCCCATCACCACTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAAcgcg                >  1:3811599/1‑71 (MQ=255)
                                                aTGCCCATCACCACTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAatat         <  1:3504329/71‑1 (MQ=255)
                                                aTGCCCATCACCACTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAatat         <  1:3508996/71‑1 (MQ=255)
                                                 tGCCCATCACCACTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGaa            <  1:4123478/67‑1 (MQ=255)
                                                 tGCCCATCACCACTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGaa            <  1:178262/67‑1 (MQ=255)
                                                 tGCCCATCACCACTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGaa            <  1:1260830/67‑1 (MQ=255)
                                                 tGCCCATCACCACTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGaa            <  1:367909/67‑1 (MQ=255)
                                                   cccATCACCACTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAata          >  1:1849925/1‑67 (MQ=255)
                                                         accacTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCgg  >  1:1041638/1‑69 (MQ=255)
                                                         accacTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCgg  >  1:51094/1‑69 (MQ=255)
                                                         accacTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCgg  >  1:116626/1‑69 (MQ=255)
                                                         accacTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCgg  >  1:3837109/1‑69 (MQ=255)
                                                         accacTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCgg  >  1:3645087/1‑69 (MQ=255)
                                                         accacTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCgg  >  1:3623016/1‑69 (MQ=255)
                                                         accacTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCgg  >  1:1852826/1‑69 (MQ=255)
                                                         accacTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCgg  >  1:2223325/1‑69 (MQ=255)
                                                         accacTAATACCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCgg  >  1:228552/1‑69 (MQ=255)
                                                                  |                                                           
CCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCGG  >  W3110S.gb/259239‑259364

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: