Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,643,992	T→C	T105A (ACA→GCA)	ydfU ←	Qin prophage; protein YdfU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,643,992	0	T	C	100.0%	72.1 / NA	21	T105A (ACA→GCA)	ydfU	Qin prophage; protein YdfU
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/11);  total (10/11)

GTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAAT  >  NC_000913/1643864‑1644119                                                                                                                                 |                                                                                                                                gTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGCTTCa                                                                                                                             <  1:384626/133‑1 (MQ=255) gTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGCTTCa                                                                                                                             <  2:469716/133‑1 (MQ=255)     aTCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGCTTCAgtgt                                                                                                                         >  2:477569/1‑133 (MQ=255)        aCCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGCTTCAGTGTGAt                                                                                                                      <  1:477569/133‑1 (MQ=255)             tGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGCTTCAGTGTGATGGTaa                                                                                                                 <  2:394729/133‑1 (MQ=255)                 gTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGCTTCAGTGTGATGGTAATCGc                                                                                                             <  2:523883/133‑1 (MQ=255)                           gttgAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGCTTCAGTGTGATGGTAATCGCCATTAGGCCa                                                                                                   <  1:437424/133‑1 (MQ=255)                                         tCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGCTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTaa                                                                                     >  1:76730/1‑133 (MQ=255)                                           cAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGCTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACa                                                                                   <  1:585585/133‑1 (MQ=255)                                             gTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGCTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATg                                                                                 <  2:579627/133‑1 (MQ=255)                                             gTGATTCCGACGTCTGGTCACGCATTTGTTTGTCGCAGTGCTAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGCTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATg                                                                                 <  2:214585/133‑1 (MQ=255)                                              tGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGCTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGa                                                                                >  2:566973/1‑133 (MQ=255)                                                                                                        ccGGTACCCTAACGGGGAGTGACTGCTTCATTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCcgttcgt                      <  1:193057/133‑1 (MQ=255)                                                                                                        ccGGTACCATAACGGTGAATGACTGCTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCcgttcgt                      >  1:46770/1‑133 (MQ=255)                                                                                                           gTACCATAACGGTGAATGACTGCTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTAc                   >  2:326879/1‑133 (MQ=255)                                                                                                                      gTGAATGACTGCTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCTCGAATCAGCCGTTCTTTACTATAAAACGgc        >  1:172072/1‑133 (MQ=255)                                                                                                                         aaTGACTGCTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGGCAGCAGCAAGAAGCACGCGTTCGTTCCTAAAAAACGgcagc     >  2:199788/1‑133 (MQ=255)                                                                                                                         aaTGACTGCTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGgcagc     >  2:257144/1‑133 (MQ=255)                                                                                                                         aaTGACTGCTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGgcagc     >  1:421302/1‑133 (MQ=255)                                                                                                                         aaTGACTGCTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGgcagc     >  2:600328/1‑133 (MQ=255)                                                                                                                         aaTGACTGCTGCAGTATGATGGTAATCGGCATGAGGCCACTCGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGgcagc     <  2:257595/133‑1 (MQ=255)                                                                                                                            gACTGCTTCACTGTGATGGTAATCCCCACTAGACCACTGGCTGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAAt  <  2:172072/133‑1 (MQ=255)                                                                                                                                 |                                                                                                                                GTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTGGTCACGCAGTTGGTTGTCGCAGTGCCAGCACAACACCATTGCGCCGGTACCATAACGGTGAATGACTGTTTCAGTGTGATGGTAATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAAT  >  NC_000913/1643864‑1644119

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