Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,743,335	T→C	8.3%	intergenic (‑168/+79)	fbp ← / ← glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,743,335	0	T	C	8.3%	131.5 / 2.3	48	intergenic (‑168/+79)	fbp/glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (19/25);  new base C (2/2);  total (21/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 = 8.31e-01

GGGGGATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATC  >  NC_002947/5743188‑5743483                                                                                                                                                     |                                                                                                                                                      gggggATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGcgctc                                                                                                                                                      >  8:250648/1‑149 (MQ=255)     gATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGTCTAGGGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCgctc‑gctc                                                                                                                                                 <  2:244665/149‑1 (MQ=255)        gAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCCC‑GCTCCca                                                                                                                                              >  7:68509/1‑149 (MQ=255)           aaaCCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAg                                                                                                                                           >  7:222923/1‑149 (MQ=255)           aaaCCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAg                                                                                                                                           >  7:280071/1‑149 (MQ=255)                         tACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGcc                                                                                                                             >  6:255937/1‑149 (MQ=255)                         tACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGc                                                                                                                              <  8:142976/148‑1 (MQ=255)                         tACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGc                                                                                                                              <  7:153476/148‑1 (MQ=255)                           cAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCtt                                                                                                                           <  8:280071/149‑1 (MQ=255)                           cAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCtt                                                                                                                           <  5:202621/149‑1 (MQ=255)                                 gagTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCt                                                                                                                     >  2:84479/1‑149 (MQ=255)                                    tGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGc                                                                                                                  <  8:176746/149‑1 (MQ=255)                                      aCGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                <  8:222923/149‑1 (MQ=255)                                      aCGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                <  3:192747/149‑1 (MQ=255)                                      aCGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                <  4:56372/149‑1 (MQ=255)                                      aCGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCg                                                                                                                 <  7:250648/148‑1 (MQ=255)                                        gggCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCgaga                                                                                                              <  1:294796/149‑1 (MQ=255)                                            tAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCgctgtacc                                                                                                          >  7:112013/1‑142 (MQ=25)                                            tAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTg                                                                                                                          >  4:247290/1‑133 (MQ=255)                                            tAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTg                                                                                                                          <  3:247290/133‑1 (MQ=255)                                            tAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTGCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGCCCATTCGCTCGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGacac                                                                                                           <  1:84479/149‑1 (MQ=255)                                                gAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgtg                                                                                                      <  5:255937/149‑1 (MQ=255)                                                gAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGT‑GTCTGGTTGTGTGTCCATTCGC‑GGGCACGCCC‑GCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                     <  2:127323/149‑1 (MQ=18)                                                     atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGc                                                                                                                              <  5:201584/120‑1 (MQ=255)                                                     atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGc                                                                                                                              >  6:201584/1‑120 (MQ=255)                                                     atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGc                                                                                                 >  3:236987/1‑149 (MQ=255)                                                      tgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGATCTGGTTGTGTGTCCATTCGC‑GGGCACGCCC‑GCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                     <  1:189265/144‑1 (MQ=18)                                                      tgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGATCTGGTTGTGTGTCCATTCGC‑GGGCACGCCC‑GCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                     >  2:189265/1‑144 (MQ=18)                                                        agaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCCGTGGGAGCga                                                                                               <  8:112013/148‑2 (MQ=21)                                                          aCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCg                                                                                            >  7:203493/1‑149 (MQ=255)                                                               ccACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGccc                                                                                       >  3:257986/1‑149 (MQ=255)                                                                 aCGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGGCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcgc                                                                                     <  7:301037/149‑1 (MQ=255)                                                                    gCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGAGAGACACCTGTGGGAGCGGGCGTACCcgaac                                                                                   >  6:92021/1‑145 (MQ=255)                                                                    gCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCAATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTACCcgaac                                                                                   >  6:96796/1‑145 (MQ=255)                                                                      tctcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTAGAACCTTGCGAGACACCTGTGGGAGCGCGCGTGCCCGCGAACg                                                                                >  6:264262/1‑149 (MQ=255)                                                                          aGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcct                                                                            >  7:299883/1‑149 (MQ=255)                                                                            aCGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctca                                                                          >  3:147991/1‑149 (MQ=255)                                                                                             tAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCtt                                                         <  4:129546/149‑1 (MQ=255)                                                                                                ttCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTCGGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTc                                                        >  1:93614/1‑148 (MQ=255)                                                                                                 tCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCa                                                      <  8:299883/148‑1 (MQ=255)                                                                                                       ttCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCg                                               >  7:161828/1‑149 (MQ=255)                                                                                                          gCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTg                                            <  2:93614/149‑1 (MQ=255)                                                                                                               gCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTAGTTGGCata                                       >  6:314382/1‑149 (MQ=255)                                                                                                                    tggtCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCAGAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGACATATTcgc                                  <  4:236987/149‑1 (MQ=255)                                                                                                                         tGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCACCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATg                             >  1:35484/1‑149 (MQ=255)                                                                                                                           gttgttTGTCCTTTGGG‑TGGCGCGGTC‑GCTCCCACAGGTACAGCTCTAGCATTGAAGTTTGTGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATgg                            <  1:134678/148‑1 (MQ=255)                                                                                                                            ttGTGTCGCCCTTCGC‑GGGCGCGCTC‑GCTCCCACCGGTACAGCCCAAGCCTTGAAGCTTGCGAGACACCTGTTGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGt                           <  4:257986/148‑1 (MQ=38)                                                                                                                                         cgc‑gGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACgag              <  2:35484/148‑1 (MQ=255)                                                                                                                                          gc‑gGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACgaga             <  8:203493/148‑1 (MQ=255)                                                                                                                                                cgcggTC‑GCTCCCAGAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCGGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGc       <  5:264262/149‑1 (MQ=255)                                                                                                                                                   gctc‑gctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCa    >  3:27610/1‑149 (MQ=255)                                                                                                                                                     tc‑gctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATc  >  6:52934/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GGGGGATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGC‑GGGCGCGCTC‑GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATC  >  NC_002947/5743188‑5743483

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑