Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,587,323	C→T	A138T (GCC→ACC)	ydeQ ←	putative fimbrial adhesin protein YdeQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,587,323	0	C	T	100.0%	83.7 / NA	25	A138T (GCC→ACC)	ydeQ	putative fimbrial adhesin protein YdeQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/13);  total (12/13)

AAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGCTGAGACAAGTC  >  NC_000913/1587077‑1587560                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              agaaTGAGAGCTTTTGTGCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGGCCGGGAGGTCCACGGTGACATTTCGACTATCAACGGGGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTAGACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCTCCGgccgcc                                                                                                                                                                                                                                           <  2:413209/249‑1 (MQ=255) aaaaTGAGAGCTTTTGTTCGCTTGAGCAATATCAGCCAAGGGGCATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGgtcgcc                                                                                                                                                                                                                                           >  2:303061/1‑251 (MQ=255)      gagaGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTACACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAAc                                                                                                                                                                                                                                      >  1:266735/1‑251 (MQ=255)         agCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAATGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTgg                                                                                                                                                                                                                                   <  2:300389/251‑1 (MQ=255)              ttGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGCGTa                                                                                                                                                                                                                               >  2:256631/1‑250 (MQ=255)                    gCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTGTATGACAACTCCACCGGTCGCCCCAACTGTGGTAATATaga                                                                                                                                                                                                                        <  2:236053/251‑1 (MQ=255)                            aaTATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATCGAGTTTTAAc                                                                                                                                                                                                                >  1:283587/1‑251 (MQ=255)                                   gCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATg                                                                                                                                                                                                         <  1:130515/251‑1 (MQ=255)                                                      cTTCCAGGACAGTCCGGGAGATCCAAGGTGACATTTCGACTATCAACGGTGCAGCCGCATGTTGGAATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGACTACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCCTCGGTGTCTTATCAcc                                                                                                                                                                                      >  2:536992/1‑251 (MQ=255)                                                          cAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAat                                                                                                                                                                                   >  1:14515/1‑250 (MQ=255)                                                             gAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATcc                                                                                                                                                                               <  1:256631/251‑1 (MQ=255)                                                              aaaGTCCGGGAGTTCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCa                                                                                                                                                                              <  2:264128/251‑1 (MQ=255)                                                               aaGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCaa                                                                                                                                                                             <  1:133401/251‑1 (MQ=255)                                                               aaGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCaa                                                                                                                                                                             <  1:276558/251‑1 (MQ=255)                                                                      ggAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACTGGCAATGACTTCCCCCGTTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTg                                                                                                                                                                      <  1:465700/251‑1 (MQ=255)                                                                        agaTCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAACTAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCATTGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGt                                                                                                                                                                    >  2:338143/1‑251 (MQ=255)                                                                                                                  gTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGt                                                                                                                           <  1:189300/250‑1 (MQ=255)                                                                                                                                      aTTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACt                                                                                                      >  1:498467/1‑251 (MQ=255)                                                                                                                                            aataatATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAggg                                                                                                >  1:577028/1‑251 (MQ=255)                                                                                                                                             ataataTTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAggg                                                                                                >  1:305323/1‑250 (MQ=255)                                                                                                                                                                                tCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTatat                                                            >  1:197913/1‑251 (MQ=255)                                                                                                                                                                                                          gAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAgggg                                                                                                                                              >  1:463948/1‑143 (MQ=255)                                                                                                                                                                                                          gAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAgggg                                                                                                                                              <  2:463948/143‑1 (MQ=255)                                                                                                                                                                                                                 ggCAATGACTTCGCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCGTCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGTAGGGGGTAGGTCACATTATTCCAGTACAGCTATCGTTTATATTACTGTAGGGAGGCTGCGAAAGCTGAGCGTTTTACGAGGTTTATATGAGAAGTGTCGTACCAGCCGCCGTAATCATTcc                           <  2:235166/251‑1 (MQ=255)                                                                                                                                                                                                                                     ttATGACAACTACACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGTAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGCTgaga        <  1:179780/250‑1 (MQ=255)                                                                                                                                                                                                                                       aTGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCAACAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGCTGAGACa      >  1:332901/1‑250 (MQ=255)                                                                                                                                                                                                                                          aCAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGTTAATGGCAGCGGTGTCTTAACACCAATATCCAATACGTTTGTTTTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCTATCCTTTATAGTACTGCAGGGAGCCTGCGAAAGCGGAGCCTGGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGCTGAGACAAGTc  <  1:584636/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              AAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGCTGAGACAAGTC  >  NC_000913/1587077‑1587560

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