Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,907,703	T→C	T222A (ACG→GCG)	pstB ←	phosphate ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,907,703	0	T	C	100.0%	117.6 / NA	34	T222A (ACG→GCG)	pstB	phosphate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/17);  total (17/17)

TCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTT  >  NC_000913/3907574‑3907831                                                                                                                                  |                                                                                                                                 tCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGGCCATGTACATAAAGGCCGcgtg                                                                                                                               >  1:136337/1‑133 (MQ=255) tCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGcgtg                                                                                                                               >  2:509912/1‑133 (MQ=255)  ccATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTgg                                                                                                                              >  2:392066/1‑133 (MQ=255)           aCTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGaacaa                                                                                                                     >  2:480194/1‑133 (MQ=255)           aCTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCAAGGTACATAAAGGCAGCGTGGGCGGaacaa                                                                                                                     >  1:152427/1‑133 (MQ=255)           aCTCCTGAATCAACCGTAACGACCGGAGATGTAGACTTCTGTTTGTCTCTTCGCTGGCTTGGTGAACAGATCCACCGTGGTGCTGCACTCAATCAATTCGCTCAGGGACATAAACTCCGCGGGGGCGGaacaa                                                                                                                     >  2:510540/1‑133 (MQ=255)                       aCCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCctgc                                                                                                         >  2:500109/1‑133 (MQ=255)                       aCCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCctgc                                                                                                         >  1:351602/1‑133 (MQ=255)                       aCCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCctgc                                                                                                         >  2:256175/1‑133 (MQ=255)                       aCCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCctg                                                                                                          >  2:54725/1‑132 (MQ=255)                       aCCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCctg                                                                                                          <  1:54725/132‑1 (MQ=255)                         cGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCctgctg                                                                                                       >  2:417345/1‑133 (MQ=255)                          gTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCtgctgc                                                                                                      <  1:312689/133‑1 (MQ=255)                                    gTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTggg                                                                                            <  2:351602/133‑1 (MQ=255)                                               ttctGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATcacc                                                                                 <  1:305143/133‑1 (MQ=255)                                               ttctGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATcacc                                                                                 >  1:392013/1‑133 (MQ=255)                                                tttgtttgttCCTTCTCTGGCTTGTTCAGCAGCTCGTCCGTTTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATcacca                                                                                <  2:235660/131‑1 (MQ=255)                                                 ctGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGATCAACGCGCTGCCCGCTACATGTTGTCGGTGTCGATcaccac                                                                               >  2:154549/1‑133 (MQ=255)                                                      tgtttCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGgtgt                                                                          >  2:108609/1‑133 (MQ=255)                                                       gtttCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTa                                                                         >  2:382972/1‑133 (MQ=255)                                                                cTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCtt                                                                <  2:276242/133‑1 (MQ=255)                                                                cTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCtt                                                                <  1:392066/133‑1 (MQ=255)                                                                       ggTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCg                                                         <  1:256175/133‑1 (MQ=255)                                                                             cAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGtt                                                           <  1:457882/125‑1 (MQ=255)                                                                             cAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGtt                                                           >  2:457882/1‑125 (MQ=255)                                                                                  cGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGctc                                              >  1:391900/1‑133 (MQ=255)                                                                                     ccGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTc                                           >  2:296798/1‑133 (MQ=255)                                                                                            gCTGAACTCAGTCAAATCGCCCAGCTACCTAATCGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGc                                    <  2:152427/133‑1 (MQ=255)                                                                                            gCTGAACTCACTCAATTCGCCCAGCTACATAACCGCCGCGTCGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGc                                    <  2:136337/133‑1 (MQ=255)                                                                                             cTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGcc                                   <  1:296798/133‑1 (MQ=255)                                                                                                 aCTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGt                               <  2:392013/133‑1 (MQ=255)                                                                                                 aCTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGt                               <  1:509912/133‑1 (MQ=255)                                                                                                 aCTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGt                               <  1:108609/133‑1 (MQ=255)                                                                                                 aCTCAATCAATTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTAACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGt                               <  1:500109/133‑1 (MQ=255)                                                                                                   tcaatcaaTTCGCCCAGGTACATAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTag                             <  1:480194/133‑1 (MQ=255)                                                                                                                        aTAAACGCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGCTAGGGTCGACCGCCGAac        >  1:236705/1‑133 (MQ=255)                                                                                                                              gCCGCGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGtt  >  2:50675/1‑133 (MQ=255)                                                                                                                                  |                                                                                                                                 TCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTT  >  NC_000913/3907574‑3907831

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