Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,845,721	C→T	8.8%	intergenic (‑88/+10)	PP_3394 ← / ← PP_3395	3‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,845,721	0	C	T	8.8%	126.6 / 2.5	46	intergenic (‑88/+10)	PP_3394/PP_3395	3‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (22/19);  new base T (2/2);  total (25/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 = 6.48e-01

ACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACC  >  NC_002947/3845575‑3845869                                                                                                                                                    |                                                                                                                                                     aCTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCTTCAGcc                                                                                                                                                     <  2:8922/149‑1 (MQ=255) aCTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGc                                                                                                                                                      >  5:367022/1‑148 (MQ=255)      cAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGa                                                                                                                                                 <  4:251206/148‑1 (MQ=255)         gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGa                                                                                                                                             <  3:322268/149‑1 (MQ=255)         gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGa                                                                                                                                             >  1:64960/1‑149 (MQ=255)             tGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAg                                                                                                                                         >  2:71281/1‑149 (MQ=255)                  gCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTGGCTGACGATGGATCa                                                                                                                                          >  8:212298/1‑143 (MQ=39)                  gCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTGGCTGACGATGGATCa                                                                                                                                          <  7:212298/143‑1 (MQ=39)                                       attgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACg                                                                                                               >  5:298152/1‑149 (MQ=255)                                         tgattgTTATTGTTGTCATGGCGGGATCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCGGTATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCAATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACgc                                                                                                              <  6:367022/148‑1 (MQ=255)                                           attgttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTGGCTGACGATGGATCAGccc                                                                                                                                      <  2:15512/122‑1 (MQ=37)                                               ttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCaa                                                                                                         >  8:216247/1‑147 (MQ=255)                                                 attgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAgg                                                                                                     >  8:103976/1‑149 (MQ=255)                                                    gttgtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATACATAGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtgg                                                                                                  >  2:328634/1‑149 (MQ=255)                                                             gcgGTGTCTCTGTCTGAAGCGTTGCGCCTGACTAATAGCTGGATGAGGGGCAGAAAAGACGTGTTTGGCGGGCTTATCCATGGTCAGCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAgg                                                                                         <  6:298152/149‑1 (MQ=255)                                                                cgctctctGTATGACGCGTTGCGCCAGAGTAAGAGCTGGATGACGGGCAGACAAGACGTGTTTCGCGGGCTGATCCAGCGTCAGCCACGACGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGgctc                                                                                      <  8:154361/147‑3 (MQ=255)                                                                     tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGAAAAAGGTGGTCAAGCAGGCgccgagc                                                                                 >  7:231417/1‑149 (MQ=255)                                                                      ctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCa                                                                                                                           >  2:53651/1‑106 (MQ=255)                                                                      ctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCa                                                                                                                           >  2:59107/1‑106 (MQ=255)                                                                      ctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCa                                                                                                                           <  1:59107/106‑1 (MQ=255)                                                                      ctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCa                                                                                                                           <  1:53651/106‑1 (MQ=255)                                                                           tgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGt                                                                                                                    >  3:150352/1‑108 (MQ=255)                                                                           tgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCAGATCCATCGTCAGCCAAGATGGATCAGCCCGCGCCCAGCCAGGCCGGTGCACCCAGAAACATGTCGTAAAGCAGGCCCACAGCCGACAg                                                                           >  1:322251/1‑149 (MQ=255)                                                                             aaGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtgg                                                                                                  >  4:33147/1‑124 (MQ=255)                                                                             aaGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtgg                                                                                                  <  3:33147/124‑1 (MQ=255)                                                                             aaGCGTTCCGCCAGACTAAGAGCTGTATGACGGGCATAAACGACGTGTTTCGCGGGCTTACCCATCGCCCCCCACTTTGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                                                                          <  6:1093/148‑1 (MQ=255)                                                                                cGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCgaga                                                                      >  6:374251/1‑149 (MQ=255)                                                                                cGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCgaga                                                                      <  7:103976/149‑1 (MQ=255)                                                                                          gACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACTATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCg                                                            <  2:64960/149‑1 (MQ=255)                                                                                                        tGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTc                                                         >  4:103127/1‑138 (MQ=255)                                                                                                        tGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTc                                                         <  3:103127/138‑1 (MQ=255)                                                                                                            gggCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgt                                                                  >  4:284229/1‑125 (MQ=255)                                                                                                            gggCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgt                                                                  <  3:284229/125‑1 (MQ=255)                                                                                                                aGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCAAGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGGCGAAt                                      >  1:105638/1‑149 (MQ=255)                                                                                                                   aaaGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTc                                    <  5:374251/148‑1 (MQ=255)                                                                                                                     aGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAc                                  >  5:343007/1‑148 (MQ=255)                                                                                                                      gACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCa                                <  7:216247/149‑1 (MQ=255)                                                                                                                         gggtTTCGCGGGCTGATCCATGGTCGGCCACTATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgc                             <  5:348496/147‑1 (MQ=255)                                                                                                                                      tGATCCATCGTGGCTGACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtcag                >  5:263411/1‑149 (MQ=255)                                                                                                                                      tGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCgag                                                                       >  1:199264/1‑94 (MQ=255)                                                                                                                                      tGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCgag                                                                       <  2:199264/94‑1 (MQ=255)                                                                                                                                       gATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtcagt               >  3:174994/1‑149 (MQ=255)                                                                                                                                        aTCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGCGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtcagt               >  6:129310/1‑148 (MQ=255)                                                                                                                                         tCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGCGCGTGAGCCCGTCGCGCATGCAGATGACGAATTCCCCAAGCGCCCAGACggcagt               >  4:11502/1‑147 (MQ=255)                                                                                                                                           cATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAgc           >  8:73759/1‑149 (MQ=255)                                                                                                                                            aTCGTCAGACACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAgcg          >  1:218161/1‑149 (MQ=255)                                                                                                                                                tCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAAt      <  7:73759/149‑1 (MQ=255)                                                                                                                                                   gCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGAc   <  1:328634/149‑1 (MQ=255)                                                                                                                                                    ccACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtcag                >  4:152969/1‑135 (MQ=255)                                                                                                                                                    ccACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGAcc  >  5:367594/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACC  >  NC_002947/3845575‑3845869

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