Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	895,210	G→C	10.9%	intergenic (+50/+14)	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,210	0	G	C	10.9%	129.7 / 7.4	46	intergenic (+50/+14)	PP_0777/PP_0778	glutathione peroxidase/group 1 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (17/24);  new base C (3/2);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

TGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGC  >  NC_002947/895063‑895344                                                                                                                                                     |                                                                                                                                       tggtggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCgcgg                                                                                                                                        >  5:337039/1‑148 (MQ=255)   gtggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCgcgg                                                                                                                                        <  4:163587/146‑1 (MQ=255)   gtggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCgcgg                                                                                                                                        <  6:337039/146‑1 (MQ=255)     gggCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGggcggc                                                                                                                                   >  4:40090/1‑149 (MQ=255)        cAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCt                                                                                                                                <  5:171041/149‑1 (MQ=255)        cAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCt                                                                                                                                <  3:175893/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                                                                                                              >  3:109818/1‑77 (MQ=255)                                                                                                  aaGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGc                                                                                                              <  4:109818/77‑1 (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)                                                                                                                                                     |                                                                                                                                       TGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGC  >  NC_002947/895063‑895344

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