Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,747,964	C→T	7.7%	intergenic (‑203/‑333)	etfB ← / → PP_4203	electron transfer flavoprotein subunit beta/ubiquinone oxidoreductase electron transfer flavoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,747,964	0	C	T	7.7%	150.9 / 2.0	52	intergenic (‑203/‑333)	etfB/PP_4203	electron transfer flavoprotein subunit beta/ubiquinone oxidoreductase electron transfer flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (25/23);  new base T (2/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

TGCGCGTGGGTACAAGGGCACCTCTAAAAACGTATCGGGAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTAACACCCATTCCGGTCAATACGGCAAATCGGTCACCCTCCAGCCGCGTT  >  NC_002947/4747817‑4748107                                                                                                                                                    |                                                                                                                                                 tGCGCGTGGGTACAAGGGCACCTCTAAAAACGTATCGGGAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGcc                                                                                                                                                 >  6:318264/1‑149 (MQ=255)    gcgTGGGTACAAGGGCACCTCTAAAAACGTATCGGGAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGc                                                                                                                                              <  4:196985/149‑1 (MQ=255)           tACAAGGGCACCTCTAAAAACGTATCGGGAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAAccc                                                                                                                                       <  3:207967/149‑1 (MQ=255)                gggCACCTCTAAAAACGTATCGGGAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCctcctc                                                                                                                                  >  2:110534/1‑149 (MQ=255)                       tctAAAAACGTATCGGGAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCt                                                                                                                          >  5:339830/1‑149 (MQ=255)                         tAAAAACGTATCGGGAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTgg                                                                                                                        >  1:132342/1‑149 (MQ=255)                          aaaaaCGTATCGGGAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGt                                                                                                                       >  3:108096/1‑149 (MQ=255)                           aaaaCGTATCGGGAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTc                                                                                                                      >  4:164509/1‑149 (MQ=255)                             aaCGTATCGGGAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCtt                                                                                                                    <  4:108096/149‑1 (MQ=255)                              aCGTATCGGGAGGCAGCAACTTCACAGTGACCAAACAGTCTAGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCttt                                                                                                                   <  5:318264/149‑1 (MQ=255)                               cGTATCGGGAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTa                                                                                                                  <  6:339830/149‑1 (MQ=255)                                   tCGGGAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTagag                                                                                                               <  3:346387/148‑1 (MQ=255)                                     gggAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGgtg                                                                                                            >  1:168372/1‑149 (MQ=255)                                       gAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTa                                                                                                          <  5:32331/149‑1 (MQ=255)                                                  ttCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCgg                                                                                               >  3:134664/1‑149 (MQ=255)                                                   tCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGAAAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAAACGACCTAAGCACGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGAGTACCTGCGCCCGGc                                                                                              >  5:277954/1‑149 (MQ=255)                                                     acaGTGACCAAACAGTCTTGTTCGGACTCTTGTGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCt                                                                                             <  3:184904/148‑1 (MQ=255)                                                        gTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCcac                                                                                                                               <  7:13256/112‑1 (MQ=255)                                                        gTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCcac                                                                                                                               >  8:13256/1‑112 (MQ=255)                                                                        gTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAAc                                                                            <  4:134664/146‑1 (MQ=255)                                                                         ttCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGTCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTa                                                                         >  6:36017/1‑148 (MQ=255)                                                                         ttCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGTCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTa                                                                         >  2:362937/1‑148 (MQ=255)                                                                         ttCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTa                                                                         >  6:91642/1‑148 (MQ=255)                                                                           cGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTaaaa                                                                      >  4:359393/1‑149 (MQ=255)                                                                              aCTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGaa                                                                             <  7:201884/139‑1 (MQ=255)                                                                              aCTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGaa                                                                             >  8:201884/1‑139 (MQ=255)                                                                              aCTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCg                                                                   >  2:197391/1‑149 (MQ=255)                                                                                  ttGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCg                                                               >  3:269301/1‑149 (MQ=255)                                                                                  ttGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCg                                                               >  8:253501/1‑149 (MQ=255)                                                                                  ttGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCg                                                               >  7:87157/1‑149 (MQ=255)                                                                                     cGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTAt                                                            <  7:253501/149‑1 (MQ=255)                                                                                            tCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGAcc                                                     <  2:132342/149‑1 (MQ=255)                                                                                            tCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGAcc                                                     >  8:148338/1‑149 (MQ=255)                                                                                              cTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCg                                                    >  3:334341/1‑148 (MQ=255)                                                                                                 aCGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTaa                                                 <  7:148338/148‑1 (MQ=255)                                                                                                    ggTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGa                                                       <  8:214652/139‑1 (MQ=255)                                                                                                    ggTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGa                                                       >  7:214652/1‑139 (MQ=255)                                                                                                             cTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTAACACCCATTCCGGt                                    <  3:359393/149‑1 (MQ=255)                                                                                                                tCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTAACACCCATTCCGGTc                                   <  1:197391/147‑1 (MQ=255)                                                                                                                tCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGATCGTATATTGACCGTAACACCCATTCCGGTCaa                                 >  3:148464/1‑149 (MQ=255)                                                                                                                          tAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACAGCCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCATACGGCGAACGTAAAACCGAACGTATCTTGACCGTAACACCCATTCCGGACAATACGGCaaa                        >  1:81677/1‑149 (MQ=255)                                                                                                                            aaCTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTAACACCCATTCCGGTCAATACGGCAAATCg                     <  2:81677/149‑1 (MQ=255)                                                                                                                              cTCCACTACAAACCGACCTAAGTCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTAACACCCATTCCGGTCAATACGGCAAATCGGt                   <  5:36017/149‑1 (MQ=255)                                                                                                                              cTCCACTACAAACCGACCTAAGTCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTAACACCCATTCCGGTCAATACGGCAAATCGGt                   <  1:362937/149‑1 (MQ=255)                                                                                                                              cTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTAACACCCATTCCGGTCAATACGGCAAATCGGt                   <  8:87157/149‑1 (MQ=255)                                                                                                                                 cACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTAACACCCATTCCGGTCAATACGGCAAATCGGTCAc                >  1:336300/1‑149 (MQ=255)                                                                                                                                      cAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCaa                                                                                         <  2:120872/71‑1 (MQ=255)                                                                                                                                      cAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCaa                                                                                         >  1:120872/1‑71 (MQ=255)                                                                                                                                      cAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTAACACCCATTCCGGTCAATACGGCAAATCGGTCACCCt             <  2:168372/147‑1 (MQ=255)                                                                                                                                       aaaCCGACCTAAGCCCGCAACACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTAACACCCATTCCGGTCAATACGGCAAATCGGTCACCCTCCa          <  6:277954/149‑1 (MQ=255)                                                                                                                                               cTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTAACACCCATTCCGGTCAATACGGCAAATCGGTCACCCTCCAGCCGCGtt  >  5:228837/1‑149 (MQ=255)                                                                                                                                                 aaGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTAACACCCATTCCGGTCAATACGGCAAATCGGTCACCCTCCAGCCGCGt   >  1:140174/1‑146 (MQ=255)                                                                                                                                                 aaGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTAACACCCATTCCGGTCAATACGGCAAATCGGTCACCCTCCAGCCGCGt   <  2:140174/146‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 TGCGCGTGGGTACAAGGGCACCTCTAAAAACGTATCGGGAGGCAGCAACTTCACAGTGACCAAACAGTCTTGTTCGGACTCTTGCGACAAATCCTCACGGGTCATTTCCTGTCGTGGCGTGTAAACTCCACTACAAACCGACCTAAGCCCGCAAAACCCTGCTCCACA‑CCTGGTCTTTAGAGGTGTACCTGCGCCCGGCTGCAAGCCTACGGCGAACGTAAAACCGCTCGTATCTTGACCGTAACACCCATTCCGGTCAATACGGCAAATCGGTCACCCTCCAGCCGCGTT  >  NC_002947/4747817‑4748107

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