Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,632,712	A→G	intergenic (‑427/‑178)	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,712	0	A	G	100.0%	50.6 / NA	22	intergenic (‑427/‑178)	ydfJ/ynfT	putative transporter YdfJ/Qin prophage; protein YnfT
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/6);  total (16/6)

CATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAAGATAAAAATT  >  NC_000913/1632599‑1632847                                                                                                                  |                                                                                                                                        cattcattCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCtt                                                                                                               <  2:340620/140‑1 (MQ=17)        tcctccGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTg                                                                                                        >  1:195168/1‑140 (MQ=18)               cATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATa                                                                                                 >  2:293105/1‑140 (MQ=18)                          tataGGCTATATCACCTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGAt                                                                                      >  1:124106/1‑140 (MQ=17)                                   tatCACCTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTtgctgc                                                                             >  2:54461/1‑140 (MQ=25)                                   tatCACCTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTtgctgc                                                                             >  2:191078/1‑140 (MQ=21)                                   tatCACCTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTtgctgc                                                                             >  2:176229/1‑140 (MQ=21)                                       aCCTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCAcc                                                                         <  2:295053/140‑1 (MQ=25)                                        ccTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCt                                                                        <  1:330835/140‑1 (MQ=25)                                         cTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTg                                                                       >  1:315948/1‑140 (MQ=21)                                          ttGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGt                                                                      >  2:9088/1‑140 (MQ=25)                                          ttGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGGCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGt                                                                      >  2:322695/1‑140 (MQ=25)                                                                aaaaaCTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGaaa                                                <  1:293105/140‑1 (MQ=18)                                                                                gTATACTTTTAAACCCAATACGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTa                                >  2:177508/1‑140 (MQ=12)                                                                                  ataCTTTTAAACCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACg                              <  2:195168/140‑1 (MQ=17)                                                                                  ataCTTTTAAACCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACg                              >  1:186826/1‑140 (MQ=17)                                                                                            aCCCAATCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGa                    <  1:54461/140‑1 (MQ=17)                                                                                                aaTCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAacac                >  1:277660/1‑140 (MQ=17)                                                                                                  tCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACAccc              >  2:286486/1‑140 (MQ=17)                                                                                                  tCCGTGTAGAGTCTCTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACAccc              >  1:192071/1‑140 (MQ=17)                                                                                                             tctcTGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAAGATAAAAAt   >  1:102264/1‑140 (MQ=21)                                                                                                              ctctGCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAAGATAAAAAtt  >  2:308618/1‑140 (MQ=21)                                                                                                                  |                                                                                                                                        CATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAAGATAAAAATT  >  NC_000913/1632599‑1632847

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