Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,883,255	C→T	D159N (GAT→AAT)	dnaA ←	chromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,883,255	0	C	T	100.0%	205.6 / NA	54	D159N (GAT→AAT)	dnaA	chromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/26);  total (28/26)

TGAACAAAGCGCTCGGAGTGCATATAAACCACTTTGGCATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATCCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCG  >  NC_000913/3883107‑3883401                                                                                                                                                     |                                                                                                                                                   tGAACAAAGCGCTCGGAGTGCATATAAACCACTTTGGCATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCg                                                                                                                                                  >  1:244019/1‑151 (MQ=255)  gAACAAAGCGCTCGGAGTGCATATAAACCACTTTGGCATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGc                                                                                                                                                 >  2:39612/1‑151 (MQ=255)     cAAAGCGCTCGGAGTGCATATAAACCACTTTGGCATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCAc                                                                                                                                              >  1:561643/1‑151 (MQ=255)     cAAAGCGCTCGGAGTGCATATAAACCACTTTGGCATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCAc                                                                                                                                              >  2:284977/1‑151 (MQ=255)      aaaGCGCTCGGAGTGCATATAAACCACTTTGGCATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCAAcc                                                                                                                                             >  2:580406/1‑151 (MQ=255)        aGCGCCCGGAGTGCATATAAACCACTTTGGCATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCCCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCt                                                                                                                                            >  1:73583/1‑150 (MQ=255)              cGGAGTGCATATAAACCACTTTGGCATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGa                                                                                                                                       <  2:705510/149‑1 (MQ=255)              cGGAGTGCATATAAACCACTTTGGCATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGa                                                                                                                                       <  1:387820/149‑1 (MQ=255)                      atatAAACCACTTTGGCATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCcgcgc                                                                                                                             >  2:555529/1‑151 (MQ=255)                           aaCCACTTTGGCATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGcc                                                                                                                        >  2:108355/1‑151 (MQ=255)                           aaCCACTTTGGCATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGcc                                                                                                                        >  2:62053/1‑151 (MQ=255)                                      cATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGtt                                                                                                                >  2:671737/1‑148 (MQ=255)                                      cATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGtt                                                                                                                >  2:401770/1‑148 (MQ=255)                                        ttCGGCTTGCGCGCCATAATGACGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGAtt                                                                                                           >  2:527712/1‑151 (MQ=255)                                           ggCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTAc                                                                                                        <  2:228430/151‑1 (MQ=255)                                               tGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCtt                                                                                                     <  1:186963/150‑1 (MQ=255)                                               tGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTc                                                                                                    >  1:199360/1‑151 (MQ=255)                                                 cgcgcCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCa                                                                                                   >  2:344760/1‑150 (MQ=255)                                                 cgcgcCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCa                                                                                                   <  1:39581/150‑1 (MQ=255)                                                 cgcgcCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCACCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCa                                                                                                   >  1:267115/1‑150 (MQ=255)                                                   cgcCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACg                                                                                                >  2:671634/1‑151 (MQ=255)                                                     ccttAATGCCTTTACCCACCGCATGCAACAGGTGAGTTTTCCCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTTTTCGCCCCCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCt                                                                                                      <  1:242642/140‑1 (MQ=255)                                                     ccATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGa                                                                                               >  1:410592/1‑150 (MQ=255)                                                       aTAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGt                                                                                            >  2:196012/1‑151 (MQ=255)                                                            gCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATc                                                                                        <  1:555477/150‑1 (MQ=255)                                                                   cccACCGCATGCAGCAGGTGAGTTTTACCCAGCCCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACgtgt                                                                                <  2:771359/151‑1 (MQ=255)                                                                     cACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACgtgtgt                                                                              >  2:763414/1‑151 (MQ=255)                                                                     cACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACgtgtgt                                                                              <  2:514656/151‑1 (MQ=255)                                                                        cGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTg                                                                           <  1:543629/151‑1 (MQ=255)                                                                        cGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTg                                                                           <  2:342112/151‑1 (MQ=255)                                                                           aTGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACg                                                                        >  2:432214/1‑151 (MQ=255)                                                                              cagcagGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGttt                                                                     <  1:195979/151‑1 (MQ=255)                                                                               agcagGTGAGTTTTACCCAGACCCGTGCTGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTa                                                                    <  2:178916/151‑1 (MQ=255)                                                                                     tGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGtt                                                                >  2:432201/1‑149 (MQ=255)                                                                                               ccAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGt                                                    <  2:95381/151‑1 (MQ=255)                                                                                                cAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTc                                                   >  2:660003/1‑151 (MQ=255)                                                                                                 aGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCg                                                  <  2:25863/151‑1 (MQ=255)                                                                                                   aCCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGt                                                >  1:280509/1‑151 (MQ=255)                                                                                                    cccGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGtt                                               <  2:570423/151‑1 (MQ=255)                                                                                                          ccgccATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTgccg                                         <  2:134320/151‑1 (MQ=255)                                                                                                          ccgccATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTgccg                                         <  1:584243/151‑1 (MQ=255)                                                                                                               aTAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTgccggggcc                                    >  2:634560/1‑151 (MQ=255)                                                                                                                tAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTgccggggccg                                   >  1:55189/1‑151 (MQ=255)                                                                                                                     gAACAACGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACg                              <  1:35111/151‑1 (MQ=255)                                                                                                                       acaacGGGTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGtt                            >  2:443251/1‑151 (MQ=255)                                                                                                                            gggTTAGAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACGTTCAACGGAGTTATCACACGTGTGTTTGACGTTTACGTTAGAACGATAGGGCGGTTCTGCCGGGCCCGGGACGCTATccc                       >  1:209850/1‑151 (MQ=255)                                                                                                                             ggTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTGGGTTCTGCCGGGGCCGGGACGTTATCCCa                      <  1:505796/151‑1 (MQ=255)                                                                                                                             ggTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATccc                       >  2:273780/1‑150 (MQ=255)                                                                                                                             ggTTATAGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCa                      <  1:659945/151‑1 (MQ=255)                                                                                                                                   aGGCACCGCCAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGa                <  2:475025/151‑1 (MQ=255)                                                                                                                                           ccAGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAg        <  1:275520/151‑1 (MQ=255)                                                                                                                                             aGGGTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGaa      <  2:569474/151‑1 (MQ=255)                                                                                                                                              gggTTATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAg     <  2:750029/151‑1 (MQ=255)                                                                                                                                                  tATTCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGgca  <  1:656905/150‑2 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TGAACAAAGCGCTCGGAGTGCATATAAACCACTTTGGCATTCGGCTTGCGCGCCATAATGCCGTTACCCACCGCATGCAGCAGGTGAGTTTTACCCAGACCCGTGCCGCCATAAAGGAACAACGGGTTATAGGCACCGCCAGGGTTATCCGCCACCTGGCGAGCCGCCGCGCGCGCCAGTTGGTTAGATTTACCTTCAACGAAGTTATCAAACGTGTGTTTGACGTTTACGTTAGAACGATAGGTCGGTTCTGCCGGGGCCGGGACGTTATCCCAACCTGAGCGCGTAGAAGGCG  >  NC_000913/3883107‑3883401

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