Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	994,960	G→A	9.1%	intergenic (+39/+86)	PP_0859 → / ← PP_0860	ketoglutaramate omega‑amidase/sulfite reductase flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	994,960	0	G	A	9.1%	155.6 / 4.4	55	intergenic (+39/+86)	PP_0859/PP_0860	ketoglutaramate omega‑amidase/sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (22/28);  new base A (2/3);  total (24/31)
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.37e-01

GGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTG  >  NC_002947/994814‑995090                                                                                                                                                    |                                                                                                                                   ggggAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggc                                                                                                                                  >  6:215955/1‑149 (MQ=255)    gAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccgg                                                                                                                               <  3:290039/149‑1 (MQ=255)      ggCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggccggcc                                                                                                                             <  8:63063/149‑1 (MQ=255)            cGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgcca                                                                                                                       <  7:57820/149‑1 (MQ=255)                 tgtTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCCCTTTCGCGGGCAAGCCCGCGAAGAAGCCGGGCCTgctgtcagt                                                                                                                  <  4:140756/149‑8 (MQ=25)                   ttCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAAtt                                                                                                                >  3:81427/1‑149 (MQ=255)                      aCTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTc                                                                                                              >  1:114678/1‑148 (MQ=255)                             tACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGTCAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGc                                                                                                       <  2:292761/148‑1 (MQ=255)                               cTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCg                                                                                                     <  5:215955/148‑1 (MQ=255)                                    cgcgGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCAcc                                                                                               <  1:45716/149‑1 (MQ=255)                                           gAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTgccg                                                                                        >  3:82436/1‑149 (MQ=255)                                            aaCTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCAGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGtt                                                                                            >  1:66836/1‑145 (MQ=255)                                            aaCTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGtt                                                                                            <  2:66836/144‑1 (MQ=255)                                               tGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTgccgccgc                                                                                    >  5:133980/1‑149 (MQ=255)                                                 gcggcgCACCGGGCCAAGTTCCCGCCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACACCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTgccgccgc                                                                                    <  4:232153/147‑1 (MQ=255)                                                         ccGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgccagc                                                                                                                     >  1:32386/1‑106 (MQ=255)                                                         ccGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTgccagc                                                                                                                     <  2:32386/106‑1 (MQ=255)                                                            ggCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGcc                                                                                                      >  6:176187/1‑118 (MQ=255)                                                            ggCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGcc                                                                                                      <  5:176187/118‑1 (MQ=255)                                                             gCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                      <  7:80260/149‑1 (MQ=255)                                                             gCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                      <  4:225960/149‑1 (MQ=255)                                                             gCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                      <  2:114678/149‑1 (MQ=255)                                                               cAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGGCCGGCCCTTTCGCGGGCAAGCCCGCGAAGAAGCCGGGCCTgct                                                                                                                        >  6:97398/1‑96 (MQ=17)                                                               cAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGGCCGGCCCTTTCGCGGGCAAGCCCGCGAAGAAGCCGGGCCTgct                                                                                                                        <  5:97398/97‑2 (MQ=17)                                                                aaGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTa                                                                   <  1:71285/149‑1 (MQ=255)                                                                      ccGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCAACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCg                                                             >  5:98096/1‑149 (MQ=255)                                                                           cAATCTGGATGCC‑GATACCTTCGAGATGCACTGACAGCAGGCACGGCTTCTTCGCGGGCTTGCCCTCGAAAGGGCCGGCCCTGACAGCCAATTTCAACCTGCCGTTCACCCTTGCCGCCGCCACCGTATTCGAGGGTAAACCCGcaca                                                         >  4:43116/1‑145 (MQ=25)                                                                            aaTCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTgccg                                                                                        <  5:184733/116‑1 (MQ=255)                                                                            aaTCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTgccg                                                                                        >  6:184733/1‑116 (MQ=255)                                                                             aTCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTac                                                      <  2:144338/149‑1 (MQ=255)                                                                             aTCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTa                                                       >  6:9998/1‑148 (MQ=255)                                                                                  gATGGC‑GATACCTTCGAGGTGCACTGACAGCAGGCCTGGCTTCTTCGCTGGGTTGCCCGAGAAAGGGCAGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGTTTAAACCCGCTCCTacacac                                                  <  6:98096/148‑1 (MQ=21)                                                                                       c‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATc                                             >  7:137146/1‑148 (MQ=255)                                                                                             ccTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATaca                                       <  5:9998/149‑1 (MQ=255)                                                                                             ccTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATac                                        <  4:82436/148‑1 (MQ=255)                                                                                                       gCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGcc                              >  6:134008/1‑148 (MQ=255)                                                                                                        cacaGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGcc                              <  6:250465/147‑1 (MQ=255)                                                                                                        cACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCaa                                           >  7:137323/1‑134 (MQ=255)                                                                                                        cACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCaa                                           <  8:137323/134‑1 (MQ=255)                                                                                                          cTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGata                          <  8:137146/149‑1 (MQ=255)                                                                                                              cagcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGc                      <  5:134008/149‑1 (MQ=255)                                                                                                              cagcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGc                      <  4:81427/149‑1 (MQ=255)                                                                                                               agcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCAACCGGCCGTTCACCGTTGCGGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGcc                     >  7:86289/1‑149 (MQ=255)                                                                                                                gcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGAGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCggg                                                                     <  8:28349/100‑1 (MQ=255)                                                                                                                gcagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGAGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCggg                                                                     >  7:28349/1‑100 (MQ=255)                                                                                                                 cagGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGccc                    >  4:15800/1‑148 (MQ=255)                                                                                                                    gCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATcc                >  5:133071/1‑149 (MQ=255)                                                                                                                      ccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTc              >  7:20884/1‑149 (MQ=255)                                                                                                                       cGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATaca                                       <  6:190836/123‑1 (MQ=255)                                                                                                                       cGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATaca                                       >  5:190836/1‑123 (MQ=255)                                                                                                                         gCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGAt           <  1:3652/149‑1 (MQ=255)                                                                                                                           ttcttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGAt           <  2:54631/147‑1 (MQ=255)                                                                                                                           ttcttcGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGAt           >  1:54631/1‑147 (MQ=255)                                                                                                                                cgcgGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActgc    >  7:263214/1‑149 (MQ=255)                                                                                                                                cgcgGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActg     <  7:55653/148‑1 (MQ=255)                                                                                                                                  cgGGCTTGCCCGAGAAAGGGCCGTCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActgctg  <  3:43116/149‑1 (MQ=255)                                                                                                                                               gAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCggg                                                                     >  2:33318/1‑69 (MQ=255)                                                                                                                                               gAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCggg                                                                     <  1:33318/69‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                   GGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCC‑GATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCTG  >  NC_002947/994814‑995090

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