Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,618,075	C→T	*234* (TAG→TAA)	hda ←	ATPase regulatory factor involved in DnaA inactivation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,618,075	0	C	T	100.0%	106.2 / NA	34	*234* (TAG→TAA)	hda	ATPase regulatory factor involved in DnaA inactivation
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/14);  total (20/14)

AGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACGCC  >  NC_000913/2617954‑2618208                                                                                                                          |                                                                                                                                      aGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAAtttc                                                                                                                      >  1:173893/1‑139 (MQ=255)  gTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAAtttct                                                                                                                     <  1:653336/139‑1 (MQ=255)      aaGCACATTGCCGGATGCGGCGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTca                                                                                                                 >  1:559034/1‑139 (MQ=255)       aGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTcac                                                                                                                >  1:140359/1‑139 (MQ=255)                ccGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAAt                                                                                                       >  1:255224/1‑139 (MQ=255)                  ggATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                     <  2:140359/139‑1 (MQ=255)                                                            tGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCTTCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGAt                                                           >  2:381297/1‑139 (MQ=255)                                                            tGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGAt                                                           >  1:204626/1‑139 (MQ=255)                                                               aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                        <  2:255224/139‑1 (MQ=255)                                                               aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                        >  1:675072/1‑139 (MQ=255)                                                                    cTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                        <  2:598023/118‑1 (MQ=255)                                                                    cTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                        >  1:598023/1‑118 (MQ=255)                                                                        gcgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATaa                                               >  2:431399/1‑139 (MQ=255)                                                                        gcgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATaa                                               >  2:434827/1‑139 (MQ=255)                                                                             gTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGc                                          >  1:82136/1‑139 (MQ=255)                                                                              ttCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCg                                         >  2:51001/1‑139 (MQ=255)                                                                                cGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTg                                       <  1:434827/139‑1 (MQ=255)                                                                                       ggCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCAttt                                <  1:51001/139‑1 (MQ=255)                                                                                        gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                     <  1:467134/69‑1 (MQ=255)                                                                                        gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                     >  2:467134/1‑69 (MQ=255)                                                                                        gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                     <  1:357249/69‑1 (MQ=255)                                                                                        gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                     >  2:357249/1‑69 (MQ=255)                                                                                        gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACg                                                                          <  2:724477/96‑1 (MQ=255)                                                                                        gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACg                                                                          >  1:724477/1‑96 (MQ=255)                                                                                        gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                        >  2:278577/1‑98 (MQ=255)                                                                                        gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                        <  1:278577/98‑1 (MQ=255)                                                                                        gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTtc                               >  1:259182/1‑139 (MQ=255)                                                                                         cGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTtct                              <  2:259182/139‑1 (MQ=255)                                                                                             cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGc                                                                         >  2:421851/1‑92 (MQ=255)                                                                                             cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGc                                                                         <  1:421851/92‑1 (MQ=255)                                                                                                    cATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCg                   <  1:815453/139‑1 (MQ=255)                                                                                                                    ccTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGt                                                                                  <  2:239414/60‑1 (MQ=255)                                                                                                                    ccTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGt                                                                                  >  1:239414/1‑60 (MQ=255)                                                                                                                     cTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAGCTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACGcc  >  2:599635/1‑139 (MQ=255)                                                                                                                          |                                                                                                                                      AGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAGAAACGCC  >  NC_000913/2617954‑2618208

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