Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,145	G→T	intergenic (+113/‑489)	cra → / → mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,145	0	G	T	100.0%	65.8 / NA	19	intergenic (+113/‑489)	cra/mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (7/12);  total (7/12)

GCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATT  >  NC_000913/88921‑89386                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  gCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACt                                                                                                                                                                                                                         <  1:52202/251‑1 (MQ=255)   gtgtGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTc                                                                                                                                                                                                                        >  2:219046/1‑250 (MQ=255)     gtgCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCAc                                                                                                                                                                                                                      <  1:539/250‑1 (MQ=255)              aTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAAt                                                                                                                                                                                                            >  2:163402/1‑251 (MQ=255)                                     aCCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCg                                                                                                                                                                                      >  2:81432/1‑250 (MQ=255)                                        gcgTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATGAAGTGATATCGACAg                                                                                                                                                                                  <  2:240646/251‑1 (MQ=255)                                                  aaaCCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTg                                                                                                                                                                        <  1:255592/251‑1 (MQ=255)                                                         gTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCattatt                                                                                                                                                                 <  1:54450/251‑1 (MQ=255)                                                                         aCGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATgcagaga                                                                                                                                                 <  1:55351/251‑6 (MQ=255)                                                                           gTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATgcagagatg                                                                                                                                               >  1:75625/1‑244 (MQ=255)                                                                                  ctATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAAt                                                                                                                                            <  1:238176/247‑1 (MQ=255)                                                                                                      gCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCt                                                                                                                    <  1:265505/251‑1 (MQ=255)                                                                                                           aGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCt                                                                                                               <  1:231891/251‑1 (MQ=255)                                                                                                                                         gTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCca                                                                                 <  1:179455/251‑1 (MQ=255)                                                                                                                                                                      ggaatGGACATCAATATCCTTATAGTAAATAAGACTATTCCTGACGATTCTTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCAAAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCACGCGCCACACAGACGCAACAGGGACTGTTAACCggg                                                    >  1:110874/5‑251 (MQ=255)                                                                                                                                                                                                aaaTAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTcg                          <  2:186960/251‑1 (MQ=255)                                                                                                                                                                                                                  cTATGATTGATAAGTTCTTTTAAACCCGCCCGGTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCAGTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAATCCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGGTAGTCGGGGAAGATAGGTCCTAAAATGCCGCTCGCGTCGCAAACTg              <  1:170350/245‑1 (MQ=255)                                                                                                                                                                                                                   tattatTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACttt       >  1:103727/1‑251 (MQ=255)                                                                                                                                                                                                                        ttGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATAtt  >  2:51746/1‑251 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  GCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATT  >  NC_000913/88921‑89386

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