Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,632,837	A→G	intergenic (‑552/‑53)	ydfJ ← / → ynfT	putative transporter YdfJ/Qin prophage; protein YnfT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,632,837	0	A	G	100.0%	83.5 / NA	36	intergenic (‑552/‑53)	ydfJ/ynfT	putative transporter YdfJ/Qin prophage; protein YnfT
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/19);  total (17/19)

AGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACAT  >  NC_000913/1632703‑1632872                                                                                                                                       |                                    agagTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATaa                                <  1:50405/140‑1 (MQ=18)      ctctACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAAtt                           >  2:159486/1‑140 (MQ=37)      ctctACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAAtt                           <  1:676217/140‑1 (MQ=18)       tctACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTg                          >  2:114960/1‑140 (MQ=21)       tctACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTg                          >  1:751402/1‑140 (MQ=25)        ctACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGt                         >  2:852322/1‑140 (MQ=25)        ctACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGt                         >  2:164329/1‑140 (MQ=25)        ctACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGt                         >  2:10015/1‑140 (MQ=25)        ctACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGt                         >  1:618205/1‑140 (MQ=25)        ctACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGt                         <  1:613064/140‑1 (MQ=18)        ctACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATA‑CTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTc                        >  2:156889/1‑140 (MQ=255)         tACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTc                        <  2:914888/140‑1 (MQ=18)         tACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTc                        >  1:15639/1‑140 (MQ=21)         tACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTc                        <  1:159486/140‑1 (MQ=18)           cATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCaa                      <  1:526276/140‑1 (MQ=14)             tAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACt                    <  1:578406/140‑1 (MQ=14)             tAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACt                    <  1:609977/140‑1 (MQ=14)             tAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACt                    <  1:678622/140‑1 (MQ=17)              aaGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACt                    <  2:838578/139‑1 (MQ=14)              aaGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACt                    >  1:838578/1‑139 (MQ=18)                gATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACtata                 <  1:10015/140‑1 (MQ=14)                gATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACtata                 >  2:529067/1‑140 (MQ=18)                gATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACtata                 <  1:212516/140‑1 (MQ=14)                gATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACtata                 >  1:727837/1‑140 (MQ=18)                  tAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATc               <  2:194481/140‑1 (MQ=14)                  tAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATc               <  1:638457/140‑1 (MQ=14)                  tAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATc               <  1:513131/140‑1 (MQ=14)                   aGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCa              >  1:556454/1‑140 (MQ=18)                    gTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCat             <  1:164308/140‑1 (MQ=14)                     tttGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACtata                 <  1:210441/135‑1 (MQ=14)                     tttGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACtata                 >  2:210441/1‑135 (MQ=18)                            ttGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAaca     <  1:452811/140‑1 (MQ=14)                            ttGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAaca     <  2:906108/140‑1 (MQ=17)                              gCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAAcaca   >  1:149707/1‑140 (MQ=17)                               ccGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACAt  >  1:452515/1‑140 (MQ=18)                               ccGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACAt  >  2:47157/1‑140 (MQ=18)                                                                                                                                       |                                    AGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACAT  >  NC_000913/1632703‑1632872

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