Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,909,781	A→G	G206G (GGT→GGC)	pstC ←	phosphate ABC transporter membrane subunit PstC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,909,781	0	A	G	100.0%	121.5 / NA	37	G206G (GGT→GGC)	pstC	phosphate ABC transporter membrane subunit PstC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/18);  total (19/18)

GTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGA  >  NC_000913/3909646‑3909917                                                                                                                                        |                                                                                                                                         gtTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGt                                                                                                                                      >  1:408649/1‑140 (MQ=255)     ccGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGc                                                                                                                                  >  2:369212/1‑140 (MQ=255)        atgatAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGa                                                                                                                               <  1:679817/140‑1 (MQ=255)         tgaaAAAGGCCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGAc                                                                                                                              <  2:720502/140‑1 (MQ=255)          gatAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGAct                                                                                                                             <  2:408649/140‑1 (MQ=255)          gatAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGAct                                                                                                                             <  1:382427/140‑1 (MQ=255)            tAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCg                                                                                                                                       <  1:520181/128‑1 (MQ=255)            tAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCg                                                                                                                                       >  2:520181/1‑128 (MQ=255)            tAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGActct                                                                                                                           <  2:830872/140‑1 (MQ=255)                 gTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTtcat                                                                                                                      <  1:705025/140‑1 (MQ=255)                    aCCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGa                                                                                                                               >  2:668073/1‑128 (MQ=255)                    aCCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGa                                                                                                                               <  1:668073/128‑1 (MQ=255)                    aCCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTtcatcat                                                                                                                   >  1:81520/1‑140 (MQ=255)                          aTGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCgg                                                                                                             >  1:135891/1‑140 (MQ=255)                              tttCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGtt                                                                                                         <  1:522056/140‑1 (MQ=255)                                   ccGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTc                                                                                                    >  2:171644/1‑140 (MQ=255)                                      aGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGaa                                                                                                 <  1:354408/140‑1 (MQ=255)                                           gcgGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACAt                                                                                            <  2:792267/140‑1 (MQ=255)                                                       cagcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAccg                                                                                <  1:171644/140‑1 (MQ=255)                                                         gcATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTAccgcc                                                                              >  2:149617/1‑140 (MQ=255)                                                           atgatgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACcgccgc                                                                            >  1:849781/1‑140 (MQ=255)                                                              atgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCa                                                                           >  1:298384/1‑138 (MQ=255)                                                              atgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCa                                                                           <  2:298384/138‑1 (MQ=255)                                                              atgCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAAt                                                                         >  2:181445/1‑140 (MQ=255)                                                                 ccgccgATACCCCCTTTTTTGGTGAACGGAAGAACGATCCGCCAGATAACTTCCCGGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTa                                                                      <  2:151100/140‑1 (MQ=255)                                                                 ccgccgATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCCCCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTa                                                                      <  1:369212/140‑1 (MQ=255)                                                                 ccgccgATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCCCGCATTACCGCCGCAATGTa                                                                      <  2:135891/140‑1 (MQ=255)                                                                  cgccgATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTAc                                                                     >  2:529325/1‑140 (MQ=255)                                                                            cccATTTTGGGTGAACGGAAGAACGATACGCAAGATACCTCCCCAGGGGGTGCAGCCAATGCCGTAGGCCGACTCTTTCACCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCa                                                           <  2:849781/139‑1 (MQ=255)                                                                                              aagaaCGATACGCCAGATACCTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTTTTCGAACACATCCCGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGc                                         <  2:780137/140‑1 (MQ=255)                                                                                                    gATCCGCCAGATACCTCCCCGGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCCCCGGGTTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCcg                                   <  2:340311/140‑1 (MQ=255)                                                                                                        cGCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAg                               >  2:324130/1‑140 (MQ=255)                                                                                                         gCCAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAgg                              <  1:324130/140‑1 (MQ=255)                                                                                                           cAGATAACTTCCCAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGgat                            >  1:562136/1‑140 (MQ=255)                                                                                                                     cccAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAg                               >  1:840212/1‑127 (MQ=255)                                                                                                                     cccAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAg                               <  2:840212/127‑1 (MQ=255)                                                                                                                     cccAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCa                  >  1:213840/1‑140 (MQ=255)                                                                                                                     cccAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCa                  <  2:562136/140‑1 (MQ=255)                                                                                                                       cAGGTGGTGCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCaaa                >  2:130886/1‑140 (MQ=255)                                                                                                                              tgCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCgggg         >  2:422873/1‑140 (MQ=255)                                                                                                                              tgCAGCCAATGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCACCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCgggg         <  1:831028/140‑1 (MQ=255)                                                                                                                                     aaTGCCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAgaga  >  2:401135/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                         GTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACACCATTTTTGGTGAACGGAAGAACGATACGCCAGATAACTTCCCAGGTGGTGCAGCCAATACCGTAGGCCGACTCTTTCATCATCACCGGGGTTTGTTCGAACACATCACGCATTACCGCCGCAATGTACGGAATAATCATGATGGCGAGGATCACGCCTGCCGCGAGGATACCGATACCAAATGCGGGGCCAGAGA  >  NC_000913/3909646‑3909917

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