Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,635,029	C→T	A273A (GCG→GCA)	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,029	0	C	T	100.0%	95.6 / NA	30	A273A (GCG→GCA)	ydfN	Qin prophage; putative side tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/17);  total (13/17)

TTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCG  >  NC_000913/1634952‑1635159                                                                              |                                                                                                                                   ttAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCtt                                                                      >  2:615360/1‑140 (MQ=255) ttAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCtt                                                                      >  1:337135/1‑140 (MQ=255)  tAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTc                                                                     <  2:524619/140‑1 (MQ=255)    aCGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTccc                                                                   <  2:455287/140‑1 (MQ=255)           tGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTa                                                            <  1:44630/140‑1 (MQ=255)                tgtCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                       <  2:381481/140‑1 (MQ=255)                tgtCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                       >  1:381481/1‑140 (MQ=255)                tgtCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                       <  1:601674/140‑1 (MQ=255)                  tCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCt                                                     >  2:290519/1‑140 (MQ=255)                         gAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACt                                              <  1:144656/140‑1 (MQ=255)                         gAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACt                                              <  2:127060/140‑1 (MQ=255)                               aTCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCt                                        >  2:461436/1‑140 (MQ=255)                                   cAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCttgtt                                    <  2:3711/140‑1 (MQ=255)                                    aaCGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCttgttg                                   <  2:337135/140‑1 (MQ=255)                                           tGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTtgtg                                 >  1:655357/1‑135 (MQ=21)                                           tGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTtgtg                                 <  2:655357/135‑1 (MQ=255)                                                  gAGCACCAATCCCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATcccc                     <  2:447187/140‑1 (MQ=25)                                                  gAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATcccc                     >  1:582635/1‑140 (MQ=21)                                                   aGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCcg                    >  1:600156/1‑140 (MQ=21)                                                   aGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCcg                    >  1:680841/1‑140 (MQ=21)                                                       ccAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGCGTGGGTATGAGCACCATCTGATGATGTCTTCCCTCCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCt                <  1:615360/140‑1 (MQ=25)                                                       ccAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCt                >  1:51584/1‑140 (MQ=25)                                                            aCCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCC‑GACAGCGAGTGGGTATTAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGAtt          >  2:720338/1‑140 (MQ=37)                                                               gACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTg        >  2:243386/1‑140 (MQ=37)                                                                aCAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGa       >  2:192383/1‑140 (MQ=255)                                                                   gTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAAt    <  1:461436/140‑1 (MQ=255)                                                                   gTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAAt    <  1:274080/140‑1 (MQ=255)                                                                    tATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCAGGCCACTACCGCTTGTTGTACTCATAATCCCCGCGCTTAGATTTGAAATc   <  1:164627/140‑1 (MQ=37)                                                                     aTGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATcg  <  1:290519/140‑1 (MQ=255)                                                                     aTGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATATCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATcg  <  2:560887/140‑1 (MQ=37)                                                                              |                                                                                                                                   TTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCG  >  NC_000913/1634952‑1635159

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