Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,743,373	A→C	13.3%	intergenic (‑206/+41)	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,373	0	A	C	13.3%	133.8 / 14.5	53	intergenic (‑206/+41)	fbp/glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (30/16);  new base C (3/4);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

GAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCT  >  NC_002947/5743235‑5743520                                                                                                                                             |                                                                                                                                                    gAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgtg                                                                                                                                           <  8:166313/149‑1 (MQ=255)  aGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCGTTCGCGGGCACG‑‑CCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                                                          <  3:365565/149‑1 (MQ=18)  aGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                                                          <  3:292091/149‑1 (MQ=21)    tgatgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCACG‑‑CCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGa                                                                                                                                        <  3:57614/149‑1 (MQ=18)      atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAg                                                                                                                                       <  6:129047/148‑1 (MQ=255)      atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGc                                                                                                                                      >  2:342238/1‑149 (MQ=255)       tgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCGTGAAGCTTGCGAGACACCTGTGGGAGCg                                                                                                                                     >  4:252593/1‑149 (MQ=255)         agaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCgg                                                                                                                                    <  6:10619/148‑1 (MQ=255)              cGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCCCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGagcgagcgc                                                                                                                                <  1:105725/149‑2 (MQ=17)                ccACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGcac                                                                                                                            >  3:121427/1‑147 (MQ=255)                 cACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcg                                                                                                                           >  4:335220/1‑149 (MQ=255)                 cACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcg                                                                                                                           >  6:94651/1‑149 (MQ=255)                                         gTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCGTTGAAGCTGGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGCGCGGCCTCTCGTCAATCCAAt                                                                                                   >  2:158704/1‑149 (MQ=38)                                             gTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCt                                                                                               >  8:105267/1‑149 (MQ=255)                                                gTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTc                                                                                             <  7:105267/148‑1 (MQ=255)                                                          cGCCTGGCGGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑AGCGCTCCCCAAGGTACAGCGCAAGCCTTTAAGCTTGCGAGACACCTGTGGGAGCGGGGGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCgg                                                                                   <  1:342238/148‑1 (MQ=255)                                                           gCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGTAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTg                                                                                 <  3:335220/149‑1 (MQ=255)                                                           gCCAGGCCGGTGGGCTGGTTGGGGGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTg                                                                                 <  6:259818/149‑1 (MQ=255)                                                               ggCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATCCCTTCCAGATCTCGGTGGCat                                                                             >  5:62783/1‑149 (MQ=255)                                                                   ggtggtCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTc                                                                         <  7:157202/149‑1 (MQ=255)                                                                       gtCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGGGAGGCACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCTTATTcgcg                                                                      >  7:76998/1‑148 (MQ=255)                                                                            gTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGt                                                                >  3:133239/1‑149 (MQ=255)                                                                                 gtgtCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGCCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGt                                                           >  8:116087/1‑149 (MQ=255)                                                                                    tCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                               <  4:26942/142‑1 (MQ=255)                                                                                    tCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                               >  3:26942/1‑142 (MQ=255)                                                                                      cATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                                                      <  5:94651/149‑1 (MQ=255)                                                                                       aTTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCAGCTCAATCCAATGCCTTCCAGATCTCGGTGGGATATTCGCGAATGGTGCGGTCcgac                                                      >  5:123946/1‑148 (MQ=255)                                                                                           gcgGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGaa                                                 >  1:40275/1‑149 (MQ=255)                                                                                             gggCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATc                                                                                                       >  4:261930/1‑93 (MQ=255)                                                                                             gggCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATc                                                                                                       <  3:261930/93‑1 (MQ=255)                                                                                              gcggcg‑‑cTCCCTCCCACACGTACACCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCa                                              <  8:76998/146‑1 (MQ=255)                                                                                               gcacg‑‑cccgctCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGAGCGAGCGCGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAg                                             >  6:161439/4‑149 (MQ=17)                                                                                                 gcg‑‑cTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAAc                                                <  6:66774/144‑1 (MQ=255)                                                                                                 gcg‑‑cTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAAc                                                >  5:66774/1‑144 (MQ=255)                                                                                                        cgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTCCGAGACACATGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATcc                                      <  3:72106/149‑1 (MQ=255)                                                                                                         gctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccg                                     >  3:107289/1‑149 (MQ=255)                                                                                                           tcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcg                                                                        <  5:320522/112‑1 (MQ=255)                                                                                                           tcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcg                                                                        >  6:320522/1‑112 (MQ=255)                                                                                                           tcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtg                                   >  3:170462/1‑149 (MQ=255)                                                                                                             ccACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgcc                                 >  5:379820/1‑149 (MQ=255)                                                                                                             ccACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgcc                                 >  8:186032/1‑149 (MQ=255)                                                                                                               acaGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccgt                               >  5:379083/1‑149 (MQ=255)                                                                                                               acaGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCACATccgtgccgt                               >  4:83179/1‑149 (MQ=255)                                                                                                                   gTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGtt                            >  3:210430/1‑148 (MQ=255)                                                                                                                    tACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTc                           >  5:306867/1‑148 (MQ=255)                                                                                                                     gtctcgcAAGCTTCAAGGCTTGCGCGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCaa                         >  7:104062/5‑149 (MQ=38)                                                                                                                     gtctcgcAAGCTTCAAGGCTTGCGCGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCaa                         >  7:98887/5‑149 (MQ=38)                                                                                                                     aCAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGCGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCACTGCCGTGTTCaa                         >  5:47/1‑149 (MQ=255)                                                                                                                         cgcAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATAcc                     >  5:58828/1‑149 (MQ=255)                                                                                                                          gcAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCg                    >  4:313838/1‑149 (MQ=255)                                                                                                                                cTTGACGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTgc              <  3:252593/149‑1 (MQ=255)                                                                                                                                      gCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGccac         >  1:215058/1‑148 (MQ=255)                                                                                                                                        ttGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCcaccac      >  3:162676/1‑149 (MQ=255)                                                                                                                                            gagaCACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGAGCCGTGTTCAATACCGCCTTGCGCCACCACTCCt  >  4:125410/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                    GAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCG‑‑CTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCT  >  NC_002947/5743235‑5743520

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