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%	137.4 / NA	44	*234* (TAG→TAA)	hda	ATPase regulatory factor involved in DnaA inactivation
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/21);  total (23/21)

CTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAG  >  NC_000913/2617952‑2618201                                                                                                                            |                                                                                                                                  cTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAAtt                                                                                                                    >  1:225053/1‑139 (MQ=255)      tGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAAtttcttt                                                                                                               <  1:709736/139‑1 (MQ=255)                        gCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                 <  2:93642/135‑1 (MQ=255)                        gCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                 >  1:93642/1‑135 (MQ=255)                           aCGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCtt                                                                                          >  1:328439/1‑139 (MQ=255)                                           aTCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATc                                                                          >  1:407749/1‑139 (MQ=255)                                                gCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGc                                                                     >  2:266873/1‑139 (MQ=255)                                                  cTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCAc                                                                   >  2:402982/1‑139 (MQ=255)                                                     cGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGAt                                                                >  2:529434/1‑139 (MQ=255)                                                        aTACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCa                                                             <  2:225053/139‑1 (MQ=255)                                                         tACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGCTCGCGCCGCATCCGACAATAAACGCCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCaa                                                            <  2:502255/139‑1 (MQ=255)                                                          aCTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGGCAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAAc                                                           <  2:569146/139‑1 (MQ=255)                                                             aTGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGa                                                        >  2:170177/1‑139 (MQ=255)                                                              tGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGAt                                                       >  1:663601/1‑139 (MQ=255)                                                               ggATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                    >  1:172965/1‑125 (MQ=255)                                                               ggATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                                                    <  2:172965/125‑1 (MQ=255)                                                               ggATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATc                                                      >  2:568598/1‑139 (MQ=255)                                                                 aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACGAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                    >  1:408912/1‑139 (MQ=255)                                                                 aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                    <  2:118955/139‑1 (MQ=255)                                                                 aTTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCa                                                    <  1:295972/139‑1 (MQ=255)                                                                   tAACGTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑c                                               <  1:73282/139‑1 (MQ=255)                                                                    aaCTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑Cg                                              >  1:674436/1‑139 (MQ=255)                                                                          gcgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATaa                                        >  2:763771/1‑139 (MQ=255)                                                                          gcgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATaa                                        >  1:364633/1‑139 (MQ=255)                                                                           cgTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATaaa                                       >  2:176593/1‑139 (MQ=255)                                                                               gTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGc                                   <  2:408912/139‑1 (MQ=255)                                                                               gTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGc                                   <  2:364633/139‑1 (MQ=255)                                                                               gTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGc                                   >  1:64996/1‑139 (MQ=255)                                                                                    gATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGCGTgcg                              <  1:630572/139‑1 (MQ=255)                                                                                    gATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGCGTgcg                              <  1:402982/139‑1 (MQ=255)                                                                                     aTAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGCGTgcgc                             <  1:72355/139‑1 (MQ=255)                                                                                          gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGCGt                                 <  2:131975/130‑1 (MQ=255)                                                                                          gCGTTCGCGCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGCGt                                 >  1:131975/1‑130 (MQ=255)                                                                                               cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAAt                                                                           <  2:711171/86‑1 (MQ=255)                                                                                               cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAAt                                                                           >  1:711171/1‑86 (MQ=255)                                                                                               cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCGTCATAAATAGCGTGCGCATTtctc                      >  2:123620/1‑139 (MQ=255)                                                                                               cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGCGTGCGCATTtct                       <  2:135891/135‑1 (MQ=255)                                                                                               cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGCGTGCGCATTtct                       >  1:135891/1‑135 (MQ=255)                                                                                               cgcgCCGCATCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGCGTGCGCATTTCTCTGt                   <  2:590384/139‑1 (MQ=255)                                                                                                        tCCGACAATAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCt          <  1:529434/139‑1 (MQ=255)                                                                                                                tAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAg  >  2:12806/1‑139 (MQ=255)                                                                                                                tAAACACCTTATTTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAg  <  1:12806/139‑1 (MQ=255)                                                                                                                         tATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                 >  2:174832/1‑38 (MQ=255)                                                                                                                         tATTTACAACTTCAGAATTTCTTTCACAAACGGAATgg                                                                                                 <  1:174832/38‑1 (MQ=255)                                                                                                                            |                                                                                                                                  CTAGTTGAAGCACATTGCCGGATGCGACGCTTGCCACGTCTTATCCGGCCTACGAATACTATGGATTAACTTCGCGTAGTTCGGATAAGGCGTTCGCGCCGCATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAA‑‑‑CGTCATAAATAGCGTGCGCATTTCTCTGTCGAGCCGCTTCAGCAAG  >  NC_000913/2617952‑2618201

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