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%	168.2 / NA	47	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 (23/24);  total (23/24)

AATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCA  >  NC_000913/3908404‑3908622                                                                                                               |                                                                                                             aaTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCa                                                                                  <  1:344807/139‑1 (MQ=255)     aGGCTTACCGTAGGCCGGTTCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCaa                                                                              >  2:464791/1‑139 (MQ=255)                 ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTGCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgtgt                                                                  >  2:290503/1‑139 (MQ=255)                 ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCa                                                                                  <  1:778421/123‑1 (MQ=255)                 ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCa                                                                                  >  2:778421/1‑123 (MQ=255)                 ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAatat                                                                           >  1:5431/1‑130 (MQ=255)                 ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAatat                                                                           <  2:5431/130‑1 (MQ=255)                 ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgtgt                                                                  >  2:208719/1‑139 (MQ=255)                   ccGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                       >  2:230878/1‑132 (MQ=255)                   ccGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                       <  1:230878/132‑1 (MQ=255)                      gATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAtt                                                             >  2:491972/1‑139 (MQ=255)                      gATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAtt                                                             >  2:289562/1‑139 (MQ=255)                      gATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAtt                                                             >  2:298796/1‑139 (MQ=255)                            gCGTTCACGCCACATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                       <  2:425773/123‑1 (MQ=255)                            gCGTTCACGCCACATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                       >  1:425773/1‑123 (MQ=255)                                        cATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcg                                           <  1:50483/139‑1 (MQ=255)                                          tCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgcg                                         <  2:302471/139‑1 (MQ=255)                                              gCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCaa                                                                              <  1:268606/98‑1 (MQ=255)                                              gCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCaa                                                                              >  2:268606/1‑98 (MQ=255)                                               cAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                                    <  1:491972/139‑1 (MQ=255)                                                  aGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCgt                                                                                      <  1:528356/86‑1 (MQ=255)                                                  aGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCgt                                                                                      >  2:528356/1‑86 (MQ=255)                                                         aCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                          >  1:625809/1‑139 (MQ=255)                                                         aCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                          >  2:673233/1‑139 (MQ=255)                                                         aCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                          <  1:298796/139‑1 (MQ=255)                                                          ccGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcag                         <  2:625809/139‑1 (MQ=255)                                                              gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGt                     >  1:412118/1‑139 (MQ=255)                                                                   gagTCAAAAAGGCGCTCGGTCGGTGCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                <  1:290503/139‑1 (MQ=255)                                                                   gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                <  1:673233/139‑1 (MQ=255)                                                                    agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                       >  1:482750/1‑83 (MQ=255)                                                                    agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                       <  2:482750/83‑1 (MQ=255)                                                                    agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                <  1:427769/138‑1 (MQ=255)                                                                    agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGAACGca                >  2:427769/1‑138 (MQ=255)                                                                     gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                          <  1:711279/127‑1 (MQ=255)                                                                     gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                          >  2:711279/1‑127 (MQ=255)                                                                       cAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGtt                                                                                                         <  1:752714/46‑1 (MQ=255)                                                                       cAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGtt                                                                                                         >  2:752714/1‑46 (MQ=255)                                                                        aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                <  1:470073/134‑1 (MQ=255)                                                                        aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                >  2:470073/1‑134 (MQ=255)                                                                          aaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTa         <  2:412118/139‑1 (MQ=255)                                                                          aaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTa         <  1:289562/139‑1 (MQ=255)                                                                           aaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTaa        <  1:208719/139‑1 (MQ=255)                                                                                gcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                >  2:206695/1‑126 (MQ=255)                                                                                gcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                <  1:206695/126‑1 (MQ=255)                                                                                 cTCGGTCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCa  >  1:528180/1‑139 (MQ=255)                                                                                        ggtcCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCgt                                                                                      >  2:473503/1‑48 (MQ=255)                                                                                        ggtcCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAACCgt                                                                                      <  1:473503/48‑1 (MQ=255)                                                                                                            gaggGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                          >  2:298288/1‑88 (MQ=255)                                                                                                            gaggGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                          <  1:298288/88‑1 (MQ=255)                                                                                                               |                                                                                                             AATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCA  >  NC_000913/3908404‑3908622

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