Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,908,514	A→G	intergenic (‑148/+35)	pstB ← / ← pstA	phosphate ABC transporter ATPase/phosphate ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,908,514	0	A	G	100.0%	107.3 / NA	32	intergenic (‑148/+35)	pstB/pstA	phosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/16);  total (16/16)

ATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAA  >  NC_000913/3908382‑3908623                                                                                                                                     |                                                                                                               aTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTa                                                                                                          >  2:325911/1‑139 (MQ=255)                         tCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                                  <  2:465584/139‑1 (MQ=255)                           aGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCaa                                                                                >  2:728990/1‑139 (MQ=255)                             gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCa                                                                                    >  2:641007/1‑133 (MQ=255)                             gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCa                                                                                    <  1:641007/133‑1 (MQ=255)                                  ccGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                         <  1:162248/139‑1 (MQ=255)                                  ccGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                         >  2:162248/1‑139 (MQ=255)                                  ccGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                         <  1:728990/139‑1 (MQ=255)                                  ccGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                         <  2:512899/139‑1 (MQ=255)                                       ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgtgt                                                                    >  2:271156/1‑139 (MQ=255)                                       ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgtgt                                                                    >  1:648369/1‑139 (MQ=255)                                            gATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCGGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAtt                                                               >  1:174752/1‑139 (MQ=255)                                                                   ggCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCa                                       <  1:271156/139‑1 (MQ=255)                                                                               aCCGTGTACAGAGTCAAAAAGGCGCTCGGCCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAGCCGGGCCGGCATTA‑CCACCCTGGTTTTTCTTCCGAAAAAACAACCCCCGCAAGAATGttaacca                           >  2:534668/1‑134 (MQ=16)                                                                                  tgtaccgaagTAAAAAGGGCGCTGGGCCGGCCCCTTCGCCCCTCTGGGG‑GGCGGTTAAGGGGAAGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAAAATGTTcagcagcag                       <  1:534668/131‑1 (MQ=255)                                                                                    gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCg                                                                                         <  2:353293/73‑1 (MQ=255)                                                                                    gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCg                                                                                         >  1:353293/1‑73 (MQ=255)                                                                                    gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCaa                                                                          <  1:244680/88‑1 (MQ=255)                                                                                    gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCaa                                                                          >  2:244680/1‑88 (MQ=255)                                                                                    gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGt                      >  1:156004/1‑139 (MQ=255)                                                                                    gtACAGAGTCAAAAAAGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGt                      >  2:741745/1‑139 (MQ=255)                                                                                         gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                 <  1:630220/139‑1 (MQ=255)                                                                                         gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                 <  2:641854/139‑1 (MQ=255)                                                                                         gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                 <  2:6630/139‑1 (MQ=255)                                                                                          agTCCAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAACAGTACGcac                >  1:151824/1‑139 (MQ=255)                                                                                          agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcac                >  2:164368/1‑139 (MQ=255)                                                                                          agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcac                >  1:647031/1‑139 (MQ=255)                                                                                           ttCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcaca               <  1:325911/138‑1 (MQ=255)                                                                                           gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcaca               <  1:164368/139‑1 (MQ=255)                                                                                                        tcggtcggTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAACAGTACGCACAGGGTAATGATCaa  <  2:151824/139‑1 (MQ=255)                                                                                                                            tctGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                 <  2:414467/104‑1 (MQ=255)                                                                                                                            tctGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                 >  1:414467/1‑104 (MQ=255)                                                                                                                                     |                                                                                                               ATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTT‑CGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAA  >  NC_000913/3908382‑3908623

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