Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,329,768	G→A	8.9%	E41K (GAG→AAG)	PP_2048 →	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,329,768	0	G	A	8.9%	122.1 / 2.5	45	E41K (GAG→AAG)	PP_2048	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (23/18);  new base A (2/2);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

TTTCGAACCGCGTGGGGAGGTTCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACGGCATGCGGGACGCCTATCGGCAATACGTCGAACATGGCTGGGCC  >  NC_002947/2329621‑2329914                                                                                                                                                      |                                                                                                                                                   tttCGAACCGCGTGGGGAGGTTCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTga                                                                                                                                                   <  2:113025/149‑1 (MQ=255)      aaCCGCGTGGGGAGGTTCAGGCATGAACTACCAGGCGCACCTTCGCGAAATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGACAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCAGCCTTGaggagg                                                                                                                                              >  4:316082/1‑149 (MQ=255)         cgcgTGGGGAGGTTCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCCTCCTCAAGGATGCCg                                                                                                                                           <  4:57800/149‑1 (MQ=255)            gTGGGGAGGTTCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCAc                                                                                                                                        <  1:71107/149‑1 (MQ=255)             tGGGGAGGTTCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCCTCCTCAAGGATGCCgtcga                                                                                                                                       <  2:205901/149‑5 (MQ=255)               gggAGGTTCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGtt                                                                                                                                     <  4:265126/149‑1 (MQ=255)                ggAGGTTCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGttt                                                                                                                                    <  3:178707/149‑1 (MQ=255)                    gtacAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGAAACTCGCTAT‑‑CGAACTGGACAGCGACACACTCGACGGAATCCTTGAGGAGGCCGAACGATTCAc                                                                                                                                 >  8:39876/4‑148 (MQ=255)                    gTTCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCTTTTCTTGCTGCATGAACTGTTCGACCTGATCGTACATTGCCATCAACTCTTTAT‑‑CTAACTGGACCGCTACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCa                                                                                                                                  <  7:39876/147‑1 (MQ=255)                    gTTCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACt                                                                                                                                >  3:34343/1‑149 (MQ=255)                     ttCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTg                                                                                                                               <  6:34321/149‑1 (MQ=255)                     ttCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTg                                                                                                                               <  4:306092/149‑1 (MQ=255)                            aTGAACTACCAGGCGCCCCTTCACGACATGCGTTTCGAGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTAGGTAT‑‑CGAACTGGAACGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACCTTTCACTGGc                                                                                                                             >  7:234920/1‑144 (MQ=255)                               aaCTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGatcat                                                                                                                     >  3:321/1‑149 (MQ=255)                                 cTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCg                                                                                                                   >  4:108654/1‑149 (MQ=255)                                    ccAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAAGTCGGTAT‑‑CGAACTGGACCGCGACTCCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCAc                                                                                                                <  3:316082/149‑1 (MQ=255)                                     cAGGCGCCCCTTCGCGTCATGCGTTTCGTGCTGCTTTAACTGTTCGACCTTTTCGGACATTGCCAGCATCTCGGTTT‑‑CGTTCTGGTCCGCGACATCTTCGACGGCATCCTTGAGGAGGTCGGATGTTTCTCTGGCCAGATCATCGCAc                                                                                                                <  8:234920/148‑1 (MQ=255)                                                  gcgACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCaa                                                                                                   <  4:331275/148‑1 (MQ=255)                                                     aCATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTATAGAGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACa                                                                                                 >  1:302555/1‑149 (MQ=255)                                                                gTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAggg                                                                                    <  3:108654/149‑1 (MQ=255)                                                                 tgctgcATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGc                                                                                   >  5:169224/1‑149 (MQ=255)                                                                       aTGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAg                                                                             <  5:151264/149‑1 (MQ=255)                                                                          atcTGTTCGACCTGAGCGGACATTGCCAGCAACTCTGTAT‑‑CGAACTTGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCa                                                                                                                                  <  4:64683/91‑1 (MQ=255)                                                                          aaCTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCa                                                                                                                                  >  3:64683/1‑93 (MQ=255)                                                                          aaCTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAgtgg                                                                          <  6:82814/149‑1 (MQ=255)                                                                               ttCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGAc                                                                                            >  8:60428/1‑126 (MQ=255)                                                                               ttCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGAc                                                                                            <  7:60428/126‑1 (MQ=255)                                                                                tCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGc                                                                    >  6:238528/1‑149 (MQ=255)                                                                                              cATTGCAAGCAACTCGTTAT‑‑CGACCTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCACCCTTGACGAACAGGGCTGTCTGTGGGTGGGCGGTGAGGTcaccac                                                      <  4:34343/149‑1 (MQ=255)                                                                                               aTTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCaccacc                                                     >  7:350096/1‑149 (MQ=255)                                                                                               aTTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCaccacc                                                     >  7:345313/1‑149 (MQ=255)                                                                                                    cAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGTGCGGTGAGGTCACCACCCCGGa                                                >  5:325802/1‑149 (MQ=255)                                                                                                     aGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGAc                                               <  2:302555/149‑1 (MQ=255)                                                                                                       cAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACgg                                             >  5:340169/1‑149 (MQ=255)                                                                                                        aaCTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACGGc                                            <  3:82189/149‑1 (MQ=255)                                                                                                               tAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACGGCATGCgg                                      >  3:257554/1‑148 (MQ=255)                                                                                                                      aaCTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTACCCACCCCGGACGGCAGGCTGGACGcc                                >  8:265641/1‑149 (MQ=255)                                                                                                                           gACCGCGACACCCTCGTCGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCTCATCATCGCTCCGCTCAATCGCAACTGTGTCGATCAGGGCTGTCAGTGGGTGGGCGGTGAGGTCAGCACCCGGGACGGCATACGGGACGCCTATc                            <  8:182780/148‑1 (MQ=255)                                                                                                                           gACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACGGCATGCGGGACGCCTATCg                           >  7:32750/1‑149 (MQ=255)                                                                                                                             ccGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACGGCATGCGGGACGCCTATCgg                          >  1:214484/1‑148 (MQ=255)                                                                                                                                      ccTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACGGCATGCGGGACGCCTATCGGCAATACGTc                 >  2:137800/1‑148 (MQ=255)                                                                                                                                      ccTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACGGCATGCGGGACGCCTATCGGAAATActac                 >  8:313437/1‑145 (MQ=255)                                                                                                                                         gtgcGGCCTCCTCAAGGATGCCGCACGTTTCACTGGGCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACGGCATGCGGGACGCCTATCGGCAATACGTCGaa              >  6:224369/4‑148 (MQ=255)                                                                                                                                         gtgcGGCCTCCTCAAGGATGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACGGCATGCGGGACGCCTATCGGCAATACGTCGaa              >  6:229505/4‑148 (MQ=255)                                                                                                                                             ggCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGATGGCATGCGGGACGCCTATCGGCAATACGTCGAACATgg         >  7:102363/1‑149 (MQ=255)                                                                                                                                                    ttGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACGGCATGCGGGACGCCTATCGGCAATACGTCGAACATGGCTGGGc   >  3:153828/1‑148 (MQ=255)                                                                                                                                                      gaggagGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACGGCATGCGGGACGCCTATCGGCAATACGTCGAACATGGCTGGGcc  <  4:257554/147‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TTTCGAACCGCGTGGGGAGGTTCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACGGCATGCGGGACGCCTATCGGCAATACGTCGAACATGGCTGGGCC  >  NC_002947/2329621‑2329914

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