Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,601,522	T→C	intergenic (‑5/+102)	dnaT ← / ← yjjB	primosomal protein DnaT/putative succinate exporter YjjB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,601,522	0	T	C	100.0%	107.3 / NA	31	intergenic (‑5/+102)	dnaT/yjjB	primosomal protein DnaT/putative succinate exporter YjjB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/14);  total (17/14)

GCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGAGGCTGTCTC  >  NC_000913/4601444‑4601596                                                                               |                                                                           gCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGctt                                                                          >  2:955564/1‑81 (MQ=255)  ccAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTtc                                                                         <  2:955533/81‑1 (MQ=255)    aGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTtctc                                                                       <  2:955534/81‑1 (MQ=255)      aaCGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCaa                                                                     <  2:955535/81‑1 (MQ=255)           tttGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATg                                                                >  1:955577/1‑81 (MQ=255)               gTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAAc                                                            <  1:955537/81‑1 (MQ=255)               gTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAAc                                                            <  1:955536/81‑1 (MQ=255)                  aTCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGAt                                                         >  1:955587/1‑81 (MQ=255)                   tCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATg                                                        <  2:955538/81‑1 (MQ=255)                     gtgtACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGaa                                                      >  1:955584/1‑81 (MQ=255)                     gtgtACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGa                                                       <  2:955539/80‑1 (MQ=255)                        tACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAACg                                                    <  1:955540/80‑1 (MQ=255)                         aCCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAAc                                                     >  2:955570/1‑78 (MQ=255)                         aCCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAAc                                                     >  1:955586/1‑78 (MQ=255)                              ggCGTCAATACCAACGACGTCCGGGGTCAAAGCTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAggg                                             <  2:955541/81‑1 (MQ=255)                                 gTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGgaga                                          <  2:955543/81‑1 (MQ=255)                                 gTCAATACCAACGACGTCCGAGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGgaga                                          <  1:955542/81‑1 (MQ=255)                                    aaTACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATa                                        <  1:955544/80‑1 (MQ=255)                                       aCCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACtt                                    >  1:955569/1‑81 (MQ=255)                                           acgacgTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTta                                 >  1:955572/1‑80 (MQ=255)                                             gacgTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAg                              >  2:955568/1‑81 (MQ=255)                                             gacgTCCGGGGTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAg                              >  1:955580/1‑81 (MQ=255)                                                    ggggTCAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAacaca                       >  1:955596/1‑81 (MQ=255)                                                         cAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACggag                  >  1:955583/1‑81 (MQ=255)                                                         cAAAACTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACggag                  >  1:955573/1‑81 (MQ=255)                                                            aaCTCTGGAAGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACggag                  >  1:955592/1‑78 (MQ=255)                                                               tctGGAAGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGgaggagag             >  1:955601/1‑80 (MQ=255)                                                                    aaGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGAGGCt       >  1:955576/1‑81 (MQ=255)                                                                     aGACATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGAGGCTg      <  2:955545/81‑1 (MQ=255)                                                                       aCATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGAGGCTGtc    <  2:955546/81‑1 (MQ=255)                                                                       aCATACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGAGGCTGtc    >  2:955579/1‑81 (MQ=255)                                                                         aTACTGCTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGAGGCTGtctc  >  2:955578/1‑81 (MQ=255)                                                                               |                                                                           GCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGAGGCTGTCTC  >  NC_000913/4601444‑4601596

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: ‑