Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,016,723	C→A	9.8%	intergenic (‑39/+188)	PP_5728 ← / ← PP_5267	rhs repeat‑associated core domain‑containing protein/cytochrome c5

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,016,723	0	C	A	9.8%	157.7 / 6.6	51	intergenic (‑39/+188)	PP_5728/PP_5267	rhs repeat‑associated core domain‑containing protein/cytochrome c5
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base A (3/2);  total (26/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.66e-01

AGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAA  >  NC_002947/6016583‑6016869                                                                                                                                             |                                                                                                                                                   aGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTATCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCt                                                                                                                                            <  5:3637/149‑1 (MQ=255)  ggCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCt                                                                                                                                            >  2:21305/1‑148 (MQ=255)       aTAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGTAGCGTCAGTGTTCTTCAGAACCCTGaaaaa                                                                                                                                      <  1:24059/149‑1 (MQ=255)           gCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATTGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGa                                                                                                                                  <  2:52689/149‑1 (MQ=255)            ccAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGCGTTCTTCAGAACCCTGAAAAATAGaa                                                                                                                                 <  5:243952/149‑1 (MQ=255)                 tGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGt                                                                                                                            <  8:27092/149‑1 (MQ=255)                  gCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGt                                                                                                                            >  4:228229/1‑148 (MQ=255)                    cGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGGGTTCTGAAGAACACTGacgctt                                                                                                                                     <  5:189471/137‑6 (MQ=255)                    cGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGGGTTCTGAAGAACACTGacgctt                                                                                                                                     >  6:189471/1‑132 (MQ=255)                                   tggATCAGCTTTTGCGTCTGTGTTCGTCCGGCTGCATGTTTCTCGGGACAGTTATCGCCCAGTTCCATTGTTATCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACACTGTAAAATAGATCTTTTAAAAATGGTATTTCTCAt                                                                                                          <  8:188138/148‑1 (MQ=255)                                     gATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATgg                                                                                                        <  5:75034/149‑1 (MQ=255)                                       tCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGGCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGaa                                                                                                      <  3:45221/149‑1 (MQ=255)                                          gcgtttGCTGCTGTGTTCGACCTTCTTCATGTTGCTCGGGACAGTTATCGCACAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATATAACTTGTAAAAATGCTAGTTCTCATGGAATCa                                                                                                   <  4:162413/146‑1 (MQ=255)                                                  ggCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGc                                                                                           >  7:225188/1‑149 (MQ=255)                                                  ggCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGc                                                                                           <  3:228229/149‑1 (MQ=255)                                                  ggCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGc                                                                                           >  4:196728/1‑149 (MQ=255)                                                   gCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGGGTTCTGAAGAACACTGac                                                                                                                                         <  3:196975/102‑2 (MQ=255)                                                   gCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGGGTTCTGAAGAACACTGac                                                                                                                                         >  4:196975/1‑101 (MQ=255)                                                         ttCGACCGACTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTg                                                                                    >  4:179154/1‑149 (MQ=255)                                                                 gCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCt                                                                            >  1:204545/1‑149 (MQ=255)                                                                  cTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTg                                                                           <  5:246161/149‑1 (MQ=255)                                                                    gCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACGCTGAAAAATAGaa                                                                                                                                 <  6:24918/93‑1 (MQ=255)                                                                    gCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGaa                                                                                                                                 >  5:24918/1‑93 (MQ=255)                                                                    gCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTg                                                                           >  8:134534/1‑147 (MQ=255)                                                                     cATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGga                                                                         <  6:246058/148‑1 (MQ=255)                                                                      aTGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGgag                                                                        >  4:42715/1‑148 (MQ=255)                                                                        gTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTCGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACa                                                                     >  6:227720/1‑149 (MQ=255)                                                                         ttGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACa                                                                     <  4:165601/148‑1 (MQ=255)                                                                         ttGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTCGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACAt                                                                    <  5:227720/149‑1 (MQ=255)                                                                             tCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCAt                                                                >  4:109892/1‑149 (MQ=255)                                                                               gggACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCAt                                                                <  3:237308/147‑1 (MQ=255)                                                                                 gACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGa                                                             >  6:58728/1‑148 (MQ=255)                                                                                      ttATCGCCCAGATACATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCAATGGACATgg                                                        >  2:71000/1‑148 (MQ=255)                                                                                        aTCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGctc                                                     >  5:172747/1‑149 (MQ=255)                                                                                           gCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTa                                                   <  8:225188/148‑1 (MQ=255)                                                                                                  tCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGggttg                                           <  3:42715/149‑1 (MQ=255)                                                                                                   ccATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGggttgg                                          >  4:68892/1‑149 (MQ=255)                                                                                                              caccacTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCt                               >  3:146018/1‑149 (MQ=255)                                                                                                               accacTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTg                              >  2:185440/1‑149 (MQ=255)                                                                                                               accacTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAACGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTg                              >  7:28203/1‑149 (MQ=255)                                                                                                                 cacTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGa                            <  3:196728/149‑1 (MQ=255)                                                                                                                    ttCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCt                         >  2:250523/1‑149 (MQ=255)                                                                                                                    ttCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCt                         >  3:111696/1‑149 (MQ=255)                                                                                                                      cAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCAt                                                                <  6:53863/108‑1 (MQ=255)                                                                                                                      cAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCAt                                                                >  5:53863/1‑108 (MQ=255)                                                                                                                            gAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCggg                 >  3:5333/1‑149 (MQ=255)                                                                                                                              aGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGca               <  3:68892/149‑1 (MQ=255)                                                                                                                                 ttCAGGGTTCTGAAGAACACTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCAATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGcacca            >  6:161463/2‑149 (MQ=255)                                                                                                                                  tCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGcaccac           <  4:111696/149‑1 (MQ=255)                                                                                                                                    aGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGca               >  3:201114/1‑143 (MQ=255)                                                                                                                                    aGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGca               <  4:201114/143‑1 (MQ=255)                                                                                                                                          cttcAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCtgtg                                 <  7:242612/119‑1 (MQ=255)                                                                                                                                          cttcAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCtgtg                                 >  8:242612/1‑119 (MQ=255)                                                                                                                                           ttcAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACaaa  <  2:204545/149‑1 (MQ=255)                                                                                                                                           ttcAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACa    >  4:214611/1‑147 (MQ=255)                                                                                                                                             |                                                                                                                                                   AGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAA  >  NC_002947/6016583‑6016869

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