Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	7,925	+A	coding (35/1431 nt)	yaaJ ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	7,922	1	.	A	100.0%	94.3 / NA	30	W13L (TGG→TTG)	yaaJ	putative transporter
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (15/15);  total (15/15)

ACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCC‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTC  >  NC_000913/7789‑8054                                                                                                                                       |                                                                                                                                     aCCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCC‑aaaag                                                                                                                                 >  2:351902/1‑137 (MQ=255)   cGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCC‑aaaagga                                                                                                                               <  2:992834/139‑5 (MQ=255)   cGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCC‑aaaagga                                                                                                                               <  1:1104286/139‑5 (MQ=255)             tgGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAAt                                                                                                                      >  1:214447/1‑139 (MQ=255)                         aTTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCC‑aaaaggac                                                                                                                              <  1:1232676/118‑6 (MQ=255)                         aTTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCC‑aaaaggac                                                                                                                              >  2:1232676/1‑113 (MQ=255)                                aGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgc                                                                                                   <  1:351902/139‑1 (MQ=255)                                    ttttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgc                                                                                               >  2:414610/1‑139 (MQ=255)                                    ttttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgc                                                                                               >  2:411376/1‑139 (MQ=255)                                         ccAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                          <  1:1001088/139‑1 (MQ=255)                                                                     caAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              <  2:632600/139‑1 (MQ=255)                                                                          ccGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCtgtg                                                         <  1:411376/139‑1 (MQ=255)                                                                                             aacaacCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGttt                                      <  2:237473/139‑1 (MQ=255)                                                                                                 acCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaaga                                  >  2:1265701/1‑139 (MQ=255)                                                                                                      cgcCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                              <  1:208529/42‑1 (MQ=255)                                                                                                      cgcCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                              >  2:208529/1‑42 (MQ=255)                                                                                                      cgcCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTAcc                                              <  2:1136649/122‑1 (MQ=255)                                                                                                      cgcCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTAcc                                              >  1:1136649/1‑122 (MQ=255)                                                                                                       gcCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGtata                            >  1:1344056/1‑139 (MQ=255)                                                                                                        ccGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGc                                                                                                                             >  2:1143950/1‑41 (MQ=255)                                                                                                        ccGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGc                                                                                                                             <  1:1143950/41‑1 (MQ=255)                                                                                                        ccGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGtatat                           <  2:23967/139‑1 (MQ=255)                                                                                                                           ttACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGAt        <  1:1265701/139‑1 (MQ=255)                                                                                                                             aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              <  2:1142929/83‑1 (MQ=255)                                                                                                                             aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              >  1:1142929/1‑83 (MQ=255)                                                                                                                             aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATaa      >  2:20726/1‑139 (MQ=255)                                                                                                                             aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATaa      >  2:939732/1‑139 (MQ=255)                                                                                                                                atcccca‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              >  2:290941/4‑79 (MQ=255)                                                                                                                                atcccca‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              <  1:290941/76‑1 (MQ=255)                                                                                                                                atcccca‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACg                   >  2:255756/4‑122 (MQ=255)                                                                                                                                 aTCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTgg                                                                                                    <  2:274809/41‑1 (MQ=255)                                                                                                                                 aTCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTgg                                                                                                    >  1:274809/1‑41 (MQ=255)                                                                                                                                 aTCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCata                                                                                            <  1:180727/49‑1 (MQ=255)                                                                                                                                 aTCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCata                                                                                            >  2:180727/1‑49 (MQ=255)                                                                                                                                 aTCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACg                   <  1:255756/122‑1 (MQ=255)                                                                                                                                 aTCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTc  >  1:406116/1‑139 (MQ=255)                                                                                                                                  tCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                          <  1:499914/50‑1 (MQ=255)                                                                                                                                  tCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                          >  2:499914/1‑50 (MQ=255)                                                                                                                                  cccca‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaag                                   <  1:1097722/103‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     ACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCC‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTC  >  NC_000913/7789‑8054

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