Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 589,167 T→C E596E (GAA→GAG nfrA ← bacteriophage N4 receptor, outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913589,1670TC100.0% 40.1 / NA 12E596E (GAA→GAGnfrAbacteriophage N4 receptor, outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (7/5);  total (7/5)

GCGCGCAAATCACTCACCGCGGCCGGGACATTATGACGTTGGCGATAAATTGTCGCCCGCGCAACGTAAGCGTTGGCAGAAGGCGCAATATTGATTGAGCGCGTGAGATCGTTCAGTGCGAGTTCCGGCTGACCAGGAATGTAACGTTGCGCATGCAGCCACCAGTAGAGGGCATTGCTTCCCAGTCCACGTTTTTCTGCCTGTTGTAGCCAGCGATCGCGAGCCGCACCATTTCCTGCCGCCTGGGCGGTATTGGCA  >  NC_000913/589045‑589302
                                                                                                                          |                                                                                                                                       
gcgcgcAAATCACTCACCGCGGCCGGGACATTATGACGTTGGCGATAAATTGTCGCCCGCGCAACGTAAGCGTTGGCAGAAGGCGCAATATTGATTGAGCGCGTGAGATCGTTCAGTGCGAGCTCCGGCTGACCAGGAAt                                                                                                                        >  2:32901/1‑140 (MQ=255)
     cAAATCACTCACCGCGGCCGGGACATTATGACGTTGGCGATAAATTGTCGCCCGCGCAACGTAAGCGTTGGCAGAAGGCGCAATATTGATTGAGCGCGTGAGATCGTTCAGTGCGAGCTCCGGCTGACCAGGAATGTAAc                                                                                                                   <  1:32901/140‑1 (MQ=255)
            ctcacCGCGGCCGGGACATTATGACGTTGGCGATAAATTGTCGCCCGCGCAACGTAAGCGTTGGCAGAAGGCGCAATATTGATTGAGCGCGTGAGATCGTTCAGTGCGAGCTCCGGCTGACCAGGAATGTAACGTTgcgc                                                                                                            <  2:129500/140‑1 (MQ=255)
                 cgcgGCCGGGACATTATGACGTTGGCGATAAATTGTCGCCCGCGCAACGTAAGCGTTGGCAGAAGGCGCAATATTGATTGAGCGCGTGAGATCGTTCAGTGCGAGCTCCGGCTGACCAGGAATGTAACGTTGCGCATGCa                                                                                                       <  2:310271/140‑1 (MQ=255)
                  gcgGCCGGGACATTATGACGTTGGCGATAAATTGTCGCCCGCGCAACGTAAGCGTTGGCAGAAGGCGCAATATTGATTGGGCGCGTGAGATCGTTCAGTGCGAGCTCCGGCTGACCAGGAATGTAACGTTGCGCATGCAg                                                                                                      >  2:375300/1‑140 (MQ=255)
                   cgGCCGGGACATTATGACGTTGGCGATAAATTGTCGCCCGCGCAACGTAAGCGTTGGCAGAAGGCGCAATATTGATTGAGCGCGTGAGATCGTTCAGTGCGAGCTCCGGCTGACCAGGAATGTAACGTTGCGCATGCAGc                                                                                                     <  1:276606/140‑1 (MQ=255)
                                                                                         attgattgAGCGCGTGAGATCGTTCAGTGCGAGCTCCGGCTGACCAGGAATGTAACGTTGCGCATGCAGCCACCAGTAGAGGGCATTGCTTCCCAGTCCACGTTTTTCTGCCTGTTGTAGCCAGCGATCGCGAGCCGCAc                               >  1:415168/1‑140 (MQ=255)
                                                                                            gattgaGCGCGTGAGATCGTTCAGTGCGAGCTCCGGCTGACCAGGAATGTAACGTTGCGCATGCAGCCACCAGTAGAGGGCATTGCTTCCCAGTCCACGTTTTTCTGCCTGTTGTAGCCAGCGATCGCGAGCCGCACCAt                            <  2:415168/140‑1 (MQ=255)
                                                                                                  gcgcGTGAGATCGTTCAGTGCGAGCTCCGGCTGACCAGGAATGTAACGTTGCGCATGCAGCCACCAGTAGAGGGCATTGCTTCCCAGTCCACGTTTTTCTGCCTGTTGTAGCCAGCGATCGCGAGCCGCACCATTTCCTg                      >  2:304502/1‑140 (MQ=255)
                                                                                                                    tGCGAGCTCCGGCTGACCAGGAATGTAACGTTGCGCATGCAGCCACCAGTAGAGGGCATTGCTTCCCAGTCCACGTTTTTCTGCCTGTTGTAGCCAGCGATCGCGAGCCGCACCATTTCCTGCCGCCTGGGCGGTATTgg    >  1:411376/1‑140 (MQ=255)
                                                                                                                      cGAGCTCCGGCTGACCAGGAATGTAACGTTGCGCATGCAGCCACCAGTAGAGGGCATTGCTTCCCAGTCCACGTTTTTCTGCCTGTTGTAGCCAGCGATCGCGAGCCGCACCATTTCCTGCCGCCTGGGCGGTATTGgca  >  1:399853/1‑140 (MQ=255)
                                                                                                                      cGAGCTCCGGCTGACCAGGAATGTAACGTTGCGCATGCAGCCACCAGTAGAGGGCATTGCTTCCCAGTCCACGTTTTTCTGCCTGTTGTAGCCAGCGATCGCGAGCCGCACCATTTCCTGCCGCCTGGGCGGTATTGgca  >  2:80430/1‑140 (MQ=255)
                                                                                                                          |                                                                                                                                       
GCGCGCAAATCACTCACCGCGGCCGGGACATTATGACGTTGGCGATAAATTGTCGCCCGCGCAACGTAAGCGTTGGCAGAAGGCGCAATATTGATTGAGCGCGTGAGATCGTTCAGTGCGAGTTCCGGCTGACCAGGAATGTAACGTTGCGCATGCAGCCACCAGTAGAGGGCATTGCTTCCCAGTCCACGTTTTTCTGCCTGTTGTAGCCAGCGATCGCGAGCCGCACCATTTCCTGCCGCCTGGGCGGTATTGGCA  >  NC_000913/589045‑589302

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: