Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,123,959	Δ1 bp	coding (273/405 nt)	glcG ←	putative heme‑binding protein GlcG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,123,959	0	T	.	100.0%	149.0 / NA	35	coding (273/405 nt)	glcG	putative heme‑binding protein GlcG
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (18/17);  total (18/17)

TTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGC  >  NC_000913/3123827‑3124081                                                                                                                                     |                                                                                                                           ttATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGCAACGCCCACGGCCCCAATAATTTGCCCATCCACCACAACCGGTCTTCCGCCTTCCAGCGCAGTTAA‑AACGGCGCAGt                                                                                                                 >  1:718258/1‑143 (MQ=255) ttATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGt                                                                                                                 >  2:426316/1‑143 (MQ=255)  tATTTCGCTAATCCCGCTGCTGCCGCGTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTc                                                                                                                <  1:307086/143‑1 (MQ=255)   aTTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATCATTTGCCCATCCACGACAACCGGTGCGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGtaa                                                                                                               >  1:87286/1‑141 (MQ=255)       cGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCCCCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCCTCCAGCGACGTTAA‑AACGGCGCAGTCACGaa                                                                                                           >  1:518752/1‑143 (MQ=255)          cAACACCGCTGCTGCCGCTTGCGCGACCAGCGCATCCGGTGCGACGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAAcgc                                                                                                        <  1:286687/143‑1 (MQ=255)              acCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACTGCAACAATAATTTGCAAAAACACAACAACCGGTACGCCGCATTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTa                                                                                                    >  1:61217/1‑143 (MQ=255)                     ctgcCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATACTTTTCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGATACGTCCgt                                                                                             >  1:54081/1‑143 (MQ=255)                       gcCGCTGTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCgttg                                                                                           <  2:840360/143‑1 (MQ=255)                                gcgACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AAAGGAGCAGTCAAGAACGCGGTACGTCCGTTGTTCACCAtc                                                                                  >  1:566549/1‑143 (MQ=255)                                     cTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACTCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTc                                                                             <  1:833047/143‑1 (MQ=255)                                               tgtgCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTa                                                                   <  1:370002/143‑1 (MQ=255)                                                     ccGGTTAAACCAGGAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCa                                                             <  1:99479/143‑1 (MQ=255)                                                     ccGGTTAAACCAGAAACGCACACGGCACCAATAATTTGCCCAGCCACAACAACCGGTACGCCGGTTTGGTCCTACTTTAA‑GTTGGCGCAATCACAAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCa                                                             <  1:620910/143‑1 (MQ=255)                                                          tAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCacgacg                                                        >  1:608284/1‑143 (MQ=255)                                                              ccAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCCTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACCCGGTACGTCCCTTGTTAACCATCTCTTCATAGCCCTTAGTTTCACGAAGcccc                                                    >  1:179283/1‑143 (MQ=255)                                                               cAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGCGTCACGACGACCCa                                                   >  2:746681/1‑143 (MQ=255)                                                               cAGAAACGCCCACGGAACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGGCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCa                                                   <  1:763609/143‑1 (MQ=255)                                                                      gCCCACGCCATCAATACTTTTCCCTTCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGTGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCgcggc                                            <  2:130011/143‑1 (MQ=255)                                                                         cACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGt                                         <  1:426316/143‑1 (MQ=255)                                                                          aCGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTGACCATCTCTTCATAGCCCAAAGTTGCACGACGCCCCAGCACGGCGTTa                                        >  2:548951/1‑143 (MQ=255)                                                                               aCCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAATTTCACAACGCACCAACGCGGCGGACcgcga                                   >  1:838660/1‑142 (MQ=255)                                                                                       tttGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTTGCCGCCCGACGCCCCAGCGCGGCGGTACGCCCTTTCTCCt                           >  1:573801/1‑143 (MQ=255)                                                                                        cggCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTg                          <  1:143084/141‑1 (MQ=255)                                                                                         tGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCGTCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTgg                         <  1:600967/143‑1 (MQ=255)                                                                                            ccATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGgag                      <  2:103545/143‑1 (MQ=255)                                                                                              aTCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGAGCAGTCACGAACGCGGTACGTCCGTTGTTCACAATCTCTTCATAGCCCTTAGTTTCACGCCGCCCCAGCGCGGCGGTACGCGCTTTATCCTGCGAGat                    >  2:274219/1‑143 (MQ=255)                                                                                                  acaacaACCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATaa                <  2:608284/143‑1 (MQ=255)                                                                                                    aacaacCGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTGCGACCGTAGTTCCCCATCACTTCATAGCCCTTAGCTTCCCGACGCCCCCGCGCGGCGGCACGCGCTTTCTCCTCGGATATCTaacc              >  2:240927/1‑141 (MQ=255)                                                                                                         ccGGTACGCCGCCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCAACTCTCCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCATTCTCCTGGGAGATATAAGACGCac         >  2:249767/1‑142 (MQ=255)                                                                                                            gtgcgccgcCTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCTCCAGGGCTGCGTTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCg      <  2:179283/140‑1 (MQ=255)                                                                                                            ggaCGCCGTCTTCCAGCGCCGTTAA‑ACCGGCGCAGTCACGCATGCGGTACGTCCGTTGTTCACCATCTCTTCATAGTCCTTAGTTTCACGACGCTCCAGCGCGGCGGTACGCGTTTTCTCCTGGGAGATATAAGCCGCAATCg      <  2:92643/141‑1 (MQ=255)                                                                                                             gactcctccTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGTGCTTTCTCCTGGGAGATATAAGCCGCAATCgg     <  2:87286/136‑1 (MQ=255)                                                                                                                gccgccTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCGCTTTCTCTTCATAGCCCTTAGTTTCACGACGCCCTAACGAGGAGATACGCGCTTTCTACTGGGAGATATAAGCAGAAATCGgcga  >  1:40369/1‑142 (MQ=255)                                                                                                                gccgccTTCCAGCGACGTTAA‑AACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGACATAACCCGCAATCGgcgc  >  2:186121/1‑143 (MQ=255)                                                                                                                                     |                                                                                                                           TTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGGCGCAGTCACGAACGCGGTACGTCCGTTGTTCACCATCTCTTCATAGCCCTTAGTTTCACGACGCCCCAGCGCGGCGGTACGCGCTTTCTCCTGGGAGATATAAGCCGCAATCGGCGC  >  NC_000913/3123827‑3124081

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