Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 260,015 C→T G29D (GGT→GAT)  phoE ← outer membrane porin PhoE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913260,0150CT100.0% 58.4 / NA 22G29D (GGT→GAT) phoEouter membrane porin PhoE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/11);  total (11/11)

CGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGTTTTTTGCTGTGCAGTATCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTACCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCGGGAGAGTCCCGTTGAAAACAGGAAAGTTTTTAACCTGAGATTGTTAAAGATATATTACAGATTAATAATATTCTTAAAATGTGGTAATTTATTAAATCTGTAATAAAAGCG  >  NC_000913/259770‑260250
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           
cGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGTTTTTTGCTGTGCAGTATCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCt                                                                                                                                                                                                                                        >  1:73733/1‑251 (MQ=255)
                                            ngcgTTTTTTGCTGTGCAGTATCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTAAAGATGCAGATGCCACAAt                                                                                                                                                                                            >  1:277725/2‑251 (MQ=255)
                                                                        ctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCAGAGTgn                                                                                                                                                                  <  1:78597/249‑2 (MQ=255)
                                                                                         aaCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCAGAGTGCTCTTTTTCATTTTCATTcc                                                                                                                                               <  1:339751/251‑1 (MQ=255)
                                                                                                         gACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCTTTATTATATATTTTTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAAcgc                                                                                                                               <  1:314229/251‑1 (MQ=255)
                                                                                                            cATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATCACTCTGGTCCCAATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACACTGCCCATCACCACTAATGCCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAAcgcgcg                                                                                                                            >  1:181867/1‑251 (MQ=255)
                                                                                                                                 tAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCGGGAGAGTcc                                                                                                        <  1:166590/250‑1 (MQ=255)
                                                                                                                                          tttCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCcac                                                                                                                                                                                               <  1:320667/154‑1 (MQ=255)
                                                                                                                                          tttCGCCTNTNNNANNAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACNCATATAATGCATGGCNTTNCCTTNGCCNTAGACATCCAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCcac                                                                                                                                                                                               >  2:320667/1‑154 (MQ=255)
                                                                                                                                              gCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCGGGAGAGTc                                                                                                         <  1:23404/236‑1 (MQ=255)
                                                                                                                                               ccTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCag                                                                                                                                                                       >  1:214998/1‑173 (MQ=255)
                                                                                                                                                       aCCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATATTGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCAGAGTGCTCTttttcattt                                                                                                                                                      <  1:403253/182‑1 (MQ=255)
                                                                                                                                                                     aaCTCTGGTCGCCATCTTTACGGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGGTTATTATCGTCGTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCGAGAGTGCTCTTATTCATTTTCATTGCTTATTTTAATGAACGCGCGAATATTCAGCGGGAGAGTTGCTTGGAAAACAGGAAAGTTTTTAACCTGAGATTGTTa                                                                   <  1:186393/251‑1 (MQ=255)
                                                                                                                                                                            ntCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCCGATGCCACAATGCCCa                                                                                                                                                                                       >  1:112530/2‑128 (MQ=255)
                                                                                                                                                                                                tATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCTTTATTATATATTTCTGGAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCAGAGTGCTCTGTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATAGTCAGCGGGAGAGTCCCGTTGAAAACAGGAAAGTTTTTAAcc                                                                              <  1:210962/213‑1 (MQ=255)
                                                                                                                                                                                                               gCATGGCTTTAACTTTGCCATCGACATACAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCAGAGTGCTCTttttcatt                                                                                                                                                       >  1:134050/1‑125 (MQ=255)
                                                                                                                                                                                                                   ggCTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCTTTCTTATATATTTCTGCAGCCTGTACAGCTGCAGCTGCCACACTGCCCa                                                                                                                                                                                       >  1:62514/1‑89 (MQ=255)
                                                                                                                                                                                                                    ncTTTAACTTTGCCATAGACATCCAGTTTATTATCGTCTTTATTATATCTTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCAGCGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCGGGAGAGTCCCGTTGAAAACAGGAAAGTTTTTAACCTGAGATTGTTAAAGATATATTACAGATTAATAATATTCTTAAAATGTGGTAATttatt                    >  1:40939/2‑251 (MQ=255)
                                                                                                                                                                                                                                      aCATCCAGGTTATTATCGTCTTTATTATATAGTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACGAATGCCAGAGTGCTCTGTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCGGGAGAGTCCCGTTGAAAACAGGAAACTGTTTAACCGGAGATTGTTAAAGAGATATTACAGAGCAATAATATGCTTAAAATGTGGTAATTTATTAAATCTGTAATAAAAGcn  <  1:44171/251‑2 (MQ=255)
                                                                                                                                                                                                                                       naTCCAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATCACACAATGCCCATCACCACTAAAGCCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGAGAATATTCAGCGGGAGAGTCCCGTTGAAAACAGGAAAGTTTTTAACCTGAGATTGTTAAAGATATATTACCGATTAATAATATTCTTAACATGTGGTAATTTATTAAATCTGTAATAAAAGCg  >  1:359599/2‑250 (MQ=255)
                                                                                                                                                                                                                                          ccAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCCGATGCCACACTGCCCATCCCCACTAATGCCCGAGTCCTCTTTTTCATTTTCATTCCTGCTTTTAATTAACGCGCGCATATTCAGCGGGagag                                                                                                           >  1:435290/1‑142 (MQ=255)
                                                                                                                                                                                                                                           cAGTTTATTATCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCAGAGTGCTCTGTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCGGGAGAGTCCCGTTGAAAACAGGAAAGTTTTTAAcc                                                                              <  1:138044/170‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           
CGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGTTTTTTGCTGTGCAGTATCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTGGCGTTATCACTCATATAATGCATGGCTTTAACTTTGCCATAGACATCCAGTTTATTACCGTCTTTATTATATATTTCTGCAGCCTGTACAGATGCAGATGCCACAATGCCCATCACCACTAATGCCAGAGTGCTCTTTTTCATTTTCATTCCTGATTTTAATTAACGCGCGAATATTCAGCGGGAGAGTCCCGTTGAAAACAGGAAAGTTTTTAACCTGAGATTGTTAAAGATATATTACAGATTAATAATATTCTTAAAATGTGGTAATTTATTAAATCTGTAATAAAAGCG  >  NC_000913/259770‑260250

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: