Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	3,787,400	Δ1 bp	coding (1235/1392 nt)	PP_3350 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,787,400	0	T	.	100.0%	82.1 / NA	19	coding (1235/1392 nt)	PP_3350	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (8/11);  total (8/11)

GGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAG  >  NC_002947/3787258‑3787522                                                                                                                                               |                                                                                                                           ggCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGAc                                                                                                                      >  7:1666/1‑148 (MQ=255)              aCAAACGACAAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAATGTACAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGgt                                                                                                         <  4:65177/148‑1 (MQ=255)                        aaTCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGGCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGc                                                                                               <  3:42234/148‑1 (MQ=255)                          tCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCAt                                                                                             <  1:16353/148‑1 (MQ=255)                                            ttGCTGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTa                                                                           <  7:137510/148‑1 (MQ=255)                                            ttGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGt                                                                            <  2:52036/147‑1 (MQ=255)                                            ttGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGt                                                                            >  4:104013/1‑147 (MQ=255)                                                      tGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCt                                                                 <  4:64043/148‑1 (MQ=255)                                                          cACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTcc                                                              >  7:129330/1‑147 (MQ=255)                                                            cgtcgtAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCa                                                           <  2:183905/148‑1 (MQ=255)                                                                            gTCAAGGTCCAGGGTCAGGTTGGGAGTGCCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACtt                                           >  8:103868/1‑148 (MQ=255)                                                                                aGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTc                                          >  5:52421/1‑145 (MQ=255)                                                                                       gggTCAGGGTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCATGCTGTTTCCACACTTCGAAGATGCCg                                <  3:53610/148‑1 (MQ=255)                                                                                                    aGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCa                   <  3:48078/148‑1 (MQ=255)                                                                                                            gCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTgcg            >  6:70348/1‑147 (MQ=255)                                                                                                                    gTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAggg      >  4:60467/1‑145 (MQ=255)                                                                                                                    gTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAggg      >  6:96445/1‑145 (MQ=255)                                                                                                                    gTGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCga   <  8:1666/148‑1 (MQ=255)                                                                                                                     tGCTGGCGCGGTAGATGGTACCGACA‑GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCgag  <  3:163931/148‑1 (MQ=255)                                                                                                                                               |                                                                                                                           GGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAG  >  NC_002947/3787258‑3787522

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