Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,601,154	C→T	A122T (GCC→ACC)	dnaT ←	primosomal protein DnaT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,601,154	0	C	T	100.0%	123.0 / NA	35	A122T (GCC→ACC)	dnaT	primosomal protein DnaT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/12);  total (23/12)

AGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGCCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTT  >  NC_000913/4601087‑4601231                                                                    |                                                                              aGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCggtg                                                                  >  2:883712/1‑81 (MQ=255)     tttGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCa                                                              >  2:883720/1‑81 (MQ=255)       tgttgCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCAcc                                                            >  1:883718/1‑81 (MQ=255)        gttgCCACTGCACATGGTGAAAGACTTTACGTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTTACCg                                                           >  2:883710/1‑81 (MQ=255)        gttgCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCg                                                           >  1:883705/1‑81 (MQ=255)            ccACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGctc                                                       >  1:883707/1‑81 (MQ=255)            ccACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGct                                                        <  1:883666/80‑1 (MQ=255)              aCTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGctctc                                                     <  2:883666/81‑1 (MQ=255)               cTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGctctct                                                    <  1:883668/81‑1 (MQ=255)                     aTGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTCAgcgc                                              >  2:883699/1‑81 (MQ=255)                       ggTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCa                                            >  1:883726/1‑81 (MQ=255)                       ggTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCa                                            >  1:883709/1‑81 (MQ=255)                             aGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGccc                                       >  1:883706/1‑80 (MQ=255)                              gACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCa                                     <  2:883668/81‑1 (MQ=255)                                cTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATa                                   <  1:883670/81‑1 (MQ=255)                                  ttACCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAgc                                 >  2:883711/1‑81 (MQ=255)                                    aCCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAgcgc                               >  1:883714/1‑81 (MQ=255)                                    aCCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAgcgc                               >  2:883724/1‑81 (MQ=255)                                    aCCTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAgcg                                <  1:883671/80‑1 (MQ=255)                                     ccTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGcgcc                              <  2:883671/81‑1 (MQ=255)                                     ccTTCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGcgcc                              <  2:883672/81‑1 (MQ=255)                                        tCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGcgccgc                            >  2:883707/1‑80 (MQ=255)                                        tCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGcgccgc                            >  2:883709/1‑80 (MQ=255)                                        tCCGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGcgccgc                            >  1:883716/1‑80 (MQ=255)                                          cGCCTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAg                         >  1:883737/1‑81 (MQ=255)                                            ccTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGc                       >  1:883719/1‑81 (MQ=255)                                            ccTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGc                       <  1:883674/81‑1 (MQ=255)                                             cTGCCAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCg                      >  2:883723/1‑81 (MQ=255)                                                ccAGTAGGCAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCgg                     <  2:883674/79‑1 (MQ=255)                                                        cAATGAATGAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCgg           >  2:883738/1‑81 (MQ=255)                                                         aatgaatgAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGTCCCATAGCGCCGCCAGGCGGATAAAATCgg           >  2:883714/1‑80 (MQ=255)                                                         aatgaatgAGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCgg           >  1:883734/1‑80 (MQ=255)                                                          atgaatgaGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCa         >  2:883719/1‑81 (MQ=255)                                                             aatgaGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATc       <  1:883676/80‑1 (MQ=255)                                                                 aGGTCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGtt  <  2:883676/81‑1 (MQ=255)                                                                    |                                                                              AGTTTTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGCCAGTTCTTCGGTGGTCACCGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTT  >  NC_000913/4601087‑4601231

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