Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	6,502	T→A	intergenic (‑43/+27)	yaaA ← / ← yaaJ	DNA binding and peroxide stress response protein YaaA/putative transporter YaaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	6,502	0	T	A	100.0%	70.5 / NA	23	intergenic (‑43/+27)	yaaA/yaaJ	DNA binding and peroxide stress response protein YaaA/putative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (8/15);  total (8/15)

CAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAATGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGA  >  NC_000913/6260‑6745                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    cAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTACTTGCAGGAAATTTATGCCGACTTTAGCAACAaaagagaatg                                                                                                                                                                                                                                              >  2:149819/1‑250 (MQ=255)             cgaccggcgTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAg                                                                                                                                                                                                                                 <  2:176321/247‑1 (MQ=255)                 gggctttGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTtgt                                                                                                                                                                                                                             <  1:159905/246‑1 (MQ=255)                  ggcgTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTtgtg                                                                                                                                                                                                                            >  1:89127/1‑251 (MQ=255)                    cgTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGAt                                                                                                                                                                                                                          <  1:207483/251‑1 (MQ=255)                             ccGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGc                                                                                                                                                                                                                 >  1:13825/1‑251 (MQ=255)                                    tttGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAAcat                                                                                                                                                                                                          >  2:78559/1‑251 (MQ=255)                                                   cattatcGTGCTAATTTGCGGAGTCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCACCGGCAGCGTATAGCGCGTGGTGGTGAACGCGCTTTGGTAATCAAGTGTTTTTGCAGGTGAAATAAGAAACAGCATATCCAGTCCTTGCAGGAAATTTATGTCGACTTTAGCAAAAAAAGAGAATGAGTCGAGCGAGAGTTGTGATTACTCTTGCGAAACATCATCCCGCGCGTCCg                                                                                                                                                                                           <  2:153797/245‑1 (MQ=255)                                                               aaTCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGc                                                                                                                                                                               <  2:194376/251‑1 (MQ=255)                                                                                                        gctgGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGgcg                                                                                                                                       <  1:232199/250‑1 (MQ=255)                                                                                                                      aaCAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATAGCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTAc                                                                                                                         <  2:89127/250‑1 (MQ=255)                                                                                                                                           cgcgTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGc                                                                                                   <  2:131583/251‑1 (MQ=255)                                                                                                                                                           ggCTTTGGTAATCAAGCGGTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCGGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTTCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACcaca                                                                                   <  1:71937/251‑1 (MQ=255)                                                                                                                                                                    aagcAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGTCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACATGCGAGAGc                                                                          <  2:103399/248‑1 (MQ=255)                                                                                                                                                                              tttCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATg                                                                 <  1:138093/250‑1 (MQ=255)                                                                                                                                                                                tCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGc                                                               <  1:86065/250‑1 (MQ=255)                                                                                                                                                                                        gAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCCAAtt                                                      >  1:139765/1‑251 (MQ=255)                                                                                                                                                                                              aGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGGAAAAAAAGAGAATGAGTGGATCGATAGTGGTGATTACTCCTGCGAAACATCAGCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACTCAACGGATCAAACACCGGTTTCATGCCGAGTGTAGTCTGTCGTAGATAATCACTGGCAATTGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAAt                                                <  1:155705/251‑1 (MQ=255)                                                                                                                                                                                                     gCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCCACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAAGAGcgga                                          >  1:8663/1‑247 (MQ=255)                                                                                                                                                                                                     gCATATCCAGTCCGTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTa                                                  <  2:8663/242‑1 (MQ=255)                                                                                                                                                                                                         atCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAg                                      <  1:26587/250‑1 (MQ=255)                                                                                                                                                                                                            cAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAACAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGAAATCCGGATAACGCAACGGATCAACCACCGGGCGCACGCCGAGCTTACGCTGCCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGAAGTAAAATGGCGGTCAAATTGGTAAAAGCCATGCCGGcc                                   >  1:135989/1‑250 (MQ=255)                                                                                                                                                                                                                                            cAAAAAAAGAGAATCAGTTGATCGATAGTTCTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTCCGCTGGCGTCGATAATCACTGGCAATGGTATGAACCACAGGCGCGAGCAGTAAAATGGCGGTCAAATTGGTAATCGCCAGGAAGGCCATTATGGTATCTGCCAGTTGCCACATCAGCgg   >  1:188198/1‑250 (MQ=255)                                                                                                                                                                                                                                             aaaaaaAGAGAATGAGTTGACAGATAGTTGTGATTACTCCTGCGCAACATCATCCCAAGCGTCCGGAGAAAGCTGACCGCCGATATCCGGATAACGCAACGGATCAAACCCCGGGAGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTCTGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTTGGAAAAGAAATGCAGGCCATTATCATATcatcccgttgcccaaacatcgga  >  2:139078/1‑228 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    CAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAATGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGA  >  NC_000913/6260‑6745

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