Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,016,729	C→A	10.0%	intergenic (‑45/+182)	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,729	0	C	A	10.0%	163.3 / 5.4	56	intergenic (‑45/+182)	PP_5728/PP_5267	rhs repeat‑associated core domain‑containing protein/cytochrome c5
Reads supporting (aligned to +/- strand):  ref base C (26/23);  new base A (3/3);  total (29/27)
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.88e-01

AGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCC  >  NC_002947/6016583‑6016872                                                                                                                                                   |                                                                                                                                                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                                                                   >  5:53863/1‑108 (MQ=255)                                                                                                                      cAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCAt                                                                   <  6:53863/108‑1 (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)                                                                                                                                              aGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGcc  >  6:17691/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                AGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCC  >  NC_002947/6016583‑6016872

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