Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,842,808	G→A	A428V (GCG→GTG)	dnaA ←	chromosomal replication initiation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,842,808	0	G	A	100.0%	106.2 / NA	29	A428V (GCG→GTG)	dnaA	chromosomal replication initiation protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/16);  total (13/16)

GTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGA  >  REL606/3842735‑3842876                                                                          |                                                                     gTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAAcaca                                                                    <  1:714107/76‑1 (MQ=255)         ctTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAAt                                                            <  1:739877/76‑1 (MQ=255)           tCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAAtct                                                          <  2:1266979/76‑1 (MQ=255)                cAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGc                                                     <  2:379351/76‑1 (MQ=255)                cAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGc                                                     <  2:128918/76‑1 (MQ=255)                  aCTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGCAg                                                   >  2:1337725/1‑76 (MQ=255)                  aCTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGCAg                                                   >  1:175428/1‑76 (MQ=255)                   cTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGCAGa                                                  >  2:1002771/1‑76 (MQ=255)                     gCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACt                                                >  1:104823/1‑76 (MQ=255)                      cTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACtg                                               >  2:1090395/1‑76 (MQ=255)                        cGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCCCGGCCACCAAACACATCGCCAATCTCCCGCAGACtggg                                             >  2:443103/1‑74 (MQ=255)                        cGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACtgtg                                             >  2:133102/1‑76 (MQ=255)                           tCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGtt                                          >  1:209852/1‑76 (MQ=255)                                    aggcaTGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGctc                                 >  2:1370425/1‑76 (MQ=255)                                       caTGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCttt                              <  2:523006/76‑1 (MQ=255)                                       caTGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCttt                              <  2:1142414/76‑1 (MQ=255)                                        aTGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTc                             <  2:104375/76‑1 (MQ=255)                                        aTGAAGCACCGTCGTGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTc                             <  2:378675/76‑1 (MQ=255)                                                 cgtcgtGTGGTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGcgcc                    >  1:940646/1‑76 (MQ=255)                                                     gtgtGGTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGcgccatcg                <  2:668239/76‑1 (MQ=255)                                                      tgtgGTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGcgccatcgc               >  2:725047/1‑76 (MQ=255)                                                       gtgGTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGcgccatcgcc              >  1:1276256/1‑76 (MQ=255)                                                         ggTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGcgccatcgccat            >  1:168655/1‑76 (MQ=255)                                                          gTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCgccatcgccatc           <  1:1090395/76‑1 (MQ=255)                                                          gTCACGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCgccatcgccatc           <  1:1498660/76‑1 (MQ=255)                                                           tCACGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCt          >  2:672727/1‑76 (MQ=255)                                                             aCGGCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTgg        <  1:495972/76‑1 (MQ=255)                                                               ggCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGgcg      <  2:239396/76‑1 (MQ=255)                                                               ggCCACCAAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGgcg      <  2:181054/76‑1 (MQ=255)                                                                 ccaccaAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGcgcg    <  2:574406/76‑1 (MQ=255)                                                                  caccaAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCgg   <  1:657424/76‑1 (MQ=255)                                                                   cccaAACACATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGa  <  2:1241431/75‑1 (MQ=255)                                                                          |                                                                     GTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGA  >  REL606/3842735‑3842876

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