Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,634,948	A→G	A300A (GCT→GCC)	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,634,948	0	A	G	100.0%	126.1 / NA	36	A300A (GCT→GCC)	ydfN	Qin prophage; putative side tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/13);  total (23/13)

TGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCT  >  NC_000913/1634888‑1635080                                                             |                                                                                                                                     tGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCAt                                                       >  2:261930/1‑140 (MQ=18)      gCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATgcgcg                                                  >  1:115989/1‑140 (MQ=38)        cTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGcc                                                >  1:863169/1‑140 (MQ=255)           acaATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTgc                                             >  2:265579/1‑140 (MQ=255)           acaATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTgc                                             >  1:25057/1‑140 (MQ=255)           acaATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTgc                                             <  2:25057/140‑1 (MQ=255)           acaATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTgc                                             >  2:571534/1‑140 (MQ=255)              atatAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCt                                          <  1:265579/140‑1 (MQ=255)                  aGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgca                                      >  2:878538/1‑140 (MQ=255)                    ttAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgcagc                                    >  2:337445/1‑140 (MQ=255)                    ttAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgcag                                     >  1:364282/1‑139 (MQ=255)                    ttAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgcag                                     >  1:541539/1‑139 (MQ=255)                    ttAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgcag                                     >  2:344604/1‑139 (MQ=255)                    ttAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgcag                                     <  1:344604/139‑1 (MQ=255)                    ttAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgcag                                     <  2:541539/139‑1 (MQ=255)                    ttAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTgcag                                     <  2:364282/139‑1 (MQ=255)                       aaTGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTg                                        <  2:71597/133‑1 (MQ=255)                       aaTGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTg                                        >  1:71597/1‑133 (MQ=255)                       aaTGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGt                                 >  1:26292/1‑140 (MQ=255)                       aaTGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGt                                 >  1:1160300/1‑140 (MQ=255)                       aaTGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGt                                 >  1:1017372/1‑140 (MQ=255)                        aTGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTg                                <  2:1053222/140‑1 (MQ=255)                           caaTGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCg                             >  2:35748/3‑140 (MQ=255)                           cGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCg                             <  2:26292/140‑1 (MQ=255)                           cGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCg                             >  2:948857/1‑140 (MQ=255)                             aTGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGa                           >  2:498102/1‑140 (MQ=255)                              tGTTTTTGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGAc                          >  1:922887/1‑140 (MQ=255)                                   ttGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATTTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGa                     >  2:954005/1‑140 (MQ=255)                                    tGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAg                    <  1:878538/140‑1 (MQ=255)                                    tGACGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAg                    <  1:35748/140‑1 (MQ=255)                                      aCGGTGTTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAGTg                  <  1:919666/140‑1 (MQ=255)                                          tgtTTTCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTa              >  2:858910/1‑140 (MQ=255)                                               tCCGCGTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGc         <  1:571534/140‑1 (MQ=255)                                                    gTTACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCAt    >  1:1158163/1‑140 (MQ=255)                                                     ttACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATc   <  2:922887/140‑1 (MQ=255)                                                      tACCAGCGGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCt  >  2:1097406/1‑140 (MQ=255)                                                             |                                                                                                                                     TGCAAGCCTCACAATATAGTTAAATGCGATGTTTTTGACGGTGTTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCT  >  NC_000913/1634888‑1635080

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