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%	123.4 / NA	35	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/19);  total (16/19)

ATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACA  >  NC_000913/3908382‑3908610                                                                                                                                     |                                                                                                 aTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTa                                                                                            >  2:606494/1‑139 (MQ=255)   cGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAgg                                                                                          >  1:534844/1‑139 (MQ=255)     tGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGt                                                                                        <  1:404807/139‑1 (MQ=255)     tGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGt                                                                                        <  2:125820/139‑1 (MQ=255)       ccTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGa                                                                                      <  2:534844/139‑1 (MQ=255)       ccTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGa                                                                                      >  2:42431/1‑139 (MQ=255)           aTGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAgggg                                                                                  >  2:102731/1‑139 (MQ=255)                         tCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                    <  1:477119/139‑1 (MQ=255)                             gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGgggg                                                                                               <  2:475986/108‑1 (MQ=255)                             gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGgggg                                                                                               >  1:475986/1‑108 (MQ=255)                             gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAata                                                                >  2:25163/1‑139 (MQ=255)                                  ccGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                           <  1:562560/139‑1 (MQ=255)                                       ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgtgt                                                      >  1:618406/1‑139 (MQ=255)                                            gATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAtt                                                 >  1:468129/1‑139 (MQ=255)                                                                tCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgcg                             >  1:607158/1‑139 (MQ=255)                                                                     cAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                        >  2:354486/1‑139 (MQ=255)                                                                          tgtTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGtt                   <  1:42431/139‑1 (MQ=255)                                                                           gtTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                  <  1:102731/139‑1 (MQ=255)                                                                                    gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGt                                                                                        >  2:86666/1‑60 (MQ=255)                                                                                    gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGt                                                                                        <  1:86666/60‑1 (MQ=255)                                                                                        agagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGc     >  2:670627/1‑139 (MQ=255)                                                                                         gagTTAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca    <  2:648289/139‑1 (MQ=255)                                                                                         gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca    <  2:618406/139‑1 (MQ=255)                                                                                         gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca    >  2:395087/1‑139 (MQ=255)                                                                                         gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca    <  1:395087/139‑1 (MQ=255)                                                                                         gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca    <  1:25163/139‑1 (MQ=255)                                                                                         gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca    <  1:257376/139‑1 (MQ=255)                                                                                          agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCa                                                                      >  1:25060/1‑72 (MQ=255)                                                                                          agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCa                                                                      <  2:25060/72‑1 (MQ=255)                                                                                          agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                    >  2:424589/1‑74 (MQ=255)                                                                                          agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                    <  1:424589/74‑1 (MQ=255)                                                                                          agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcac   >  1:68112/1‑139 (MQ=255)                                                                                           gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcaca  <  1:354486/139‑1 (MQ=255)                                                                                           gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcaca  <  1:426557/139‑1 (MQ=255)                                                                                           gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcaca  <  1:606494/139‑1 (MQ=255)                                                                                           gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca    >  2:60724/1‑137 (MQ=255)                                                                                           gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca    <  1:60724/137‑1 (MQ=255)                                                                                                                                     |                                                                                                 ATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACA  >  NC_000913/3908382‑3908610

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