Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,913,888	T→C	9.4%	T271A (ACC→GCC)	ccmF ←	holocytochrome c synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,913,888	0	T	C	9.4%	112.3 / 2.1	42	T271A (ACC→GCC)	ccmF	holocytochrome c synthetase
Reads supporting (aligned to +/- strand):  ref base T (21/17);  new base C (2/2);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGCGACGGGTCGGCGGCAAACGCGTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCGTAGTAGGCCCACCACGACCCCAAGGTGATGCCAACACCCAGGAAAGCC  >  NC_002947/4913750‑4914036                                                                                                                                           |                                                                                                                                                     cccGCGACGGGTCGGCGGCAAACGCGTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGc                                                                                                                                            >  2:156384/1‑149 (MQ=255)        cGGGTCGGCGGCAAACGCGTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCa                                                                                                                                              >  1:250098/1‑140 (MQ=255)        cGGGTCGGCGGCAAACGCGTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCa                                                                                                                                              <  2:250098/140‑1 (MQ=255)        cGGGTCGGCGGCAAACGCGTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGa                                                                                                                                     >  1:46297/1‑149 (MQ=255)              ggcggcAAACGCGTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTg                                                                                                                               >  2:195403/1‑149 (MQ=255)                  gcAAACGCGTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGt                                                                                                                           >  3:289389/1‑149 (MQ=255)                   cAAACGCGTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTg                                                                                                                          <  4:289389/149‑1 (MQ=255)                        gcgTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCccac                                                                                                                     <  1:156384/149‑1 (MQ=255)                                 gAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCCCGGTCACCGCCAGCGAGTGGATCa                                                                                                                                  <  4:183718/127‑1 (MQ=255)                                 gAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCCCGGTCACCGCCAGCGAGTGGATCa                                                                                                                                  >  3:183718/1‑127 (MQ=255)                                    gTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCTTTTCTCGGTCACCGCCAGCGAGTGTATCAGTGCGGTGCCCACCAGCCAAGGCa                                                                                                          <  8:44452/148‑1 (MQ=255)                                     tCAGCACTCCGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCAt                                                                                                         >  8:137209/1‑148 (MQ=255)                                           cTCCGGAACTGACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGGCAGCTGGACGGTCCAGCTCTTGAACAGCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAg                                                                                                  <  2:301725/149‑1 (MQ=255)                                              cGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCg                                                                                               >  8:99579/1‑149 (MQ=255)                                                      aCCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGTAATTGCCAGCAGCACGGTCCAGCTCTTGTACACCCCGCGTTTCTCGGCCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTcc                                                                                         <  6:187960/147‑1 (MQ=255)                                                      aCCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGTAATTGCCAGCAGCACGGTCCAGCTCTTGTACACCCCGCGTTTCTCGGCCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTcc                                                                                         >  5:187960/1‑147 (MQ=255)                                                             aGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTcc                                                                                <  2:46297/149‑1 (MQ=255)                                                                  cccAGCAGGCTCAGGGAGAATGCGGGAATTGCCAGCAGCACGGTCCTGCTCTTGAACACCCCGCGTTTCTCGGTGACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGaa                                                                           <  5:112667/149‑1 (MQ=255)                                                                     agcagGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATg                                                                                                        <  7:157964/117‑1 (MQ=255)                                                                     agcagGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCCCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATg                                                                                                        >  8:157964/1‑117 (MQ=255)                                                                     agcagGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGCCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAacca                                                                        >  8:197351/1‑149 (MQ=255)                                                                            tCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGCACCACCAGccc                                                                 >  2:61185/1‑148 (MQ=255)                                                                              aGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCAcc                                                               >  1:325834/1‑149 (MQ=255)                                                                                 gagaATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAacca                                                                        <  5:319667/137‑1 (MQ=255)                                                                                 gagaATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAacca                                                                        >  6:319667/1‑137 (MQ=255)                                                                                          gCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTc                                                   >  4:283644/1‑149 (MQ=255)                                                                                           cAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGCCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCg                                                  <  7:197351/149‑1 (MQ=255)                                                                                            aaTTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCTCAGAACCACCAGCCACcccagcccactacg                                                  >  1:130868/1‑143 (MQ=255)                                                                                             aTTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCgta                                                <  7:99579/149‑1 (MQ=255)                                                                                                 ccAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCgtagta                                             <  3:283644/148‑1 (MQ=255)                                                                                                 ccAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCGtagtag                                            <  7:137209/149‑1 (MQ=255)                                                                                                    gcagcaCGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCGTAGTAGGcc                                         >  7:159059/1‑149 (MQ=255)                                                                                                                aGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCGTAGTAGGCCCACCACGAccc                              >  2:246051/1‑148 (MQ=255)                                                                                                                 gCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGt                                                                                              <  3:6530/83‑1 (MQ=255)                                                                                                                 gCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGt                                                                                              >  4:6530/1‑83 (MQ=255)                                                                                                                   tctTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCGTAGTAGGCCCACCACGACCCCAAg                          >  1:171865/1‑149 (MQ=255)                                                                                                                       gAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCGTAGTAGGCCCACCACGACCCCAAGGt                        <  1:246051/147‑1 (MQ=255)                                                                                                                       gAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCGTAGTAGGCCCACCACGACCCCAAGGTGa                      >  6:283181/1‑149 (MQ=255)                                                                                                                           acCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGAGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCGTAGTAGGCCCACCACGACCCCAAGGTGATGc                   >  1:63016/1‑148 (MQ=255)                                                                                                                              ccGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCGTAGTAgg                                           >  6:55708/1‑121 (MQ=255)                                                                                                                              ccGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCGTAGTAgg                                           <  5:55708/121‑1 (MQ=255)                                                                                                                              ccGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCGTAGTAGGCCCACCACGACCCCAAGGTGATGCCAac               <  1:195403/149‑1 (MQ=255)                                                                                                                                           tCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCGTAGTAGGCCCACCACGACCCCAAGGTGATGCCAACACCCAGGAAAGcc  <  8:159059/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                     CCCGCGACGGGTCGGCGGCAAACGCGTGCACCGAGGTCAGCACTCCGGAACGCACCAGGAAGGTACCCAGCAGGCTCAGGGAGAATGCGGCAATTGCCAGCAGCACGGTCCAGCTCTTGAACACCCCGCGTTTCTCGGTCACCGCCAGCGAGTGGATCAGTGCGGTGCCCACCAGCCAAGGCATGAACGAGGCGTTTTCCACCGGGTCCCAGAACCACCAGCCACCCCAGCCCAGTTCGTAGTAGGCCCACCACGACCCCAAGGTGATGCCAACACCCAGGAAAGCC  >  NC_002947/4913750‑4914036

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