Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,534,793	G→A	7.4%	intergenic (+80/‑75)	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,793	0	G	A	7.4%	154.5 / 2.5	54	intergenic (+80/‑75)	atoB/PP_3124	3‑oxoacid CoA‑transferase subunit B/short chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base G (21/29);  new base A (2/2);  total (23/31)
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.54e-01

AATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCGCTGCTCCAACTGGG  >  NC_002947/3534650‑3534928                                                                                                                                                 |                                                                                                                                        actcatcGAGAAAACCGCCCGCCAACTGATCGGGCCGGATGATGGGAAGGAAAAGACCTGCTAAGGCACGAC‑GTGATCTCCCCCTGTAGGCGCGGGATTACCCGCGAAGGCGTCGGTACATAGACCGCCGCATTCGCGGGTAAGCCTGt                                                                                                                                    <  3:42860/147‑1 (MQ=255) aataatCGAGAAAACCGCCGGCAAACTGTTCGTGCCGGATTATTTGAAGGAAATTACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTATGCGGGGGTTTACCCGCGAATGCGGCGGTCAATGTACCGACGCATTCGCGGGTAAAcccgc                                                                                                                                    <  4:148500/149‑4 (MQ=21) aaTCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGt                                                                                                                                    <  4:69436/149‑1 (MQ=255)      tcGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTa                                                                                                                               <  5:62714/149‑1 (MQ=255)      tcGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTa                                                                                                                               >  8:33017/1‑149 (MQ=255)         agaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAg                                                                                                                            <  6:103074/149‑1 (MQ=255)         agaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAg                                                                                                                            >  4:269445/1‑149 (MQ=255)         agaAAACCGCCGGCAAACTGATCGGGCCGGATTATTTTAAGGAAATTACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAg                                                                                                                            <  4:190024/149‑1 (MQ=38)                 gccgGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGtttttgc                                                                                                                    <  1:47308/149‑3 (MQ=255)                 gccgGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTg                                                                                                                    <  1:185103/149‑1 (MQ=255)                        aaCTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGcc                                                                                                             <  1:100872/149‑1 (MQ=255)                             aTCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGc                                                                                                         >  4:247281/1‑148 (MQ=255)                                       atgatgTGAAGGAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCtt                                                                                               >  1:63947/1‑148 (MQ=255)                                       atgatgTGAAGGAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTc                                                                                              <  2:347808/149‑1 (MQ=255)                                             tgAAGGAAATGACCTTCTAAGGCATGACGGTGATATCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATcc                                                                                          >  1:28490/1‑148 (MQ=255)                                              gAAGGAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCa                                                                                        >  8:300686/1‑148 (MQ=255)                                                  gAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATaaaa                                                                                   >  2:59866/1‑149 (MQ=255)                                                  gAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATaaaa                                                                                   >  7:197166/1‑149 (MQ=255)                                                   aaaTGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAAc                                                                                  <  4:125439/149‑1 (MQ=255)                                                      tGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGc                                                                                                         >  3:229836/1‑123 (MQ=255)                                                      tGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGc                                                                                                         <  4:229836/123‑1 (MQ=255)                                                      tGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCaa                                                                               >  7:247393/1‑149 (MQ=255)                                                           ttCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAACGCCGCATTCGCCGGTAAGCCAGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCAAAAAACCAATaa                                                                            >  6:225057/1‑147 (MQ=255)                                                             cTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCt                                                                                                                                <  2:121570/93‑1 (MQ=255)                                                             cTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCt                                                                                                                                >  1:121570/1‑93 (MQ=255)                                                             cTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATaaaa                                                                          <  7:33017/147‑1 (MQ=255)                                                             cTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAg                                                                         >  5:206720/1‑148 (MQ=255)                                                                aGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAAc                                                                                  <  1:347407/136‑1 (MQ=255)                                                                aGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAAc                                                                                  >  2:347407/1‑136 (MQ=255)                                                                aGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAgg                                                                        <  3:247281/146‑1 (MQ=255)                                                                aGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGaa                                                                      <  3:269445/148‑1 (MQ=255)                                                                  gCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCt                                                                                                                                >  8:294287/1‑88 (MQ=255)                                                                  gCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCt                                                                                                                                <  7:294287/88‑1 (MQ=255)                                                                     tGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAAc                                                                <  7:327108/149‑1 (MQ=255)                                                                       aC‑GTGATCTCACCCTGTAGGCTCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAAcc                                                               <  4:153525/148‑1 (MQ=255)                                                                                             gcgGGTTTACCCGCGAATGCGTCGGTTAATTGAACGCCGCATTAGCGGGTAAGCCTGTTGCTAAAGGTTTTTTGCGTTGCCTGGCAATa                                                                                                     >  8:53103/1‑89 (MQ=25)                                                                                             gcgGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATa                                                                                                     <  7:53103/89‑1 (MQ=255)                                                                                               gggTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGc                                       <  7:300686/149‑1 (MQ=255)                                                                                               gggTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGc                                       <  8:247393/149‑1 (MQ=255)                                                                                               gggTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGc                                       <  6:206720/149‑1 (MQ=255)                                                                                                     aCCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAgg                                                                        >  7:174036/1‑110 (MQ=255)                                                                                                     aCCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAgg                                                                        <  8:174036/110‑1 (MQ=255)                                                                                                       ccGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGaa                                                                      <  6:252482/110‑1 (MQ=255)                                                                                                       ccGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGaa                                                                      >  5:252482/1‑110 (MQ=255)                                                                                                         gcgAATGCGGCGGTCAATGTACCGACGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgcgc                             >  5:264269/1‑149 (MQ=38)                                                                                                         gcgAATGCGGCGGTCAATGTACCGACGCATTCGCGGGTAAACCCGCGCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgcgc                             >  7:311188/1‑149 (MQ=25)                                                                                                         gcgAATGCGGCGGTCAATGTACCGACGCATTCGCGGGTAAACCCGCGCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAAAAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgcgc                             >  2:263812/1‑149 (MQ=25)                                                                                                           gAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCgcgcgc                           <  2:28490/149‑1 (MQ=255)                                                                                                            aaTGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCt                          <  8:197166/149‑1 (MQ=255)                                                                                                                         aTTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgc                               <  5:111148/131‑1 (MQ=255)                                                                                                                         aTTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgc                               >  6:111148/1‑131 (MQ=255)                                                                                                                         aTTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCgct              >  8:221052/1‑148 (MQ=255)                                                                                                                                   gcATTCGCGGGTAAACCCGCGCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCGCTGCTCCAACTgg   <  8:311188/149‑1 (MQ=255)                                                                                                                                    cATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCGCTGCTCCAACTggg  <  1:59866/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                        AATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGAC‑GTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCGCTGCTCCAACTGGG  >  NC_002947/3534650‑3534928

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