Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	895,216	G→A	11.2%	intergenic (+56/+8)	PP_0777 → / ← PP_0778	glutathione peroxidase/group 1 family glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	895,216	0	G	A	11.2%	113.8 / 5.0	46	intergenic (+56/+8)	PP_0777/PP_0778	glutathione peroxidase/group 1 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (17/23);  new base A (3/2);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.08e-01

CAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAGGCGCA  >  NC_002947/895070‑895363                                                                                                                                                    |                                                                                                                                                    cAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCt                                                                                                                                                   <  3:175893/149‑1 (MQ=255) cAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCt                                                                                                                                                   <  5:171041/149‑1 (MQ=255) cAAGGATGGCCGGGTGCGGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCt                                                                                                                                                   <  3:235173/149‑1 (MQ=255)  aaGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTg                                                                                                                                                  >  1:69051/1‑149 (MQ=255)   aGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGGGGTGCAGGCCATCGAGAAAGCCCGGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTgg                                                                                                                                                 <  2:58613/149‑1 (MQ=255)      aTGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGAtc                                                                                                                                              <  3:167631/149‑1 (MQ=39)      aTGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGCATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGAtc                                                                                                                                              <  3:243808/149‑1 (MQ=255)       tGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGAtct                                                                                                                                             <  5:76153/149‑1 (MQ=255)         gCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTca                                                                                                                                           <  6:196202/149‑1 (MQ=255)                  aggCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTg                                                                                                                                  >  6:358876/2‑149 (MQ=255)                   ggCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGc                                                                                                                                 <  1:38094/149‑1 (MQ=255)                   ggCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTCACAGgtgcca                                                                                                                                  >  2:244194/1‑146 (MQ=37)                    gcgcGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTg                                                                                                                                  <  2:69051/147‑1 (MQ=255)                     cgcgTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCCCTGGATCTCACAGGCGCTGCaa                                                                                                                               <  7:54804/149‑1 (MQ=255)                              cgccgcGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgcg                                                                                                                       >  6:263301/1‑148 (MQ=255)                               gccgcGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTGCAGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGAtc                                                                                                                                              <  6:123684/125‑1 (MQ=25)                               gccgcGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTGCAGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGAtc                                                                                                                                              >  5:123684/1‑125 (MQ=25)                                        ggCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCa                                                                                                            <  5:114933/149‑1 (MQ=255)                                               gATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAg                                                                                                     >  2:127252/1‑149 (MQ=255)                                                     cctgCGGTGGTGCAGGCCATCGATAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCGGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAg                                                                                               <  8:72952/146‑1 (MQ=255)                                                     ccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAg                                                                                               >  8:93506/1‑149 (MQ=255)                                                          ggtggtGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGc                                                                                          <  1:235716/149‑1 (MQ=255)                                                             ggtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCgg                                                                                       >  7:208100/1‑149 (MQ=255)                                                                    gCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATgcgc                                                                                >  7:256323/1‑149 (MQ=255)                                                                      cATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTg                                                                              >  3:19647/1‑149 (MQ=255)                                                                       aTCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGt                                                                             >  1:217695/1‑149 (MQ=255)                                                                       aTCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGt                                                                             <  3:23/149‑1 (MQ=255)                                                                        tCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGtt                                                                            <  8:208100/149‑1 (MQ=255)                                                                          gagaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTc                                                                           <  5:358876/148‑1 (MQ=255)                                                                               aGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAc                                                                     <  7:93506/149‑1 (MQ=255)                                                                                gCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAcc                                                                    >  3:106479/1‑149 (MQ=255)                                                                                 cccTTGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCCGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAccc                                                                   <  1:234357/149‑1 (MQ=255)                                                                                 cccTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAccc                                                                   <  4:19647/149‑1 (MQ=255)                                                                                    tGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTc                                                                >  6:261021/1‑149 (MQ=255)                                                                                         tGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACt                                                           >  8:139044/1‑149 (MQ=255)                                                                                           aaGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGc                                                                                                                                 <  4:109818/77‑1 (MQ=255)                                                                                           aaGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGc                                                                                                                                 >  3:109818/1‑77 (MQ=255)                                                                                               gTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCgt                                                     <  2:217695/149‑2 (MQ=255)                                                                                               gTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGa                                                     <  1:127252/149‑1 (MQ=255)                                                                                                              ttt‑tGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGc                                      >  4:349227/1‑149 (MQ=255)                                                                                                                    gCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTgg                                 <  8:256323/149‑1 (MQ=255)                                                                                                                        cTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTTCAAACATTGCGCAGAACGCCCAGACGTAGCACCAg                                                                                               >  7:224570/1‑83 (MQ=255)                                                                                                                        cTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggcggc                             >  1:25709/1‑149 (MQ=255)                                                                                                                        cTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggcggc                             <  5:263301/149‑1 (MQ=255)                                                                                                                         tgtgAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCt                            >  4:39128/1‑149 (MQ=255)                                                                                                                                tccagcgccgccCGCGCGGCACTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGc                     >  6:117581/7‑149 (MQ=39)                                                                                                                                              cgcggcaCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAg       <  5:117581/148‑1 (MQ=255)                                                                                                                                               gcggcgCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAgg      <  3:40090/149‑1 (MQ=255)                                                                                                                                                 ggcgCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAGgc     >  2:6313/1‑148 (MQ=255)                                                                                                                                                 ggcgCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAGgc     <  1:6313/148‑1 (MQ=255)                                                                                                                                                   cgcTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAGGCGCa  <  4:106479/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAGGCGCA  >  NC_002947/895070‑895363

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