Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,296,292	C→T	R162C (CGC→TGC)	diaA →	DnaA initiator‑associating factor for replication initiation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,296,292	0	C	T	100.0%	194.7 / NA	52	R162C (CGC→TGC)	diaA	DnaA initiator‑associating factor for replication initiation
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/30);  total (22/30)

TTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCT  >  NC_000913/3296145‑3296437                                                                                                                                                    |                                                                                                                                                  ttGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGc                                                                                                                                                 >  1:814750/1‑150 (MQ=255)   gTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCAt                                                                                                                                               <  1:601344/150‑1 (MQ=255)   gTTAGCCATTCCCCCCCGTGGCACCAGCCGGGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCa                                                                                                                                                <  1:188545/149‑1 (MQ=255)    ttAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCAtt                                                                                                                                              <  1:393257/150‑1 (MQ=255)         cATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCg                                                                                                                                        >  1:179134/1‑151 (MQ=255)             tCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATc                                                                                                                                    >  2:128386/1‑151 (MQ=255)              ccccccGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCAt                                                                                                                                     <  1:626945/146‑1 (MQ=255)               cACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGt                                                                                                                                  >  2:299273/1‑151 (MQ=255)                  ccGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAg                                                                                                                                >  2:629230/1‑150 (MQ=255)                    gTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGc                                                                                                                             >  1:206147/1‑151 (MQ=255)                    gTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGc                                                                                                                             >  1:376305/1‑151 (MQ=255)                         aaCAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGc                                                                                                                         >  1:609297/1‑150 (MQ=255)                                    tatTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATg                                                                                                             >  2:201611/1‑151 (MQ=255)                                      ttGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGc                                                                                                            <  2:831323/150‑1 (MQ=255)                                           aaaGCATTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATg                                                                                                       <  1:173678/150‑1 (MQ=255)                                           aaaGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATg                                                                                                       <  1:127199/150‑1 (MQ=255)                                            aaGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCt                                                                                                     <  1:348670/151‑1 (MQ=255)                                            aaGCAGTTAAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCat                                                                                                          <  1:450282/146‑1 (MQ=255)                                                         cgcggTTACGCGTGATATGACCATTGTGGCATTGGCCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGc                                                                                        <  2:341359/151‑1 (MQ=255)                                                         cgccgTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGc                                                                                        <  1:643953/151‑1 (MQ=255)                                                          gccgTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGcc                                                                                       >  2:366632/1‑151 (MQ=255)                                                           ccgTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCt                                                                                      >  2:687216/1‑151 (MQ=255)                                                           ccgTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCt                                                                                      <  2:191975/151‑1 (MQ=255)                                                            cgTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCtg                                                                                     >  1:781773/1‑151 (MQ=255)                                                                aCGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCg                                                                                 <  1:240305/151‑1 (MQ=255)                                                                         atGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATc                                                                         <  2:236926/150‑1 (MQ=255)                                                                         atGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATc                                                                         <  2:107396/150‑1 (MQ=255)                                                                         atGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCg                                                                        >  1:636319/1‑151 (MQ=255)                                                                         atGACCAGTGTGGTATTGACCGGCGATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACACGACGTGGTGATCTGCATTCCCTCGCATCGTAGTGCTCGCATTCACGAAATGCATATGCTGACGGTAAATTGCCTGTTCGATCTGaac                                                                         >  1:611458/1‑148 (MQ=255)                                                                                  gtgGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACg                                                                <  1:85743/150‑1 (MQ=255)                                                                                       aTTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTc                                                          >  2:76741/1‑151 (MQ=255)                                                                                             cGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCAc                                                    >  1:332823/1‑151 (MQ=255)                                                                                                   tGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGgat                                              >  1:142772/1‑151 (MQ=255)                                                                                                     aCGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGgatgt                                            <  1:848162/151‑2 (MQ=255)                                                                                                       ggcggcGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGAtt                                          <  1:652710/151‑1 (MQ=255)                                                                                                        gcggcgAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTTGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACTCGCTTTTCCCTCACCAGGAttatta                                         <  1:741412/151‑1 (MQ=255)                                                                                                            cgAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGga                                     <  1:723362/151‑1 (MQ=255)                                                                                                              aaCTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAgg                                      >  1:631466/1‑148 (MQ=255)                                                                                                              aaCTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAgg                                      >  2:753736/1‑148 (MQ=255)                                                                                                                   gCAGGTTTGTTAGGGCCGCAGGATGTTGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATAc                               <  1:86919/150‑1 (MQ=255)                                                                                                                   gCAGGTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACa                              <  2:631506/151‑1 (MQ=255)                                                                                                                      ggTTTGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATg                            >  1:830724/1‑150 (MQ=255)                                                                                                                         ttGTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAg                         <  1:657604/150‑1 (MQ=255)                                                                                                                           gTTAGGGCCACAGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCa                      <  1:264274/151‑1 (MQ=255)                                                                                                                                   cacaGGATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGcc              >  2:339970/1‑151 (MQ=255)                                                                                                                                       ggATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATc          <  1:422236/151‑1 (MQ=255)                                                                                                                                        gATGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCg         <  2:581711/151‑1 (MQ=255)                                                                                                                                         aTGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATc          >  1:339022/1‑149 (MQ=255)                                                                                                                                         aTGTGGAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGc        <  1:596439/151‑1 (MQ=255)                                                                                                                                             ggAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTc    <  2:206170/151‑1 (MQ=255)                                                                                                                                             ggAGATCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTc    <  2:390870/151‑1 (MQ=255)                                                                                                                                               agaTCTGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCt  <  1:164774/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCT  >  NC_000913/3296145‑3296437

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