Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,949,604	A→C	6.9%	L478R (CTG→CGG)	PP_3486 ←	cytochrome c

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,949,604	0	A	C	6.9%	178.3 / 2.7	57	L478R (CTG→CGG)	PP_3486	cytochrome c
Reads supporting (aligned to +/- strand):  ref base A (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 = 9.97e-01

ACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTC  >  NC_002947/3949460‑3949742                                                                                                                                                 |                                                                                                                                           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                                                                                                                                     <  3:278729/127‑1 (MQ=255)                          aGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGgccac                                                                                                                                     >  4:278729/1‑127 (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                                                                                                     <  7:33210/145‑1 (MQ=255)                                        aGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTggc                                                                                                     >  8:33210/1‑145 (MQ=255)                                           aCCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGc                                                                                                           <  5:38902/136‑1 (MQ=255)                                           aCCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGc                                                                                                           >  6:38902/1‑136 (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                                                                                                                      <  5:186570/117‑1 (MQ=255)                                                   cAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGc                                                                                                                      >  6:186570/1‑117 (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                                                                        >  5:126337/1‑136 (MQ=255)                                                                              aaCCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCggg                                                                        <  6:126337/136‑1 (MQ=255)                                                                              aaCCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTgg                                                           <  8:85450/149‑1 (MQ=255)                                                                              aaCCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTgg                                                           >  5:22908/1‑149 (MQ=255)                                                                               aCCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTggg                                                          <  1:52910/149‑1 (MQ=255)                                                                                 cTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCa                                                                                <  5:226152/125‑1 (MQ=255)                                                                                 cTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCa                                                                                >  6:226152/1‑125 (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             <  4:281937/134‑1 (MQ=255)                                                                                                                                           ccGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCt             >  3:281937/1‑134 (MQ=255)                                                                                                                                           ccGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTg            <  5:205910/135‑1 (MQ=255)                                                                                                                                           ccGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTg            >  6:205910/1‑135 (MQ=255)                                                                                                                                                 |                                                                                                                                           ACCCGCGCCGCCGGCACATCGGCCCAGCTCACGCTGGTCAGCACCAGCAACAGGCACAACCACACGATCCGCATGACAACCTCCGCTCAAGGGCTGGCCGGCACTGGCCGGTAGGGGCTGACCGAAGTAAAGCGTGGCCCGGGCAGGGCCACCGCCACCACATGGGCACCTGCCATGCCCTGGCGGCGACCAGGGCCAAAGCCCAGTGCCGGGGTCAGGCCGGTGGGCACATCGTGCAGGCCTTCGAGCGCAGCAATCAGTTGCTCGCGGCTGGCATCCCGTC  >  NC_002947/3949460‑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: ‑