Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,908,805	Δ2 bp	coding (634‑635/891 nt)	pstA ←	phosphate ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,908,805	0	G	.	100.0%	232.2 / NA	48	coding (635/891 nt)	pstA	phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (22/26);  total (22/26)
*	NC_000913	3,908,806	0	T	.	100.0%	232.2 / NA	48	coding (634/891 nt)	pstA	phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (22/26);  total (22/26)

CGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGTCATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAGGCTGTACGGCACCAGTTTCAGCATGTTCTCGGTGGTGCGGATAACAATCGGCAC  >  NC_000913/3908669‑3908947                                                                                                                                         ||                                                                                                                                              cGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGaca                                                                                                                                 <  1:894880/150‑1 (MQ=255)    aaaCTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCg                                                                                                                              >  2:914087/1‑150 (MQ=255)     aaCTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCCGCAGGATACCG‑‑CATAATCCCGGACACCGAc                                                                                                                            <  2:606783/151‑1 (MQ=255)          aaaGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCttt                                                                                                                       <  1:940125/151‑1 (MQ=255)          aaaGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCttt                                                                                                                       >  2:180120/1‑151 (MQ=255)           aaGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTc                                                                                                                      <  1:335348/151‑1 (MQ=255)             gATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTc                                                                                                                      >  2:689536/1‑149 (MQ=255)             gATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTc                                                                                                                      >  1:353618/1‑149 (MQ=255)                   cACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAAt                                                                                                              >  2:453798/1‑151 (MQ=255)                    aCCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATc                                                                                                             <  2:24862/151‑1 (MQ=255)                         cAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCaga                                                                                                        <  1:180128/151‑1 (MQ=255)                          aGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCa                                                                                                          >  1:974642/1‑148 (MQ=255)                            gTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGat                                                                                                     <  2:630907/151‑1 (MQ=255)                             ttGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCag                                                                                                         >  1:72745/1‑146 (MQ=255)                                gCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGatcatc                                                                                                 >  2:883813/1‑151 (MQ=255)                                 cGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCt                                                                                                >  1:327311/1‑151 (MQ=255)                                  gATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCtt                                                                                               >  1:939296/1‑151 (MQ=255)                                  gATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCtt                                                                                               >  1:291813/1‑151 (MQ=255)                                    tCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTcc                                                                                             >  2:89021/1‑151 (MQ=255)                                    tCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTc                                                                                              <  2:923981/150‑1 (MQ=255)                                     cGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTcc                                                                                             <  1:826612/150‑1 (MQ=255)                                     cGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTcc                                                                                             >  1:25290/1‑150 (MQ=255)                                     cGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCa                                                                                            <  1:53017/151‑1 (MQ=255)                                       gCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACt                                                                                          >  2:669231/1‑151 (MQ=255)                                          gCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACtt                                                                                         <  2:72386/149‑1 (MQ=255)                                          gCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACtt                                                                                         <  2:1045185/149‑1 (MQ=255)                                          gCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCg                                                                                       >  1:551348/1‑151 (MQ=255)                                           catcatGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCgg                                                                                      <  2:415834/151‑1 (MQ=255)                                               atGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGgtgt                                                                                   >  1:741004/1‑150 (MQ=255)                                                   ccGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCa                                                                              <  2:664872/151‑1 (MQ=255)                                                    cGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAg                                                                             <  1:228843/151‑1 (MQ=255)                                                             gAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATaa                                                                      >  2:52342/1‑149 (MQ=255)                                                                 tGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCt                                                                >  1:199355/1‑151 (MQ=255)                                                                       gagaGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGc                                                          <  1:851479/151‑1 (MQ=255)                                                                         gagCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAg                                                        <  1:548538/151‑1 (MQ=255)                                                                         gagCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAg                                                        <  1:883647/151‑1 (MQ=255)                                                                            cgcgGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAGGCt                                                     <  1:631618/151‑1 (MQ=255)                                                                                  aaaCAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAGGCTGTACg                                                <  1:749187/150‑1 (MQ=255)                                                                                            ggcgCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAGGCTGTACGGCACCAGTTTc                                     <  1:1059243/151‑1 (MQ=255)                                                                                               gcgGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAGGCTGTACGGCACCAGTTTCAGc                                  <  1:6691/151‑1 (MQ=255)                                                                                                cgGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAGGCTGTACGGCACCAGTTTCAGc                                  >  2:728319/1‑150 (MQ=255)                                                                                                     tCACCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAGGCTGTACGGCACCAGTTTCAGCATGTtc                            >  2:736310/1‑151 (MQ=255)                                                                                                       aCCAGCAATACGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAGGCTGTACGGCACCAGTTTCAGCATGTtctc                          <  1:567960/151‑1 (MQ=255)                                                                                                                 cGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAGGCTGTACGGCACCAGTTTCAGCATGTTCTCGGTGGTGCgg                <  2:668051/151‑1 (MQ=255)                                                                                                                 cGGGCAATCGCCAGCAGGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAGGCTGTACGGCACCAGTTTCAGCATGTTCTCGGTGGTGCgg                >  1:952884/1‑151 (MQ=255)                                                                                                                      aaTCGCCAGCAGGAGACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAGGCTGTACGGCACCAGTTGCAGCATGTTCTCGGTGGTGCGGATaa            <  2:347585/150‑1 (MQ=255)                                                                                                                             agcagGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAGGCTGTACGGCACCAGTTTCAGCATGTTCTCGGTGGTGCGGATAACAATCGGc    <  1:296056/151‑1 (MQ=255)                                                                                                                               cagGATACCG‑‑CATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAGGCTGTACGGCACCAGTTTCAGCATGTTCTCGGTGGTGCGGATAACAATCGGCAc  >  1:740716/1‑151 (MQ=255)                                                                                                                                         ||                                                                                                                                              CGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGTCATAATCCCGGACACCGACGCTTTCAGCGTAATCGCAGAGATCATCTTCCACTTCGGTGTACCCAGCGCATAAGCCGCTTCACGCAGGCTGTACGGCACCAGTTTCAGCATGTTCTCGGTGGTGCGGATAACAATCGGCAC  >  NC_000913/3908669‑3908947

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