Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	3,788,529	Δ1 bp	coding (106/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,788,526	0	C	.	100.0%	100.5 / NA	23	coding (109/1392 nt)	PP_3350	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (7/16);  total (7/16)

CCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTG  >  NC_002947/3788382‑3788665                                                                                                                                                 |                                                                                                                                            ccGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCCTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGccct                                                                                                                                          <  3:96541/148‑2 (MQ=255)   gAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGccctgg                                                                                                                                        >  2:58859/1‑145 (MQ=255)              tgcaactcGGGGCTGAGGCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCccggcc                                                                                                                            >  1:136505/8‑148 (MQ=255)                     cGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGctctc                                                                                                                     <  5:62442/148‑1 (MQ=255)                              gCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGc                                                                                                              <  5:42174/146‑1 (MQ=255)                              gCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCg                                                                                                             <  1:153125/147‑1 (MQ=255)                              gCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCg                                                                                                             <  1:200482/147‑1 (MQ=255)                                           gCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAg                                                                                               <  3:155066/148‑1 (MQ=255)                                              aGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCg                                                                                            <  3:16064/148‑1 (MQ=255)                                                  cGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAAt                                                                                        <  2:59455/148‑1 (MQ=255)                                                   gACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATg                                                                                       <  1:172225/148‑1 (MQ=255)                                                          aCAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAATGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCAGTCTGCAAGAGGGCTGGCCTCGTAGGCGGACTGGCCCAg                                                                                 >  8:67859/1‑147 (MQ=255)                                                           cAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGcaca                                                                                  <  4:52693/145‑1 (MQ=255)                                                                gATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGc                                                                          <  8:49956/148‑1 (MQ=255)                                                                               cATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGggcag                                                            >  4:62839/1‑147 (MQ=255)                                                                                 tCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGggcaggc                                                          >  6:160964/1‑147 (MQ=255)                                                                                 tCGAAAGGGGCAGTGGCCTTGGCGGGGGCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGggcaggca                                                         <  1:63521/148‑1 (MQ=255)                                                                                                            tcgtcgAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGAtcatca                               >  5:189583/1‑147 (MQ=255)                                                                                                               tcgAGAAAGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAAt                           >  5:20068/1‑148 (MQ=255)                                                                                                                    aaaGCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTccc                      <  6:189583/148‑1 (MQ=255)                                                                                                                       gCGGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTc                   <  8:22902/148‑1 (MQ=255)                                                                                                                              tCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAACGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGtt            >  7:40546/1‑148 (MQ=255)                                                                                                                                        aGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAAtt   <  4:27291/147‑1 (MQ=255)                                                                                                                                        aGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCGGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTg  <  8:41587/148‑1 (MQ=255)                                                                                                                                        aGCGAGGGG‑CCCTGGCCGCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCGTGATCATCAGAATTTCCCTTCACTTGTTTTTATAAtt   <  5:94804/147‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                            CCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTG  >  NC_002947/3788382‑3788665

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