Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,534,771	T→G	14.0%	intergenic (+58/‑97)	atoB → / → PP_3124	3‑oxoacid CoA‑transferase subunit B/short chain fatty acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,534,771	0	T	G	14.0%	109.1 / 11.3	43	intergenic (+58/‑97)	atoB/PP_3124	3‑oxoacid CoA‑transferase subunit B/short chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base T (18/19);  new base G (2/4);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

CGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCGCTGC  >  NC_002947/3534626‑3534918                                                                                                                                                  |                                                                                                                                                    cgcgCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAc                                                                                                                                                  >  5:236137/1‑149 (MQ=255)              aGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGGCGGTCAATGTACCGACGCATTcgcg                                                                                                                                     <  6:122394/149‑1 (MQ=38)              aGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGCTTACCCGCGAATGCGGCGGTCAATGTACCGACGCATTcgc                                                                                                                                      <  5:108171/148‑1 (MQ=38)                      ggAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGcc                                                                                                                             <  6:94160/149‑1 (MQ=255)                             atcGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACACTGTAGGCGCGGGTTTACACGCGAATGCGTCGGTACATTTACCGCCGCATTCGCGGGAAAGCCTGTTcc                                                                                                                       >  4:19154/1‑148 (MQ=255)                               cGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTa                                                                                                                     >  5:113026/1‑148 (MQ=255)                                gagaAAACCGCCGGCAAACTGATCGTGCCTGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTAc                                                                                                                    <  8:95735/148‑1 (MQ=255)                                gagaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACa                                                                                                                   <  7:175535/149‑1 (MQ=255)                                gagaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACa                                                                                                                   <  2:94033/149‑1 (MQ=255)                                 agaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAg                                                                                                                  <  6:236137/149‑1 (MQ=255)                                 agaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGAACAGGCTTACCCGCGAATGCGGCGGTCAATGTACCGACGCATTcgc                                                                                                                                      <  6:213656/129‑1 (MQ=21)                                 agaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGAACAGGCTTACCCGCGAATGCGGCGGTCAATGTACCGACGCATTcgc                                                                                                                                      >  5:213656/1‑129 (MQ=21)                                              cAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTg                                                                                                     >  3:137254/1‑149 (MQ=255)                                               aaaCTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGtttttt                                                                                                           >  3:63769/1‑142 (MQ=255)                                               aaaCTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGtttttt                                                                                                           <  4:63769/142‑1 (MQ=255)                                               aaaCTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGc                                                                                                    <  6:113026/149‑1 (MQ=255)                                                    gATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTgg                                                                                                >  5:158675/1‑148 (MQ=255)                                                     aTCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCa                                                                                              <  4:137254/149‑1 (MQ=255)                                                                   tgtgAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATcc                                                                                >  3:204261/1‑149 (MQ=255)                                                                       aaGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCAt                                                                             >  3:28979/1‑148 (MQ=255)                                                                        aGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATaa                                                                           <  6:158675/149‑1 (MQ=255)                                                                            aaTGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAAc                                                                        >  3:205153/1‑148 (MQ=255)                                                                                    tCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAg                                                               >  3:12441/1‑149 (MQ=255)                                                                                     cTAAGGCATGACGTGATCTCACCCTGTAGGCGCTGGTTTACCCGCGAATGCGTCGGTACATTGACCTCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATaaaa                                                                <  3:19154/147‑1 (MQ=255)                                                                                     cTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATaaaa                                                                <  4:28979/147‑1 (MQ=255)                                                                                                      ctcACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGaa                                                            <  8:246906/134‑1 (MQ=255)                                                                                                      ctcACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGaa                                                            >  7:246906/1‑134 (MQ=255)                                                                                                        cACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGGGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCt                                            >  3:57607/1‑148 (MQ=255)                                                                                                           ccTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGaaa                                        >  6:236471/1‑149 (MQ=255)                                                                                                             tGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTGCAATGACCGCCGCCGTCGCGGGTTAGCCTGGTCCTACAGGTTTTTTGCCTTGCCTGGCAATAACGCGTCATCCCATAAAACCAATaaaa                                                                >  5:176587/1‑123 (MQ=37)                                                                                                             tGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATaaaa                                                                <  6:176587/123‑1 (MQ=255)                                                                                                                 ggCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGTGTAAGCCTTTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCaa                                   <  1:51239/148‑1 (MQ=255)                                                                                                                 ggCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCaa                                   <  4:57607/148‑1 (MQ=255)                                                                                                                    gcgGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGa                                 >  1:72074/1‑147 (MQ=255)                                                                                                                              ccGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGt                                                                                                       >  3:135847/1‑67 (MQ=255)                                                                                                                              ccGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGt                                                                                                       <  4:135847/67‑1 (MQ=255)                                                                                                                              ccGCGAATGCGGCGGTCAATGTACCGACGCATTCGCGGGTAAACCCGCGCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCaa                                   <  5:269076/135‑1 (MQ=21)                                                                                                                              ccGCGAATGCGGCGGTCAATGTACCGACGCATTCGCGGGTAAACCCGCGCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCaa                                   >  6:269076/1‑135 (MQ=21)                                                                                                                               cgcgAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAg                                                               >  2:265997/1‑106 (MQ=255)                                                                                                                               cgcgAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAg                                                               <  1:265997/106‑1 (MQ=255)                                                                                                                                gcgAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgcgc                   >  7:63598/1‑149 (MQ=255)                                                                                                                                        gTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGcc           <  8:63598/149‑1 (MQ=255)                                                                                                                                            gTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTgc       <  2:72074/149‑1 (MQ=255)                                                                                                                                                 ttGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCgctgc  >  8:128961/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCGCTGC  >  NC_002947/3534626‑3534918

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