Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,909,590	G→A	A270V (GCG→GTG)	pstC ←	phosphate ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,909,590	0	G	A	100.0%	174.0 / NA	54	A270V (GCG→GTG)	pstC	phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/26);  total (28/26)

AGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACAC  >  NC_000913/3909460‑3909722                                                                                                                                   |                                                                                                                                     aGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTg                                                                                                                              >  1:985670/1‑139 (MQ=255)       gcgcATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTg                                                                                                                        >  2:910373/1‑139 (MQ=255)       gcgcATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTg                                                                                                                        >  1:1593525/1‑139 (MQ=255)         gcATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGtt                                                                                                                      <  2:1385799/139‑1 (MQ=255)           aTAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGc                                                                                                                    <  2:1300878/139‑1 (MQ=255)                       ttCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACa                                                                                                        >  2:1525110/1‑139 (MQ=255)                         cGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCGGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGc                                                                                                      >  2:928931/1‑139 (MQ=255)                                  gAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTg                                                                                             <  2:1538619/139‑1 (MQ=255)                                          tGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGc                                                                                                               >  1:67236/1‑113 (MQ=255)                                          tGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGc                                                                                                               <  2:67236/113‑1 (MQ=255)                                            aaGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGt                                                                                   <  2:892800/139‑1 (MQ=255)                                               gTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAgg                                                                                >  1:1245613/1‑139 (MQ=255)                                                  aTCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGgtgt                                                                             >  2:1512900/1‑139 (MQ=255)                                                      caAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTAcc                                                                         <  2:1245613/139‑1 (MQ=255)                                                      caAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTAcc                                                                         <  1:470735/139‑1 (MQ=255)                                                      caAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTAcc                                                                         <  1:1512900/139‑1 (MQ=255)                                                         aCAGGATCAGGCCCAGTTGCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGa                                                                                         >  2:1059069/1‑120 (MQ=255)                                                         aCAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGa                                                                                         <  1:1059069/120‑1 (MQ=255)                                                         aCAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCgat                                                                      >  2:570967/1‑139 (MQ=255)                                                             gATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCATACCGGATTCCGAGTCCGCAAATTCGTTCGCCAGCACAGAGGTAATAATGTATCCCGGCATATACAGCAAGGCGCAGTCGAGCTGATAGGTGTTACAGataata                                                                  >  1:487706/1‑139 (MQ=255)                                                                     ccAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCAc                                                          >  2:1410298/1‑139 (MQ=255)                                                                        gTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGc                                                                                                               <  1:720903/83‑1 (MQ=255)                                                                        gTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGc                                                                                                               >  2:720903/1‑83 (MQ=255)                                                                         ttCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGc                                                                                                                    >  2:520995/1‑77 (MQ=255)                                                                         ttCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGc                                                                                                                    <  1:520995/77‑1 (MQ=255)                                                                                 gTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACag                                                                                                                                   >  2:341834/1‑54 (MQ=255)                                                                                 gTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACag                                                                                                                                   <  1:341834/54‑1 (MQ=255)                                                                                  tGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACag                                                                                                                                   <  2:344622/53‑1 (MQ=255)                                                                                  tGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACag                                                                                                                                   >  1:344622/1‑53 (MQ=255)                                                                                  tGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCa                                             >  1:1081126/1‑139 (MQ=255)                                                                                   gCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGt                                                                                   >  1:1067922/1‑100 (MQ=255)                                                                                   gCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGt                                                                                   <  2:1067922/100‑1 (MQ=255)                                                                                   gCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCAc                                            >  2:256982/1‑139 (MQ=255)                                                                                   gCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCAc                                            >  1:355479/1‑139 (MQ=255)                                                                                    cGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCAcc                                           >  1:231789/1‑139 (MQ=255)                                                                                           gTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAgcgcg                                    <  2:1577208/139‑1 (MQ=255)                                                                                            tGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAgcgcgc                                   >  1:779239/1‑139 (MQ=255)                                                                                                 aCCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGccc                              <  1:171153/139‑1 (MQ=255)                                                                                                   cGGATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcccag                            >  2:586423/1‑139 (MQ=255)                                                                                                    ggATTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcccagc                           <  2:985670/139‑1 (MQ=255)                                                                                                      aTTCCGCTTCCGCAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGcccagccc                         >  2:1431267/1‑139 (MQ=255)                                                                                                            cttccgcAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCat                   >  2:785162/1‑139 (MQ=255)                                                                                                             ttccgcAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGc                                                       >  2:806711/1‑102 (MQ=255)                                                                                                             ttccgcAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGc                                                       <  1:806711/102‑1 (MQ=255)                                                                                                               ccgcAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGc                                                                                                                    >  2:1361658/1‑39 (MQ=38)                                                                                                               ccgcAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGc                                                                                                                    <  1:1361658/39‑1 (MQ=255)                                                                                                                cgcAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCatgatg               <  1:256982/139‑1 (MQ=255)                                                                                                                cgcAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCatgatg               <  2:101809/139‑1 (MQ=255)                                                                                                                cgcAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCatgatg               >  1:1296820/1‑139 (MQ=255)                                                                                                                  cAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATgcc             <  2:355479/139‑1 (MQ=255)                                                                                                                  cAAATTCGTTCGCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATgcc             <  2:1081126/139‑1 (MQ=255)                                                                                                                       tcgttcgCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGa        <  2:325936/139‑1 (MQ=255)                                                                                                                           tcgCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAac    <  2:1296820/139‑1 (MQ=255)                                                                                                                             gCCAGCACAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAacac  <  1:570967/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                     AGCCAGGCGCATAATCATAAACTTCGATGCGGCGAGGACGATGAAGGTAATCACAAACAGGATCAGGCCCAGTTCCATCAGTGCGGCAACGTGCAGACCGGATTCCGCTTCCGCAAATTCGTTCGCCAGCGCAGAGGTGATACTGTTGCCCGGCATATACAGCGAGGCGCTGTCGAGCTGGTAGGTGTTACCGATGATAAAGGTCACCGCCATGGTTTCACCGAGCGCGCGGCCCAGCCCCAGCATGATGCCGCCGATAACAC  >  NC_000913/3909460‑3909722

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