Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,904,996	A→C	7.5%	intergenic (‑163/+72)	PP_4312 ← / ← PP_4313	hypothetical protein/peptidylprolyl isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,904,996	0	A	C	7.5%	164.0 / 4.8	53	intergenic (‑163/+72)	PP_4312/PP_4313	hypothetical protein/peptidylprolyl isomerase
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base C (2/2);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

GAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCC  >  NC_002947/4904854‑4905139                                                                                                                                                |                                                                                                                                                gagaATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGcc                                                                                                                                            >  3:114448/1‑148 (MQ=255)        tCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCTTTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCCTTGCCTAgggaac                                                                                                                                    <  1:263040/149‑5 (MQ=255)        tCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAg                                                                                                                                    >  7:187093/1‑149 (MQ=255)          gggTTTACCCGCCACATCGCCCTGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGATACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAgg                                                                                                                                   <  1:86775/149‑1 (MQ=255)          gggTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCCTTGCCTAgggaacg                                                                                                                                   <  7:56105/148‑6 (MQ=255)           ggTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAgggg                                                                                                                                 <  6:97052/149‑1 (MQ=255)                aCCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTgg                                                                                                                             <  3:191797/148‑1 (MQ=255)                 cccGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGgt                                                                                                                           >  2:252405/1‑149 (MQ=255)                    gCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGtt                                                                                                                        >  3:4523/1‑149 (MQ=255)                    gCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGtt                                                                                                                        >  4:202172/1‑149 (MQ=255)                     ccACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTg                                                                                                                       <  8:187093/149‑1 (MQ=255)                     cacgggcGCCCTGCAACTAGCCGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGc                                                                                                                                             >  8:220192/7‑127 (MQ=255)                     cacgggcGCCCTGCAACTAGCCGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGc                                                                                                                                             <  7:220192/121‑1 (MQ=255)                         aTCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTgcgc                                                                                                                    >  5:294824/1‑148 (MQ=255)                            gCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAAtt                                                                                                                <  5:5045/149‑1 (MQ=255)                            gCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAAtt                                                                                                                <  4:114448/149‑1 (MQ=255)                            gCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAAtt                                                                                                                >  6:5045/1‑149 (MQ=255)                                     gTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAg                                                                                                         <  3:202172/147‑1 (MQ=255)                                      ttGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTgtgt                                                                                                            >  4:81587/1‑143 (MQ=255)                                      ttGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTgtgt                                                                                                            <  3:81587/143‑1 (MQ=255)                                      ttGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTg                                                                                                      <  1:126544/149‑1 (MQ=255)                                        gCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGc                                                                                                     >  1:106671/1‑148 (MQ=255)                                         cAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGTGGCTGGGTGTTGCGCAATTGTGTAAGGTGCaa                                                                                                   <  1:36996/149‑1 (MQ=255)                                                                cATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGaa                                                                             >  5:168119/1‑148 (MQ=255)                                                                 aTCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAaca                                                                           <  6:294824/149‑1 (MQ=255)                                                                    cTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTc                                                                        <  1:252405/149‑1 (MQ=255)                                                                     ttGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTc                                                                        >  1:11412/1‑148 (MQ=255)                                                                       ggCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTc                                                                        <  2:11412/146‑1 (MQ=255)                                                                       ggCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCaa                                                                      <  6:303237/148‑1 (MQ=255)                                                                           aCGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAatcatc                                                                 <  1:306152/149‑1 (MQ=255)                                                                             gcgTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAg                                                                                                  <  8:221187/114‑1 (MQ=255)                                                                             gcgTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAg                                                                                                  >  7:221187/1‑114 (MQ=255)                                                                                 gAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGt                                                           >  7:131866/1‑149 (MQ=255)                                                                                  aGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGt                                                           <  8:131866/148‑1 (MQ=255)                                                                                   gggTCAGAGTGCCCCGA‑TCCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGa                                                         <  4:4523/149‑1 (MQ=255)                                                                                    ggTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAg                                                        >  1:153422/1‑149 (MQ=255)                                                                                           gTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTc                                                  >  3:77959/1‑148 (MQ=255)                                                                                            tGCCCCGA‑ACCCGTTTGGTTCCTCTCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATccc                                                               <  1:105744/134‑1 (MQ=255)                                                                                            tGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATccc                                                               >  2:105744/1‑134 (MQ=255)                                                                                              ccccGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCcag                                               <  2:106671/148‑1 (MQ=255)                                                                                                       ccGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGa                                       >  3:290253/1‑148 (MQ=255)                                                                                                              gTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAag                               <  4:290253/149‑1 (MQ=255)                                                                                                                tcctccCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAaggt                             >  6:308239/1‑149 (MQ=255)                                                                                                                  ctccCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAg                                                                                                         >  2:108962/1‑71 (MQ=255)                                                                                                                  ctccCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAg                                                                                                         <  1:108962/71‑1 (MQ=255)                                                                                                                   tccCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAaggtcag                          >  2:91386/1‑149 (MQ=255)                                                                                                                     ccTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAaggtcag                          <  1:91386/147‑1 (MQ=255)                                                                                                                               aaTGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACGCTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTcgg              >  6:13044/1‑149 (MQ=255)                                                                                                                                aTGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTcggc             <  4:77959/149‑1 (MQ=255)                                                                                                                                   ccGTTCCCTAGGCAAGGCCTAGGTTAGGGGCGGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTcggcgg           >  1:194404/1‑148 (MQ=255)                                                                                                                                       aaccTAGGCCTTGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCg      >  3:248014/3‑149 (MQ=255)                                                                                                                                        accTAGGCCGTGCCTAGGTTAGGGGCGGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGCACACTCAATCCTCCCGCGGGCGCGCTTCCAGCTGTTCGGGCGCGACGGTCCGGGCCGAGGTCGGGGGGCTCGc     >  4:61482/2‑149 (MQ=255)                                                                                                                                           tAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGcccc  >  5:146912/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                GAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGA‑ACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCC  >  NC_002947/4904854‑4905139

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