Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,787,536:1	+C	100%	coding (1099/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,536	1	.	C	100.0%	183.7 / NA	54	P367A (CCG→GCG)	PP_3350	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (31/23);  total (31/23)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGG‑CCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGC  >  NC_002947/3787401‑3787674                                                                                                                                         |                                                                                                                                           ggCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACgat                                                                                                                                >  6:103375/1‑149 (MQ=255)                   tgtTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGggtggt                                                                                                              <  5:196180/149‑1 (MQ=255)                        gCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGTTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAggc                                                                                                         >  7:299453/1‑149 (MQ=255)                        gCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAgg                                                                                                          <  6:239911/148‑1 (MQ=255)                        gCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAgg                                                                                                          >  5:239911/1‑148 (MQ=255)                                gCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCg                                                                                                 >  3:131636/1‑149 (MQ=255)                                 cTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGc                                                                                                <  7:86372/149‑1 (MQ=255)                                 cTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGc                                                                                                <  4:131636/149‑1 (MQ=255)                                 cTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGc                                                                                                <  4:130291/149‑1 (MQ=255)                                 cTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGc                                                                                                <  4:118185/149‑1 (MQ=255)                                  tCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGc                                                                                                >  1:43808/1‑148 (MQ=255)                                   cATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGG‑ccccc                                                                                                                                       >  3:110047/1‑106 (MQ=255)                                   cATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGG‑ccccc                                                                                                                                       <  4:110047/107‑2 (MQ=255)                                    aTGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCCGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCAt                                                                                             <  8:252345/149‑1 (MQ=255)                                    aTGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAggctg                                                                                                       >  3:160075/1‑139 (MQ=255)                                    aTGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAggctg                                                                                                       <  4:160075/139‑1 (MQ=255)                                         gggCATGTACACGCCGTCTTCCTCACGCTGTTTCGACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGaca                                                                                        >  7:138847/1‑149 (MQ=255)                                         gggCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAggct                                                                                                        <  8:27479/133‑1 (MQ=255)                                         gggCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAggct                                                                                                        >  7:27479/1‑133 (MQ=255)                                            cATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCa                                                                                     <  2:43808/149‑1 (MQ=255)                                                             ccTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCg                                                                     <  8:138847/148‑1 (MQ=255)                                                             ccTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGATTCGCCATAGCAACAATTGCGTGGGTACAGctg                                                                    >  4:253328/1‑147 (MQ=255)                                                               tCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTgg                                                                             >  2:299011/1‑138 (MQ=255)                                                               tCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTgg                                                                             <  1:299011/138‑1 (MQ=255)                                                                      gTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATGGCGGGGGTAAAGCGGGTCGAAGGt                                                           >  4:256411/1‑149 (MQ=255)                                                                        ttCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTgg                                                                             <  2:260289/129‑1 (MQ=255)                                                                        ttCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGATCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTgg                                                                             >  1:260289/1‑129 (MQ=255)                                                                           cacaCTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGt                                                                                                                           <  1:285917/80‑1 (MQ=255)                                                                           cacaCTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGt                                                                                                                           >  2:285917/1‑80 (MQ=255)                                                                           cacaCTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGa                                                      <  5:103375/149‑1 (MQ=255)                                                                           cacaCTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGa                                                      >  7:437/1‑149 (MQ=255)                                                                            acacTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGTCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGAc                                                     >  8:260413/1‑149 (MQ=255)                                                                            acacTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGAc                                                     >  7:82853/1‑149 (MQ=255)                                                                            acacTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGAGGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGAc                                                     >  6:81125/1‑149 (MQ=255)                                                                              acTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACtt                                                   >  6:143875/1‑149 (MQ=255)                                                                                ttggAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGaa                                                              >  8:264940/4‑136 (MQ=255)                                                                                ttggAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGaa                                                              <  7:264940/133‑1 (MQ=255)                                                                                ttCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGGAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACtttt                                                 >  4:230573/1‑149 (MQ=255)                                                                                ttCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGCAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGACGGTGCCGACgtt                                                  >  4:189222/1‑148 (MQ=255)                                                                                ttCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACtttt                                                 >  2:321265/1‑149 (MQ=255)                                                                                ttCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACtttt                                                 >  2:196114/1‑149 (MQ=255)                                                                                ttCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACtttt                                                 >  7:6798/1‑149 (MQ=255)                                                                                   gAAGATGACGGGCTCGATGCTCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGTGTAAAGCGGGTCGAAGGTGCCGACTTTTcc                                               <  8:299453/148‑1 (MQ=255)                                                                                     aGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTc                                            >  2:112066/1‑149 (MQ=255)                                                                                                tCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCAc                                 >  6:295338/1‑149 (MQ=255)                                                                                                tCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCAc                                 >  6:267360/1‑149 (MQ=255)                                                                                                              gcgcCAGGGCGAGAAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCg                   <  8:82853/149‑1 (MQ=255)                                                                                                                gcCAGGGCGAGGAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGc                 <  3:230573/149‑1 (MQ=255)                                                                                                                  cAGGGCGAGAAGTCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTGCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCAt               <  7:260413/149‑1 (MQ=255)                                                                                                                        gagaAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATgc         >  7:96016/1‑149 (MQ=255)                                                                                                                        gagaAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATgc         >  2:122936/1‑149 (MQ=255)                                                                                                                        gagaAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATgc         <  1:196114/149‑1 (MQ=255)                                                                                                                        gagaAGCCGAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATgc         <  1:112066/149‑1 (MQ=255)                                                                                                                          gaAGCCGAAGGTCGGCCCCCACCATGATGATGTTTCTCTGGGTGGTGAGGCTGGATTCGCCATAGACACCATTGCGGGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGAGGCTTGCGGCATGGATgcgc       <  8:108084/149‑1 (MQ=255)                                                                                                                                gAAGGTCGGCCCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGc  >  5:273020/1‑148 (MQ=255)                                                                                                                                gAAGGTCGGCCCCCACCACGATGATGTTCCTCAGGGTCGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGACGAAGGTGCCGCCTTTTCCGTCGCCCAGGTCACTGTCGCCGATT‑CGGCAGCGATGCGCATGGc  >  5:257890/1‑147 (MQ=255)                                                                                                                                         |                                                                                                                                           GGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGCCGGGCTCGATGCGCAGCGTGCGCCAGGGCGAGAAGCCGAAGGTCGG‑CCCCACCACGATGATGTTGCTCAGGGTGGTGAGGCTGGCTTCGCCATAGACACCATTGCGTGGGTAAAGCGGGTCGAAGGTGCCGACTTTTCCGTCGCCCAGGTCACTGTCGCCGCTTGCGGCATCGATGCGCATGGC  >  NC_002947/3787401‑3787674

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