Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,635,059	C→T	S263S (TCG→TCA)	ydfN ←	Qin prophage; putative side tail fiber assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,635,059	0	C	T	100.0%	82.1 / NA	31	S263S (TCG→TCA)	ydfN	Qin prophage; putative side tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/17);  total (14/17)

TGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCG  >  NC_000913/1634967‑1635192                                                                                             |                                                                                                                                      tgtCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                                                        <  1:68669/140‑1 (MQ=14) tgtCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                                                        <  1:320609/140‑1 (MQ=14)   tCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                                                        <  2:415536/138‑1 (MQ=14)   tCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                                                        >  1:415536/1‑138 (MQ=14)   tCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCt                                                                                      >  2:138943/1‑140 (MQ=14)     cATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                                                        <  1:159175/136‑1 (MQ=14)     cATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                                                        >  2:159175/1‑136 (MQ=14)         tgAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGC‑GTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACt                                                                               >  2:113274/1‑140 (MQ=17)          gAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACt                                                                               <  1:138943/140‑1 (MQ=17)          gAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACt                                                                               <  2:394400/140‑1 (MQ=17)             ccAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTAcc                                                                            >  2:329887/1‑140 (MQ=17)                  cGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCttg                                                                       <  1:303599/140‑1 (MQ=18)                     aaCGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCttgttg                                                                    <  2:21441/140‑1 (MQ=18)                       cGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTtgtg                                                                  <  2:6090/140‑1 (MQ=18)                                 aTGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATcc                                                        >  2:400598/1‑140 (MQ=21)                                      cACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCttgtt                                                                     <  1:107557/122‑1 (MQ=18)                                      cACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCttgtt                                                                     >  2:107557/1‑122 (MQ=18)                                               cGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGAttt                                          >  1:25357/1‑140 (MQ=25)                                                   aGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGaaa                                      >  1:178304/1‑140 (MQ=37)                                                     tATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATc                                    <  2:27865/140‑1 (MQ=37)                                                             tgtgcaccTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGCTGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATcg                                   >  1:34007/7‑133 (MQ=255)                                                             tgtgcaccTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATcg                                   <  2:34007/127‑1 (MQ=255)                                                              gtgcaccTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATaa                           <  2:92458/135‑1 (MQ=255)                                                               tgcaccTGCGCTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAAc                          >  2:166178/5‑140 (MQ=255)                                                                      gcgcTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAAt                                     >  2:227730/1‑122 (MQ=255)                                                                      gcgcTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAAt                                     <  1:227730/122‑1 (MQ=255)                                                                      gcgcTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATcgcg                                 <  2:153741/126‑1 (MQ=255)                                                                      gcgcTTGCAGCAGTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATcgcg                                 >  1:153741/1‑126 (MQ=255)                                                                               gcagTGCCGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCt          <  1:329887/140‑1 (MQ=255)                                                                                     ccGGACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGtt    >  2:96418/1‑140 (MQ=255)                                                                                       ggACAGTGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTACTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCg  <  1:75024/140‑1 (MQ=255)                                                                                             |                                                                                                                                      TGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCGCGGTATAACCATTAGGGAAAATGCTCGTGTTCG  >  NC_000913/1634967‑1635192

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: ‑