Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,534,762	T→G	8.9%	intergenic (+49/‑106)	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,762	0	T	G	8.9%	170.9 / 6.8	56	intergenic (+49/‑106)	atoB/PP_3124	3‑oxoacid CoA‑transferase subunit B/short chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base T (29/22);  new base G (3/2);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

GCGCGAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCC  >  NC_002947/3534617‑3534909                                                                                                                                                  |                                                                                                                                                     gcgcGAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCgg                                                                                                                                                   >  7:306785/1‑149 (MQ=255)           gcgcCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAcc                                                                                                                                         >  6:174408/1‑149 (MQ=255)               cTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAccgccg                                                                                                                                     <  8:306785/149‑1 (MQ=255)                  gggTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCAt                                                                                                                                  >  5:289646/1‑149 (MQ=255)                       aGCGTCGAGGAAATCATCGAGAACACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTcgcg                                                                                                                             <  5:118788/149‑1 (MQ=255)                        gCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTcgc                                                                                                                              >  6:361564/1‑147 (MQ=255)                        gCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTcgc                                                                                                                              <  5:361564/147‑1 (MQ=255)                             gAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAg                                                                                                                      >  2:231946/1‑149 (MQ=255)                             gAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAg                                                                                                                      >  2:330395/1‑149 (MQ=255)                                    tcatcGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTc                                                                                                               <  2:148582/149‑1 (MQ=255)                                      atcGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTcc                                                                                                              >  6:111620/1‑148 (MQ=255)                                       tcGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTa                                                                                                            >  1:265005/1‑149 (MQ=255)                                       tcGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTa                                                                                                            >  1:363432/1‑149 (MQ=255)                                        cGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTAc                                                                                                           <  5:111620/149‑1 (MQ=255)                                         gagaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTAc                                                                                                           <  6:353981/148‑1 (MQ=255)                                         gagaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACa                                                                                                          <  2:265005/149‑1 (MQ=255)                                              aaCCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGttt                                                                                                     >  8:155101/1‑149 (MQ=255)                                                  gccgGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTg                                                                                                 <  1:152015/149‑1 (MQ=255)                                                     ggCAAACTGATCGTGCCGGATGATGTTAAGGAAATGACCTTCTAAGGCATGACGTGATCTAACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATAGAACCCCGCATTCGCGGGT‑AACCCTGTTCCTACAGGTTTTTTGCGt                                                                                              >  5:312631/1‑149 (MQ=38)                                                              aTCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTgg                                                                                       <  1:231946/148‑1 (MQ=255)                                                               tCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTACTACAGGTTTTTTGCGTTGCCTGGCaa                                                                                    >  6:94428/1‑149 (MQ=255)                                                                        atgatgTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTc                                                                           <  5:174408/149‑1 (MQ=255)                                                                                aaGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGAGTTTACCCCCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCAt                                                                    >  7:317574/1‑148 (MQ=255)                                                                                    aaaTGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAAc                                                               <  5:94428/149‑1 (MQ=255)                                                                                     aaTGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATTAGTCGGTAAAATTACCCCCGCAATAGCGGGT‑AAGCCTGTTAGAAAAAGTTTTTTTAGTTGACAGGCCATAAACCTTAATACCaaaaaaa                                                               >  2:349792/1‑147 (MQ=255)                                                                                      aTGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCa                                                             <  6:289646/149‑1 (MQ=255)                                                                                              cTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAACCCAATAAAAg                                                      >  7:273019/1‑148 (MQ=255)                                                                                              cTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAg                                                      >  2:301765/1‑148 (MQ=255)                                                                                              cTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAg                                                      >  1:232748/1‑148 (MQ=255)                                                                                                 aGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATc                                                                        >  1:259745/1‑127 (MQ=255)                                                                                                 aGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATc                                                                        <  2:259745/127‑1 (MQ=255)                                                                                                     aTGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCaa                                                            <  6:175643/135‑1 (MQ=255)                                                                                                     aTGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCaa                                                            >  5:175643/1‑135 (MQ=255)                                                                                                          gTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTg                                         >  4:256419/1‑149 (MQ=255)                                                                                                          gTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTg                                         <  1:255453/149‑1 (MQ=255)                                                                                                            gATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTgg                                        <  2:363432/148‑1 (MQ=255)                                                                                                             aTCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGc                                       <  3:256419/148‑1 (MQ=255)                                                                                                                tcACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCt                                   >  6:166472/1‑149 (MQ=255)                                                                                                                tcACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCt                                   >  3:164517/1‑149 (MQ=255)                                                                                                                  aCCCTGTAGGCGCGGGTTTACCCGCGAATGTGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTg                                  >  3:358649/1‑148 (MQ=255)                                                                                                                  aCCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTg                                  >  4:313838/1‑148 (MQ=255)                                                                                                                    ccTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAAc                                                               >  5:289596/1‑117 (MQ=255)                                                                                                                    ccTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAAc                                                               <  6:289596/117‑1 (MQ=255)                                                                                                                        tAGGCGCGGGTTTACCCGCGAATGTGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCa                           <  4:358649/149‑1 (MQ=255)                                                                                                                         attcgcgGGTAAACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCACATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCa                           >  8:155766/4‑148 (MQ=255)                                                                                                                                      cccGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcg             >  3:81175/1‑149 (MQ=255)                                                                                                                                       ccGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgc            >  8:249900/1‑149 (MQ=255)                                                                                                                                       ccGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCagc            >  2:44774/1‑149 (MQ=255)                                                                                                                                       ccGCGAATGCGGCGGTCAATGTACCGACGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATaaaa                                                       <  4:231867/106‑1 (MQ=25)                                                                                                                                       ccGCGAATGCGGCGGTCAATGTACCGACGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATaaaa                                                       >  3:231867/1‑106 (MQ=25)                                                                                                                                       ccGCGAATGCGGCGGTCAATGTACCGACGCATTCGCGGGT‑AAACCCGCGCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCaa                          <  2:291953/135‑1 (MQ=21)                                                                                                                                       ccGCGAATGCGGCGGTCAATGTACAGACGCATTCGCGGGT‑AAACCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTAAAACGTGGCAGGTGAAATTCaa                          >  7:54125/1‑135 (MQ=25)                                                                                                                                             aTGCGGCGGTCAATGTACCGACGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTGCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCtt      >  5:344348/1‑149 (MQ=38)                                                                                                                                              tGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCttt     <  3:313838/149‑1 (MQ=255)                                                                                                                                              tGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCttt     <  4:81175/149‑1 (MQ=255)                                                                                                                                                cGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGc   <  5:83283/149‑1 (MQ=255)                                                                                                                                                 ttCGGTACATTGACCTCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACTGCCGCTCCGCGCGCTTTGcc  <  6:312631/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GCGCGAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCC  >  NC_002947/3534617‑3534909

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