Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,671,122	T→A	8.3%	Q399L (CAG→CTG)	acnB ←	aconitate hydratase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,671,122	0	T	A	8.3%	128.1 / 3.6	48	Q399L (CAG→CTG)	acnB	aconitate hydratase B
Reads supporting (aligned to +/- strand):  ref base T (21/23);  new base A (2/2);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

TGTGGCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTCGGTCAGGCCGCGGCC  >  NC_002947/2670976‑2671270                                                                                                                                                    |                                                                                                                                                     tgtgGCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTgg                                                                                                                                                    <  4:307902/149‑1 (MQ=255)   tgGCAGAAGGACTGCATGACCAGGTCAGCAGAGAAGGCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGTGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGcc                                                                                                                                                  <  4:317083/149‑1 (MQ=255)     gCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCtt                                                                                                                                                <  6:39877/149‑1 (MQ=255)     gCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCtt                                                                                                                                                >  6:113632/1‑149 (MQ=255)     gCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCAGGGCGTGCGCCCAGGc                                                                                                                                                   <  8:43399/146‑1 (MQ=255)     gCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCAGGGCGTGCGCCCAGGc                                                                                                                                                   >  7:43399/1‑146 (MQ=255)                 cATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCca                                                                                                                                     >  7:64380/1‑148 (MQ=255)                 cATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCca                                                                                                                                     >  4:191309/1‑148 (MQ=255)                 cATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCACTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCca                                                                                                                                     >  1:236128/1‑148 (MQ=255)                  aTGACCAGGTCAGCAGAGAAGCCCAGGCATGCCAGGGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCca                                                                                                                                     >  1:231585/1‑148 (MQ=255)                  aTGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCca                                                                                                                                     >  6:195026/1‑147 (MQ=255)                      ccAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCa                                                                                                                               >  1:163854/1‑149 (MQ=255)                         ggTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACgg                                                                                                                            >  8:332725/1‑149 (MQ=255)                             agcagAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCt                                                                                                                                              <  8:96251/127‑1 (MQ=255)                             agcagAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCt                                                                                                                                              >  7:96251/1‑127 (MQ=255)                               cagAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAg                                                                                                                                         <  1:212100/130‑1 (MQ=255)                               cagAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAg                                                                                                                                         >  2:212100/1‑130 (MQ=255)                                  agaAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCAtt                                                                                                                   <  1:49463/149‑1 (MQ=255)                                         cAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTgcg                                                                                                            >  3:338465/1‑149 (MQ=255)                                                        cTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGccc                                                                                              >  4:335537/1‑149 (MQ=255)                                                         tttCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCt                                                                                             >  5:276460/1‑149 (MQ=255)                                                              gTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCg                                                                                        <  3:351912/149‑1 (MQ=255)                                                                        cgTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAg                                                                              <  1:30118/149‑1 (MQ=255)                                                                        cgTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAg                                                                              <  2:231585/149‑1 (MQ=255)                                                                         gTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGA‑CATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCGGgc                                                                            >  3:338187/1‑149 (MQ=255)                                                                                   cAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGt                                                                   <  4:338465/149‑1 (MQ=255)                                                                                   cAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGt                                                                   <  7:332725/149‑1 (MQ=255)                                                                                     gtggtgTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGttt                                                                 <  2:163854/149‑1 (MQ=255)                                                                                        gtgtCCTGGGAGCCGAGGGTGGTCATTTTTGGCTCGCTGTAGGCGGCTGGGCGGACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCtt                                                              <  2:236128/149‑1 (MQ=255)                                                                                         tgtCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTg                                                             <  6:276460/149‑1 (MQ=255)                                                                                          gtCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTg                                                             >  4:123218/1‑148 (MQ=255)                                                                                           tCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGaa                                                           <  8:64380/149‑1 (MQ=255)                                                                                              tGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAg                                                        >  1:108839/1‑149 (MQ=255)                                                                                                 gAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTc                                                     >  5:196499/1‑149 (MQ=255)                                                                                                 gAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTc                                                     <  5:113632/149‑1 (MQ=255)                                                                                                 gAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTc                                                     <  3:191309/149‑1 (MQ=255)                                                                                                                   tttGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAg                                   >  7:220685/1‑149 (MQ=255)                                                                                                                       gCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCg                               <  6:196499/149‑1 (MQ=255)                                                                                                                         tCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGGTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGc                             >  7:75316/1‑149 (MQ=255)                                                                                                                           gCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGc                             <  2:108839/147‑1 (MQ=255)                                                                                                                                    gcCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCttttt                  >  4:310745/1‑149 (MQ=255)                                                                                                                                    gcCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCttttt                  <  3:335537/149‑1 (MQ=255)                                                                                                                                    gcCAGGGCGTGCGCCCAGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGc                                                                                <  5:298977/87‑1 (MQ=255)                                                                                                                                    gcCAGGGCGTGCGCCCAGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGc                                                                                >  6:298977/1‑87 (MQ=255)                                                                                                                                     ggatctcgcACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTc                 <  3:343403/143‑1 (MQ=255)                                                                                                                                     ccTGGGCGCTCGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCAGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTc                 <  4:282294/149‑1 (MQ=255)                                                                                                                                         ggCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTCGGt              >  2:311820/1‑148 (MQ=255)                                                                                                                                                 cccTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTCGGTCAGGCcgcg     <  3:123218/149‑1 (MQ=255)                                                                                                                                                    tGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTCGGTCAGGCCGCGGcc  >  4:257545/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TGTGGCAGAACGACTGCATGACCAGGTCAGCAGAGAAGCCCAGGCATGCCA‑GGTCTTTCAGTTCGTCGCGCGTCATCGGGCCAGTGGTGTCCTGGGAGCCGACGGTGGTCATTTTTGGCTCGCAGTAGGCGCCTGGGCGCACGCCCTGGCCTTCTGGCAGGCCACAGGCACGGCCGACCATTTTCTGCGCCAGGGTGAAGCCCTTGCCGGAATCGGCAGGCTGCTCTGGTTTCTTGAACAGGTCGGAAGCACCCAGGCCCAGCTCGGCACGGGCTTTTTCGGTCAGGCCGCGGCC  >  NC_002947/2670976‑2671270

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