Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,329,767	T→C	8.9%	L40L (CTT→CTC)	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,767	0	T	C	8.9%	119.2 / 3.0	45	L40L (CTT→CTC)	PP_2048	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base C (2/2);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-01

CGTTTCGAACCGCGTGGGGAGGTTCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACGGCATGCGGGACGCCTATCGGCAATACGTCGAACATGGCTGGGC  >  NC_002947/2329619‑2329913                                                                                                                                                       |                                                                                                                                                   cGTTTCGAACCGCGTGGGGAGGTTCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCtt                                                                                                                                                    <  4:135162/149‑1 (MQ=255)   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)                                                                                                                                                       |                                                                                                                                                   CGTTTCGAACCGCGTGGGGAGGTTCAGGCATGAACTACCAGGCGCCCCTTCGCGACATGCGTTTCGTGCTGCATGAACTGTTCGACCTGAGCGGACATTGCCAGCAACTCGGTAT‑‑CGAACTGGACCGCGACACCCTCGACGGCATCCTTGAGGAGGCCGCACGTTTCACTGGCCAGATCATCGCACCGCTCAATCGCAACAGTGACGAACAGGGCTGTCAGTGGGTGGGCGGTGAGGTCACCACCCCGGACGGCATGCGGGACGCCTATCGGCAATACGTCGAACATGGCTGGGC  >  NC_002947/2329619‑2329913

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