Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,459,294 A→G E319G (GAA→GGA)  gspE → Type II secretion system protein GspE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,459,2940AG100.0% 76.7 / NA 22E319G (GAA→GGA) gspEType II secretion system protein GspE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/9);  total (13/9)

CGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCAT  >  NC_000913/3459162‑3459425
                                                                                                                                    |                                                                                                                                   
cGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGt                                                                                                                              <  2:176104/140‑1 (MQ=255)
  gTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGa                                                                                                                            <  1:382532/140‑1 (MQ=255)
    agaagaTCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATa                                                                                                                          >  1:154247/1‑140 (MQ=255)
                                              gCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGcc                                                                                <  2:205732/140‑1 (MQ=255)
                                                            gtgtGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGcc                                                                                >  1:188487/1‑126 (MQ=255)
                                                            gtgtGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGcc                                                                                <  2:188487/126‑1 (MQ=255)
                                                            gtgtGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCa                                                                  >  2:302684/1‑140 (MQ=255)
                                                                 gATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACt                                                             >  2:79240/1‑140 (MQ=255)
                                                                           tcgctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAAc                                                   >  2:9859/1‑140 (MQ=255)
                                                                            cgctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACa                                                  <  1:9859/140‑1 (MQ=255)
                                                                                cgcgGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGc                                              <  1:302684/140‑1 (MQ=255)
                                                                                 gcgGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCa                                             >  1:254932/1‑140 (MQ=255)
                                                                                       tGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGgc                                       >  2:102575/1‑140 (MQ=255)
                                                                                         cgcgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGc                                              >  2:154598/1‑131 (MQ=255)
                                                                                         cgcgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGc                                              <  1:154598/131‑1 (MQ=255)
                                                                                           cgcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAg                                   <  1:79240/140‑1 (MQ=255)
                                                                                                        cAGGACCCGGATGTCGTCATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCcg                      >  2:338622/1‑140 (MQ=255)
                                                                                                                     tcgtcATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTc         >  1:374110/1‑140 (MQ=255)
                                                                                                                      cgtcATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTGCTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTc                        <  1:351364/124‑1 (MQ=255)
                                                                                                                      cgtcATGGTGGGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTc                        >  2:351364/1‑124 (MQ=255)
                                                                                                                            ggtggGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCAt  >  2:250927/1‑140 (MQ=255)
                                                                                                                            ggtggGGGGAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCAt  >  1:239449/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                   
CGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCAT  >  NC_000913/3459162‑3459425

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