Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,534,767	A→C	13.9%	intergenic (+54/‑101)	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,767	0	A	C	13.9%	108.0 / 9.9	45	intergenic (+54/‑101)	atoB/PP_3124	3‑oxoacid CoA‑transferase subunit B/short chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base A (17/21);  new base C (2/4);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

CGAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGC  >  NC_002947/3534620‑3534913                                                                                                                                                    |                                                                                                                                                   cgAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTAc                                                                                                                                                   <  6:118727/149‑1 (MQ=255)   agaCCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACAt                                                                                                                                                 <  6:192204/149‑1 (MQ=255)       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                                                                                                              <  2:94033/149‑1 (MQ=255)                                      gagaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACa                                                                                                              <  7:175535/149‑1 (MQ=255)                                       agaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAg                                                                                                             <  6:236137/149‑1 (MQ=255)                                       agaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGAACAGGCTTACCCGCGAATGCGGCGGTCAATGTACCGACGCATTcgc                                                                                                                                 >  5:213656/1‑129 (MQ=21)                                       agaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGAACAGGCTTACCCGCGAATGCGGCGGTCAATGTACCGACGCATTcgc                                                                                                                                 <  6:213656/129‑1 (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)                                                                                                                                                    |                                                                                                                                                   CGAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGC  >  NC_002947/3534620‑3534913

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