Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,910,466	T→C	intergenic (‑68/+19)	pstC ← / ← pstS	phosphate ABC transporter membrane subunit PstC/phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,910,466	0	T	C	100.0%	91.5 / NA	28	intergenic (‑68/+19)	pstC/pstS	phosphate ABC transporter membrane subunit PstC/phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/18);  total (10/18)

ACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAG  >  NC_000913/3910337‑3910580                                                                                                                                  |                                                                                                                   aCGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGTGGCGTAAAACACCGCAcc                                                                                                                 <  2:301732/133‑1 (MQ=255) aCGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCAcc                                                                                                                 <  1:448706/133‑1 (MQ=255)      aaaaaTAATGTC‑CCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCt                                                                                                           <  1:11301/132‑1 (MQ=255)       aaaaTTATGTCGCCCTTTCTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCt                                                                                                           <  1:225914/133‑1 (MQ=255)          aTTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGa                                                                                                        <  1:307364/133‑1 (MQ=255)                 cGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGt                                                                                                      <  1:13881/128‑1 (MQ=255)                 cGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGt                                                                                                      >  2:13881/1‑128 (MQ=255)                                gTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTa                                                                                  >  1:2301/1‑133 (MQ=255)                                gTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTa                                                                                  >  1:95166/1‑133 (MQ=255)                                      tAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTa                                                                            <  1:22706/133‑1 (MQ=255)                                           gCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTc                                                                       >  2:460560/1‑133 (MQ=255)                                               gCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGCCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTa                                                                   >  1:85698/1‑133 (MQ=255)                                                   ggTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAatat                                                               >  1:312620/1‑133 (MQ=255)                                                     ttGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTg                                                             <  2:493795/133‑1 (MQ=255)                                                     ttGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTg                                                             <  2:500488/133‑1 (MQ=255)                                                       gCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGt                                                           <  1:460560/133‑1 (MQ=255)                                                               aaGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACg                                                   <  1:289926/133‑1 (MQ=255)                                                               aaGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACg                                                   <  2:2301/133‑1 (MQ=255)                                                                          cAGTTAAACGTGTTTACAAAGTCATTTTTTAATGCCGGATGAGGCATAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAAcc                                        <  2:209153/133‑1 (MQ=255)                                                                          cAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAAcc                                        >  1:125686/1‑133 (MQ=255)                                                                          cAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATGGGTCTTACACGCATCGCGAAcc                                        >  1:209153/1‑133 (MQ=255)                                                                                aaCGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCa                                  <  1:5129/133‑1 (MQ=255)                                                                                 aCGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCaa                                 <  2:10328/133‑1 (MQ=255)                                                                                 aCGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCaa                                 <  1:404545/133‑1 (MQ=255)                                                                                  cgcgTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAAc                                >  1:87672/1‑133 (MQ=255)                                                                                                   tttGTAATGCCGGATGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGc               <  2:85698/133‑1 (MQ=255)                                                                                                       tAATGCCGGATGCGGGGTTAAAGACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGc           <  1:173171/133‑1 (MQ=255)                                                                                                                aTGCGGCGTAAAACACCGCACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAAGTACACTAGCCGGCAAGCTGGCGTAATCCAg  >  1:133130/1‑133 (MQ=255)                                                                                                                                  |                                                                                                                   ACGCTGAAAATTATGTCGCCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAG  >  NC_000913/3910337‑3910580

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