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%	205.9 / NA	57	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 (26/31);  total (26/31)

TCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGC  >  NC_000913/3908380‑3908638                                                                                                                                       |                                                                                                                             tCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGt                                                                                                                          <  1:1148621/139‑1 (MQ=255)  cATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGtt                                                                                                                         <  2:673654/139‑1 (MQ=255)   aTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTa                                                                                                                        >  2:957626/1‑139 (MQ=255)                  cTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAAc                                                                                                         >  2:355873/1‑139 (MQ=255)                   ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAAcc                                                                                                        <  2:398666/139‑1 (MQ=255)                   ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAAcc                                                                                                        >  2:154252/1‑139 (MQ=255)                   ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAAcc                                                                                                        >  1:543340/1‑139 (MQ=255)                          aTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAg                                                                                                 >  2:658857/1‑139 (MQ=255)                           tCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                                                >  2:604826/1‑139 (MQ=255)                           tCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                                                <  1:355873/139‑1 (MQ=255)                             aGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCaa                                                                                              <  2:1060913/139‑1 (MQ=255)                               gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAAc                                                                                                         >  1:1053472/1‑126 (MQ=255)                               gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAAc                                                                                                         <  2:1053472/126‑1 (MQ=255)                                    ccGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACAGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                                                <  1:1169562/130‑1 (MQ=255)                                    ccGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACAGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                                                >  2:1169562/1‑130 (MQ=255)                                         ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                                                >  1:1077744/1‑125 (MQ=255)                                         ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                                                <  2:1077744/125‑1 (MQ=255)                                              gATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAtt                                                                             >  2:53346/1‑139 (MQ=255)                                               aTCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttc                                                                            >  1:248134/1‑139 (MQ=255)                                                                               tAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcag                                            <  1:604826/139‑1 (MQ=255)                                                                                 aCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                                          <  1:129226/139‑1 (MQ=255)                                                                                      gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGt                                     <  2:964330/139‑1 (MQ=255)                                                                                      gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGt                                     >  2:611155/1‑139 (MQ=255)                                                                                      gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGt                                     >  1:1151548/1‑139 (MQ=255)                                                                                      gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGt                                     <  2:950326/139‑1 (MQ=255)                                                                                           gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                <  2:889135/139‑1 (MQ=255)                                                                                           gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                <  2:1151548/139‑1 (MQ=255)                                                                                           gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                <  1:658857/139‑1 (MQ=255)                                                                                           gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                <  1:957626/139‑1 (MQ=255)                                                                                           gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                <  1:53346/139‑1 (MQ=255)                                                                                            agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                                       >  1:1145182/1‑83 (MQ=255)                                                                                            agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                                       <  2:1145182/83‑1 (MQ=255)                                                                                            agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCa                                                      >  2:788207/1‑116 (MQ=255)                                                                                            agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCa                                                      <  1:788207/116‑1 (MQ=255)                                                                                            agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcac                               >  2:382932/1‑139 (MQ=255)                                                                                            agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcac                               >  2:836014/1‑139 (MQ=255)                                                                                            agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcac                               >  1:875507/1‑139 (MQ=255)                                                                                            agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCT‑GGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcaca                              <  2:936357/139‑1 (MQ=255)                                                                                             gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcg                                                           <  1:645560/110‑1 (MQ=255)                                                                                             gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcg                                                           >  2:645560/1‑110 (MQ=255)                                                                                             gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGt                                     <  2:231868/132‑1 (MQ=255)                                                                                             gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGt                                     >  1:231868/1‑132 (MQ=255)                                                                                             gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcaca                              <  2:528182/139‑1 (MQ=255)                                                                                             gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                >  1:7156/1‑137 (MQ=255)                                                                                             gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                >  1:666392/1‑137 (MQ=255)                                                                                             gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                <  2:666392/137‑1 (MQ=255)                                                                                             gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                <  2:7156/137‑1 (MQ=255)                                                                                                  aaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTa                         <  1:966034/139‑1 (MQ=255)                                                                                                      gcgcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGa                     >  2:1000261/1‑139 (MQ=255)                                                                                                               cggtcCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAccc            <  1:1000261/139‑1 (MQ=255)                                                                                                                        cGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgcg                                                         >  2:248193/1‑85 (MQ=255)                                                                                                                        cGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgcg                                                         <  1:248193/85‑1 (MQ=255)                                                                                                                        cGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCaa                 <  2:866828/125‑1 (MQ=255)                                                                                                                        cGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCaa                 >  1:866828/1‑125 (MQ=255)                                                                                                                         gCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCATTACGCACAGGGTAATGATCAATACCCCGGCCCAGGc  <  1:611155/139‑1 (MQ=255)                                                                                                                                tGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                <  2:592101/102‑1 (MQ=255)                                                                                                                                tGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                >  1:592101/1‑102 (MQ=255)                                                                                                                                       |                                                                                                                             TCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGC  >  NC_000913/3908380‑3908638

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