Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,450,781	T→C	T99A (ACT→GCT)	yfcQ ←	putative fimbrial protein YfcQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,450,781	0	T	C	100.0%	68.9 / NA	21	T99A (ACT→GCT)	yfcQ	putative fimbrial protein YfcQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/5);  total (16/5)

GATAGCAAGCGGAGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGT  >  NC_000913/2450707‑2450856                                                                           |                                                                            gATAGCAAGCGGAGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGCTTCAAt                                                                       >  1:387848/1‑81 (MQ=255)     gCAAGCGGAGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGCTTCAATtgct                                                                   <  2:387837/81‑1 (MQ=255)        aGCGGAGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGCTTCAATTGCTGca                                                                >  2:387861/1‑81 (MQ=255)         gCGGAGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGCTTCAATTGCTGcat                                                               >  1:387870/1‑81 (MQ=255)             aGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGCTTCAATTGCTGcatcat                                                            >  2:387856/1‑80 (MQ=255)                 tttaGTTTAAAACGCTGACCATCATCCTGAAGCTCGATACCAAAACTGGGTACCTCAGCTTCTATTGCTGCATCATTaaac                                                       >  2:895309/4‑80 (MQ=255)                      tttAAAACGCTGACCTACGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGCTTCAATTGCTGCATCATTAAAACTGGt                                                  <  2:387838/81‑1 (MQ=255)                           aaCGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGCCTCAATTGCTGCATCATTAAAACTGGTTtgtg                                             >  2:387858/1‑81 (MQ=255)                              gCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGCTTCAATTGCTGCATCATTAAAACTGGTTTGTGTAc                                          >  2:387857/1‑81 (MQ=255)                                tGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGCTTCAATTGCTGCATCATTAAAACTGGTTTGTGTAc                                          >  1:387854/1‑79 (MQ=255)                                   ccATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGCTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTc                                     <  2:387839/81‑1 (MQ=255)                                   ccATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGCTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTc                                     >  2:387866/1‑81 (MQ=255)                                     aTCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGCTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTc                                     >  1:387862/1‑79 (MQ=255)                                        gTGCTGAAGCTCGATACCAAAACCGGGTACATCAGCTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTCCAgg                                 <  1:387840/80‑1 (MQ=255)                                        gTGCTGAAGCTCGATACCAAAACCGGGTACATCAGCTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTCCAGGt                                >  1:387869/1‑81 (MQ=255)                                                    gATACCAAAACCGGGTACATCAGCTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTCCAGGTGAGGGTCAtc                      <  1:387841/79‑1 (MQ=255)                                                       aCCAAAACCGGGTACATCAGCTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTCCAGGTGAGGGTCATCTccc                  >  2:387867/1‑80 (MQ=255)                                                           aaaCCGGGTACATCAGCTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTCCAGGTGAGGGTCATCTCCCAATCt             >  1:387859/1‑81 (MQ=255)                                                                    aCATCAGCTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACgg    >  1:387871/1‑81 (MQ=255)                                                                    aCATCAGCTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACgg    >  2:387868/1‑81 (MQ=255)                                                                      aTCAGCTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGgt  >  1:387872/1‑81 (MQ=255)                                                                           |                                                                            GATAGCAAGCGGAGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGT  >  NC_000913/2450707‑2450856

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