Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	306,891	A→G	G814G (GGT→GGC)	ecpC ←	putative fimbrial usher protein EcpC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	306,891	0	A	G	100.0%	130.3 / NA	38	G814G (GGT→GGC)	ecpC	putative fimbrial usher protein EcpC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/17);  total (21/17)

TGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATAT  >  NC_000913/306754‑307024                                                                                                                                          |                                                                                                                                      tGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGcc                                                                                                                                     >  2:487702/1‑140 (MQ=255)   aGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGcgccgc                                                                                                                                   >  2:209492/1‑140 (MQ=255)        aGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGcc                                                                                                                              <  1:209492/140‑1 (MQ=255)        aGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGcc                                                                                                                              >  2:342814/1‑140 (MQ=255)           tACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTgg                                                                                                                           >  2:173978/1‑140 (MQ=255)                gCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAg                                                                                                                      <  2:214652/140‑1 (MQ=255)                    gATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTcc                                                                                                                  >  1:377801/1‑140 (MQ=255)                     aTTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCa                                                                                                                 >  2:658469/1‑140 (MQ=255)                        aGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCagag                                                                                                              <  2:785442/140‑1 (MQ=255)                               aaaTTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTc                                                                                                       >  1:754188/1‑140 (MQ=255)                                    atcatcGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTcc                                                                                                                  >  1:284300/1‑124 (MQ=255)                                    atcatcGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTcc                                                                                                                  <  2:284300/124‑1 (MQ=255)                                                   agGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGt                                                                                   >  2:288980/1‑140 (MQ=255)                                                    ggTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTa                                                                                  >  1:432441/1‑140 (MQ=255)                                                        aaCTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGa                                                                              >  1:678348/1‑140 (MQ=255)                                                         aCTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGaa                                                                             >  2:795023/1‑140 (MQ=255)                                                               aTTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGa                                                                       <  2:432441/140‑1 (MQ=255)                                                                  ctcttctcCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAg                                                                    <  1:658469/140‑1 (MQ=255)                                                                  ctcttctcCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAg                                                                    <  2:678348/140‑1 (MQ=255)                                                                       ctcCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTgggg                                                               >  2:389954/1‑140 (MQ=255)                                                                       ctcCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTgggg                                                               >  2:804630/1‑140 (MQ=255)                                                                        tcCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGt                                                              <  1:342814/140‑1 (MQ=255)                                                                               tGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCtt                                                       >  2:10777/1‑140 (MQ=255)                                                                               tGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCtt                                                       >  1:208320/1‑140 (MQ=255)                                                                                 cGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTAt                                                     >  1:263675/1‑140 (MQ=255)                                                                                  gTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATc                                                    >  2:716152/1‑140 (MQ=255)                                                                                   tCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATcc                                                   <  1:804630/140‑1 (MQ=255)                                                                                       cgccgcCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGt                                               >  2:321001/1‑140 (MQ=255)                                                                                                     agagGCCGCTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCg                                 <  1:10777/140‑1 (MQ=255)                                                                                                             cTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTc                         <  2:263675/140‑1 (MQ=255)                                                                                                             cTGCAGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTc                         >  2:651089/1‑140 (MQ=255)                                                                                                                 aGACCACATCACCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCg                     >  2:485994/1‑140 (MQ=255)                                                                                                                           aCCGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGTAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCggttcgggttc           <  2:754188/140‑1 (MQ=255)                                                                                                                            ccGACCCAGACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCggttcgggttcg          >  1:152700/1‑140 (MQ=255)                                                                                                                                    gACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGatat  <  1:389954/140‑1 (MQ=255)                                                                                                                                    gACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGatat  <  1:487702/140‑1 (MQ=255)                                                                                                                                    gACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGatat  <  1:651089/140‑1 (MQ=255)                                                                                                                                    gACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGatat  <  2:377801/140‑1 (MQ=255)                                                                                                                                    gACGGCGCCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGatat  <  2:208320/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                      TGAGCGCAGTTACTGGCCAGATTAGCGCGAAAATTATCATCGTTATTAGTAGGTTAACTCTCATTCTCTTCTCCTGTCTGCGTCACCGCCGCCCACGATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATAT  >  NC_000913/306754‑307024

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