Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,948,026	C→T	intergenic (‑59/‑60)	hdfR ← / → maoP	DNA‑binding transcriptional dual regulator HdfR/macrodomain Ori protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,948,026	0	C	T	100.0%	108.8 / NA	30	intergenic (‑59/‑60)	hdfR/maoP	DNA‑binding transcriptional dual regulator HdfR/macrodomain Ori protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/14);  total (16/14)

TTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGAAAGCTTT  >  NC_000913/3947949‑3948100                                                                              |                                                                           ttAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAg                                                                         >  1:772436/1‑81 (MQ=255)      aaTTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTaa                                                                    >  2:772443/1‑81 (MQ=255)            gTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATa                                                              >  2:772452/1‑81 (MQ=255)              aTCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACt                                                            >  2:772432/1‑81 (MQ=255)                ccACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTcc                                                          >  1:772439/1‑81 (MQ=255)                    aGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCgcgc                                                      >  1:772450/1‑81 (MQ=255)                        ctctCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATa                                                  <  1:772408/81‑1 (MQ=255)                             cGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACt                                               <  1:772409/79‑1 (MQ=255)                              gCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGc                                            >  1:772463/1‑81 (MQ=255)                                aaaaTTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTc                                          <  1:772410/81‑1 (MQ=255)                                aaaaTTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTc                                          <  2:772411/81‑1 (MQ=255)                                aaaaTTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTc                                          >  1:772441/1‑81 (MQ=255)                                 aaaTTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCg                                         >  1:772461/1‑81 (MQ=255)                                 aaaTTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCg                                         <  1:772412/81‑1 (MQ=255)                                  aaTTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCgg                                        >  2:772455/1‑81 (MQ=255)                                    ttattTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGt                                       >  1:772464/1‑80 (MQ=255)                                      attTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCa                                     >  1:772449/1‑80 (MQ=255)                                             gTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAAt                               >  1:772445/1‑79 (MQ=255)                                             gTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAAt                               >  1:772444/1‑79 (MQ=255)                                             gTTATGATTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTa                             <  2:772413/81‑1 (MQ=255)                                                   aTTTAAATGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGc                        <  2:772414/80‑1 (MQ=255)                                                       aaaTGTTTTGTTTTACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGc                    <  1:772415/80‑1 (MQ=255)                                                            ttttgttttACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGAt               <  1:772416/80‑1 (MQ=255)                                                               tgttttACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTTGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCg           >  2:772453/1‑81 (MQ=255)                                                                 ttttACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCgg          <  1:772417/80‑1 (MQ=255)                                                                   ttACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGa         >  2:772454/1‑79 (MQ=255)                                                                    tACACTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGaaa       <  2:772418/80‑1 (MQ=255)                                                                     acacTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGAAAg      <  2:772419/80‑1 (MQ=255)                                                                      cacTCTGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGAAAGc     <  1:772420/80‑1 (MQ=255)                                                                        ctctGTTAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGAAAGCttt  <  2:772421/81‑1 (MQ=255)                                                                              |                                                                           TTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGAAAGCTTT  >  NC_000913/3947949‑3948100

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