Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,189,993 A→G E38G (GAG→GGG)  yqiH → putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,189,9930AG100.0% 103.2 / NA 28E38G (GAG→GGG) yqiHputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/13);  total (15/13)

AAAATTAATTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAA  >  NC_000913/3189858‑3190118
                                                                                                                                       |                                                                                                                             
aaaaTTAATTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTc                                                                                                                            >  1:611304/1‑139 (MQ=255)
                                 tAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCa                                                                                           <  2:611304/139‑1 (MQ=255)
                                                      cagcagGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCt                                                                      >  1:265750/1‑139 (MQ=255)
                                                          agGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGcc                                                                  <  2:502886/139‑1 (MQ=255)
                                                          agGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGcc                                                                  <  1:66952/139‑1 (MQ=255)
                                                                    ctctCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGa                                                        <  2:340032/139‑1 (MQ=255)
                                                                       tcCTGTATGAGTAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATag                                                     >  1:621630/1‑139 (MQ=255)
                                                                       tcCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATag                                                     >  2:284193/1‑139 (MQ=255)
                                                                       tcCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATag                                                     <  2:265750/139‑1 (MQ=255)
                                                                                 gcagcaTTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCa                                                                                           >  1:332962/1‑91 (MQ=255)
                                                                                 gcagcaTTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCa                                                                                           <  2:332962/91‑1 (MQ=255)
                                                                                   agcaTTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAg                                         >  1:523090/1‑139 (MQ=255)
                                                                                     caTTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGc                                       >  2:277679/1‑139 (MQ=255)
                                                                                      aTTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTaaaaa                                                                                  <  1:108209/95‑1 (MQ=255)
                                                                                      aTTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTaaaaa                                                                                  >  2:108209/1‑95 (MQ=255)
                                                                                                    catcatCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCaaa                                                              <  2:435252/101‑1 (MQ=255)
                                                                                                    catcatCCCGGATCGTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATAAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCaaa                                                              >  1:435252/1‑101 (MQ=255)
                                                                                                                  gTACCCGTATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGt          <  1:284193/139‑1 (MQ=255)
                                                                                                                         tATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGcc                                                                  <  1:527374/76‑1 (MQ=255)
                                                                                                                         tATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGcc                                                                  >  2:527374/1‑76 (MQ=255)
                                                                                                                         tATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCa   >  1:604655/1‑139 (MQ=255)
                                                                                                                         tATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCa   >  2:54665/1‑139 (MQ=255)
                                                                                                                         tATCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCa   <  2:549864/139‑1 (MQ=255)
                                                                                                                          aTCATAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCaa  >  2:535832/1‑139 (MQ=255)
                                                                                                                             aTAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAGAAACCTGCCCTATCTGGCCCAAACa                                                            <  2:479595/78‑1 (MQ=255)
                                                                                                                             aTAATGAATGGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTTGCCCAAACa                                                            >  1:479595/1‑78 (MQ=255)
                                                                                                                               aatgaatgGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCa                                                                                           <  1:564919/45‑1 (MQ=255)
                                                                                                                               aatgaatgGGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCa                                                                                           >  2:564919/1‑45 (MQ=255)
                                                                                                                                       |                                                                                                                             
AAAATTAATTCAGGATGAAACTAATGCGCTACTTAAATACAAAGAATATCATTGCAGCAGGAGTATTACTCTCCTGTATGAGCAGCATTGCCTGGGGTGCCATCATCCCGGATCGTACCCGTATCATAATGAATGAGTCAGATAAAGGCGAAGCATTAAAACTTACTAATCAAAGTAAAAACCTGCCCTATCTGGCCCAAACATGGATAGAAGATACGAAAGGCAATAAATCGCGTGACTTTATTGTCACTGTTCCTCCAA  >  NC_000913/3189858‑3190118

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