Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,949,602	G→C	6.7%	P479A (CCC→GCC) ‡	PP_3486 ←	cytochrome c

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,949,602	0	G	C	6.7%	184.4 / 2.1	57	P479A (CCC→GCC) ‡	PP_3486	cytochrome c
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base C (2/2);  total (26/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 = 1.00e+00

CCATTGACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTC  >  NC_002947/3949454‑3949742                                                                                                                                                     |                                                                                                                                             ccATTGACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCggg                                                                                                                                              <  6:351395/149‑1 (MQ=255) ccATTGACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCggg                                                                                                                                              <  6:43023/149‑1 (MQ=255)       aCCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGgc                                                                                                                                        <  3:30415/149‑1 (MQ=255)           gcgcCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGgccacc                                                                                                                                    >  7:164764/1‑149 (MQ=255)                gccgGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCTGCCCGGgccacgcttta                                                                                                                               <  4:240206/149‑7 (MQ=255)                   ggCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCTGCCCGGgccacgctttac                                                                                                                              <  8:61470/147‑8 (MQ=255)                   ggCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCTGCCCGGgccacgctttac                                                                                                                              >  7:61470/1‑140 (MQ=255)                       caTCGTCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACTATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCTCCACATgg                                                                                                                        <  6:260102/149‑1 (MQ=255)                        aTCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGgccac                                                                                                                                     >  6:151314/1‑135 (MQ=255)                        aTCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGgc                                                                                                                                        >  5:56885/1‑132 (MQ=255)                        aTCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGgc                                                                                                                                        <  6:56885/132‑1 (MQ=255)                        aTCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCACCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGgccac                                                                                                                                     <  5:151314/135‑1 (MQ=255)                         tCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATggg                                                                                                                       >  6:224004/1‑148 (MQ=255)                                aGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGgccac                                                                                                                                     >  4:278729/1‑127 (MQ=255)                                aGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGgccac                                                                                                                                     <  3:278729/127‑1 (MQ=255)                                           gTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTggcg                                                                                                    >  6:190558/1‑149 (MQ=255)                                            tCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCTCTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCt                                                                                                        <  5:251009/144‑1 (MQ=255)                                            tCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCt                                                                                                        >  6:251009/1‑144 (MQ=255)                                            tCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTggcgg                                                                                                   >  5:216965/1‑149 (MQ=255)                                              aGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTggc                                                                                                     >  8:33210/1‑145 (MQ=255)                                              aGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTggc                                                                                                     <  7:33210/145‑1 (MQ=255)                                                 aCCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGc                                                                                                           >  6:38902/1‑136 (MQ=255)                                                 aCCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGc                                                                                                           <  5:38902/136‑1 (MQ=255)                                                 aCCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGa                                                                                                <  8:38710/147‑1 (MQ=255)                                                   cAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATg                                                                                                                         <  4:315462/120‑1 (MQ=255)                                                   cAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATg                                                                                                                         >  3:315462/1‑120 (MQ=255)                                                       aaCAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGc                                                                                         >  2:52910/1‑148 (MQ=255)                                                         cAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGc                                                                                                                      >  6:186570/1‑117 (MQ=255)                                                         cAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGc                                                                                                                      <  5:186570/117‑1 (MQ=255)                                                          aGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCaaa                                                                                     <  5:224004/149‑1 (MQ=255)                                                                  ccACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTg                                                                             <  6:216965/149‑1 (MQ=255)                                                                   cacaCGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGc                                                                            >  7:85450/1‑149 (MQ=255)                                                                            cGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCa                                                                    <  5:67747/148‑1 (MQ=255)                                                                             gCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAgg                                                                  >  3:80780/1‑149 (MQ=255)                                                                                   cAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTg                                                            <  8:164764/149‑1 (MQ=255)                                                                                    aaCCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCggg                                                                        <  6:126337/136‑1 (MQ=255)                                                                                    aaCCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCggg                                                                        >  5:126337/1‑136 (MQ=255)                                                                                    aaCCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTgg                                                           <  8:85450/149‑1 (MQ=255)                                                                                    aaCCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTgg                                                           >  5:22908/1‑149 (MQ=255)                                                                                     aCCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTggg                                                          <  1:52910/149‑1 (MQ=255)                                                                                       cTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCa                                                                                >  6:226152/1‑125 (MQ=255)                                                                                       cTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCa                                                                                <  5:226152/125‑1 (MQ=255)                                                                                           gCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATc                                                    <  1:306981/149‑1 (MQ=255)                                                                                            cTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGc                                                         >  1:233218/1‑143 (MQ=255)                                                                                            cTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGc                                                         <  2:233218/143‑1 (MQ=255)                                                                                                     tGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTg                                                            >  1:69969/1‑131 (MQ=255)                                                                                                     tGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTg                                                            <  2:69969/131‑1 (MQ=255)                                                                                                     tGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCt                                          <  4:80780/149‑1 (MQ=255)                                                                                                         cGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGa                                      <  4:40697/149‑1 (MQ=255)                                                                                                                                        aaGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATc       <  6:22908/149‑1 (MQ=255)                                                                                                                                          gCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATccc     >  2:34754/1‑149 (MQ=255)                                                                                                                                           ggTGGCCCTGCCCGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCg    >  6:137420/2‑149 (MQ=255)                                                                                                                                             tGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTc  <  5:190558/149‑1 (MQ=255)                                                                                                                                                 ccGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCa                                               >  7:184020/1‑100 (MQ=255)                                                                                                                                                 ccGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCa                                               <  8:184020/100‑1 (MQ=255)                                                                                                                                                 ccGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCt             >  3:281937/1‑134 (MQ=255)                                                                                                                                                 ccGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCt             <  4:281937/134‑1 (MQ=255)                                                                                                                                                 ccGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTg            <  5:205910/135‑1 (MQ=255)                                                                                                                                                 ccGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTg            >  6:205910/1‑135 (MQ=255)                                                                                                                                                     |                                                                                                                                             CCATTGACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTC  >  NC_002947/3949454‑3949742

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