Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	3,787,503	Δ2 bp	coding (1131‑1132/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,501	0	G	.	100.0%	96.2 / NA	22	coding (1134/1392 nt)	PP_3350	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (12/10);  total (12/10)
*	NC_002947	3,787,502	0	C	.	100.0%	96.2 / NA	22	coding (1133/1392 nt)	PP_3350	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (12/10);  total (12/10)

AAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCC  >  NC_002947/3787363‑3787629                                                                                                                                           ||                                                                                                                                gcgCCAGCACGCGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTgcg                                                                                                                       >  7:12238/3‑148 (MQ=255)                              tACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATg                                                                                                                                  <  4:37865/110‑1 (MQ=255)                                                tGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCa                                                                        <  1:124089/148‑1 (MQ=255)                                                            tGGCCAGCATGCTCATGCCGGGCATGTACACGCCGGCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAggc                                                            <  8:52680/148‑1 (MQ=255)                                                              gCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGGggct                                                           <  4:168090/147‑1 (MQ=255)                                                              gCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAggct                                                           <  3:39703/147‑1 (MQ=255)                                                              gCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTGCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAggctg                                                          >  4:123434/1‑148 (MQ=255)                                                               ccAGCATGCTCATTCCGGGCTTGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGGCGAAGGGCGGCCCCAGCTCGTTGATGTTGCTCAGGGTGGTGAggctgg                                                         <  8:163187/148‑1 (MQ=255)                                                                cAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAggctggc                                                        >  6:109948/1‑148 (MQ=255)                                                                        tCATGCCGGGCAGGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGAGGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCa                                                 <  3:161276/147‑1 (MQ=255)                                                                           tGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATa                                               <  4:20263/146‑1 (MQ=255)                                                                               gggCATGTACACACCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACAcc                                         >  3:60745/1‑148 (MQ=255)                                                                                         acGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCCCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTgg                                >  6:152940/1‑147 (MQ=255)                                                                                          cGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGt                              <  8:9001/148‑1 (MQ=255)                                                                                           gCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGt                              >  2:128576/1‑147 (MQ=255)                                                                                                 ttCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCgg                       >  7:41881/1‑148 (MQ=255)                                                                                                    cTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGt                     >  4:29125/1‑147 (MQ=255)                                                                                                     tCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTc                    >  6:27553/1‑147 (MQ=255)                                                                                                       aCGCTGTTTCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGaa                 >  1:192257/1‑148 (MQ=255)                                                                                                          cTGTATCCACACTTCGAAGAGGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGGTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGt              <  8:149781/148‑1 (MQ=255)                                                                                                             tttCCACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGc            >  8:119344/1‑147 (MQ=255)                                                                                                                ccACACTTCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCTCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCgcc        >  3:41681/1‑146 (MQ=255)                                                                                                                      ttCGAAGATGCCGGGCTCGAT‑‑GCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTcc  <  7:119344/148‑1 (MQ=255)                                                                                                                                           ||                                                                                                                                AAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCC  >  NC_002947/3787363‑3787629

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