Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,907,945	T→C	N141S (AAC→AGC)	pstB ←	phosphate ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,907,945	0	T	C	100.0%	133.6 / NA	36	N141S (AAC→AGC)	pstB	phosphate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/17);  total (19/17)

CGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCAC  >  NC_000913/3907799‑3908090                                                                                                                                                   |                                                                                                                                                  cGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTcc                                                                                                                                                >  2:425760/1‑150 (MQ=255)     agagATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAAt                                                                                                                                           >  1:538439/1‑151 (MQ=255)      gagaTAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATg                                                                                                                                          <  1:523912/151‑1 (MQ=255)                    cgcCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAAt                                                                                                                            <  1:283004/151‑1 (MQ=255)                        gAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGccc                                                                                                                        >  1:384499/1‑151 (MQ=255)                         aaCACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGccc                                                                                                                        >  1:71132/1‑150 (MQ=255)                            acGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTg                                                                                                                    >  1:500615/1‑151 (MQ=255)                                        gcagcagcaCTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTc                                                                                                        >  1:430643/1‑151 (MQ=255)                                         cagcagcaCTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTcc                                                                                                       <  1:290816/151‑1 (MQ=255)                                           gcagcaCTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCa                                                                                                      <  2:525633/150‑1 (MQ=255)                                                   tCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCa                                                                                              <  1:598437/150‑1 (MQ=255)                                                       ggCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACggg                                                                                          >  2:438631/1‑150 (MQ=255)                                                       ggCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACggg                                                                                          >  2:442951/1‑150 (MQ=255)                                                         cGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACggg                                                                                          >  2:126783/1‑148 (MQ=255)                                                          gAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGgag                                                                                        >  1:55693/1‑149 (MQ=255)                                                                        cgcgCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGAc                                                                        <  1:231517/151‑1 (MQ=255)                                                                            cAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAAc                                                                    >  2:198323/1‑151 (MQ=255)                                                                               tACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGcc                                                                 >  2:587211/1‑151 (MQ=255)                                                                                   caGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCaaa                                                              >  1:176311/1‑150 (MQ=255)                                                                                   caGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCCTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCaaa                                                              >  1:591891/1‑150 (MQ=255)                                                                                                 gACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCg                                                 <  1:203206/149‑1 (MQ=255)                                                                                                 gACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGt                                                <  2:167596/150‑1 (MQ=255)                                                                                                    caccaGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGt                                                >  1:593548/1‑147 (MQ=255)                                                                                                         gagagagagTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGAc                                        >  2:589223/1‑150 (MQ=255)                                                                                                         gagagagagTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACa                                       <  2:543560/151‑1 (MQ=255)                                                                                                          agagagagTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACAt                                      <  2:350956/151‑1 (MQ=255)                                                                                                              agagTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCgg                                   >  2:330887/1‑150 (MQ=255)                                                                                                               gagTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGaa                                 <  2:497384/151‑1 (MQ=255)                                                                                                                 gTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAAc                               <  1:442946/151‑1 (MQ=255)                                                                                                                 gTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAAc                               <  1:438626/151‑1 (MQ=255)                                                                                                                       gCTCTGGTGCAATTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTc                         <  1:198322/151‑1 (MQ=255)                                                                                                                                  aTTTATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGaa              <  2:625515/151‑1 (MQ=255)                                                                                                                                     tATCTTTGGTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAacac           >  2:480254/1‑151 (MQ=255)                                                                                                                                             gTTTCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCa   <  1:83549/151‑1 (MQ=255)                                                                                                                                              tttCGCTCCACAATGCGTCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCa   <  1:465409/150‑1 (MQ=255)                                                                                                                                              tttCGCTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCAc  >  1:40560/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCAC  >  NC_000913/3907799‑3908090

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