Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	6,898	A→G	L354L (CTT→CTC)	yaaJ ←	putative transporter YaaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	6,898	0	A	G	100.0%	49.0 / NA	18	L354L (CTT→CTC)	yaaJ	putative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/5);  total (13/5)

AGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGC  >  NC_000913/6669‑7126                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                     agagCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGAAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTAAGCAccc                                                                                                                                                                                                                 >  2:169806/1‑251 (MQ=255)                  gTCAAATGGGTAATAGCCAGGCAGGCCATTATGAGATCTGCCAGTTGCCACATCAGCGGAAGACGTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAacac                                                                                                                                                                                                <  2:13825/251‑1 (MQ=255)                                  ccATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTgg                                                                                                                                                                                >  1:222525/1‑251 (MQ=255)                                            cATTATGATATCTGACAGTTGCCACATCAGCGGAAGACTAAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGCTTTCGGCGTAAATGGAGTTGGCAACGATGGCGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGa                                                                                                                                                                       >  1:53750/1‑250 (MQ=255)                                               tATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATAc                                                                                                                                                                   >  1:38053/1‑251 (MQ=255)                                                 tGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGCACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCt                                                                                                                                                                 >  1:114388/1‑251 (MQ=255)                                                      gcTGCCAGTTGCCACCTCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGAACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTcc                                                                                                                                                              >  2:154318/2‑249 (MQ=255)                                                                                                              gAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAAt                                                                                                                                                                     <  2:48520/186‑1 (MQ=255)                                                                                                              gAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAAt                                                                                                                                                                     >  1:48520/1‑186 (MQ=255)                                                                                                                            gCAAACACAAGATCGTTTTGGGTTTTTTCAGGCAAAAAAAGAAGAGATTGTTTTCGGGTAAAATGTAGTGGACAACGATGGAGCTGAAGGCAAACAAAATAACCACGAGGTTAAAAACCCCAGCACCCCAGGAACCCATAAACCCCCCCATCGCCTTCGGGATAAGCTGAATCCCTTCCAGCGGCATGTGGTTTGGGCGTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGgtg                                                                                      <  1:169806/251‑1 (MQ=255)                                                                                                                                acacAAGATCGATTTCGGGGTGTACAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTCAATGTAGTTGGCACCGATGGAGCTGAAGGCAAACAGAATACACACGAGGGTAACAAACTCAGCACCCCAGGAGCACCTTAACACACGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGATGGTGCCGTTACCCGCCAATAATATCAGCATGTCGCTTGCCGTCCAGATGACCAggcggtcgc                                                                                  >  2:244551/1‑244 (MQ=255)                                                                                                                                                gggTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCAAGGCGCTTGCCGTACAGATGACCCGGGTGTCGATAAAAATGCCAatcat                                                                  >  2:162903/1‑251 (MQ=255)                                                                                                                                                                aaaaGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTCCAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTgc                                                  >  1:11383/1‑251 (MQ=255)                                                                                                                                                                      agagATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGcccccc                                                                                                                                                                                                                >  2:46659/1‑87 (MQ=255)                                                                                                                                                                      agagATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCAcccc                                                                                                                                                                                                                <  1:46659/87‑1 (MQ=255)                                                                                                                                                                       gagaTTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGAGCCGTTACCCGCCAg                                                                                                                              >  1:175523/1‑168 (MQ=255)                                                                                                                                                                       gagaTTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGGTGTGCCGTAACCCGCCAg                                                                                                                              <  2:175523/168‑1 (MQ=255)                                                                                                                                                                                                                tGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCCGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGc  >  2:101614/1‑251 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                     AGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACGCCGCTGCCGCTGCCGCGTTTGGC  >  NC_000913/6669‑7126

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