Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,331,722	G→A	A407A (GCG→GCA)	proP →	osmolyte:H(+) symporter ProP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,331,722	0	G	A	100.0%	81.6 / NA	25	A407A (GCG→GCA)	proP	osmolyte:H(+) symporter ProP
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/14);  total (11/14)

CTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTAGTGGC‑GGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAA  >  NC_000913/4331602‑4331841                                                                                                                         |                                                                                                                           cTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCt                                                                                                         >  2:519933/1‑140 (MQ=255)    tACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCACTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGc                                                                                                      >  1:96889/1‑140 (MQ=255)       gggCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGa                                                                                                   >  2:291234/1‑140 (MQ=255)             tATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCtga                                                                                             <  1:80018/140‑1 (MQ=255)                ggCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCtgatga                                                                                          <  1:519933/140‑1 (MQ=255)                    tCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcc                                                                                      <  1:168935/140‑1 (MQ=255)                    tCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCACTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcc                                                                                      >  1:467254/1‑140 (MQ=255)                             ccAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTa                                                                             >  1:208780/1‑140 (MQ=255)                                   aTGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGAt                                                                     <  2:538083/140‑1 (MQ=255)                                       tCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTa                                                                 <  1:397196/140‑1 (MQ=255)                                                 tatCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTAGTGGC‑ggtgg                                                      <  1:161503/140‑1 (MQ=255)                                                        tACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTGGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTAGTGGC‑GGTTGCTTGttt                                               >  1:531365/1‑140 (MQ=255)                                                           aGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑TCTGATGGTAGTGGC‑GGTGGTTGGTTTAAt                                            >  1:493141/1‑140 (MQ=255)                                                             cgcgcTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTAGTGGC‑GGTGGTTGGTTTAATCa                                          <  2:437177/140‑1 (MQ=255)                                                               cgcTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTAGTGGCGGGGGGTTGGTTTatcatc                                         >  2:322432/1‑135 (MQ=255)                                                               cgcTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTAGTGGC‑GGGGGGTGGTTTTAttaac                                        >  1:419733/1‑136 (MQ=255)                                                               cgcTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTAGTGGC‑GGGGGGTGGGTTTAttaac                                        >  2:92044/1‑136 (MQ=255)                                                                                      tCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTAGTGGC‑GGTGGGTGGGTTTATTACCCGCGTGACCAATAAAAAGAAggg                 >  2:80447/1‑139 (MQ=255)                                                                                           gCTGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTAGTGGC‑GGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATc            <  1:92044/140‑1 (MQ=255)                                                                                             tGGTTGCCGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTAGTGGC‑GGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGt          <  1:322432/140‑1 (MQ=255)                                                                                                   ccGGTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCCCTGATGGTAGTGGC‑GGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGt      <  1:332568/140‑1 (MQ=255)                                                                                                     ggTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑TCTGATGGTAGTGGC‑GGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGaa  <  2:493141/140‑1 (MQ=255)                                                                                                     ggTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTAGTGGC‑GGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGaa  <  2:208780/140‑1 (MQ=255)                                                                                                     ggTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTAGTGGC‑GGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGaa  <  1:291234/140‑1 (MQ=255)                                                                                                     ggTCTGACGCCAACGCTGGCAGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTAGTGGC‑GGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGaa  <  1:22947/140‑1 (MQ=255)                                                                                                                         |                                                                                                                           CTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA‑‑CCTGATGGTAGTGGC‑GGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAA  >  NC_000913/4331602‑4331841

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