Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,016,716 T→G 9.4% intergenic (‑32/+195) 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_0029476,016,7160TG 9.4% 155.3 / 8.4 54intergenic (‑32/+195)PP_5728/PP_5267rhs repeat‑associated core domain‑containing protein/cytochrome c5
Reads supporting (aligned to +/- strand):  ref base T (26/22);  new base G (3/2);  total (29/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 = 6.21e-01

AGAGATAGCCATAAGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCAC  >  NC_002947/6016570‑6016860
                                                                                                                                                  |                                                                                                                                                
agagaTAGCCATAAGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAgtgt                                                                                                                                                >  5:246058/1‑149 (MQ=255)
    aTAGCCATAAGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGttctt                                                                                                                                            >  3:165601/1‑149 (MQ=255)
      aGCCATAAGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGttctt                                                                                                                                            >  7:27092/1‑147 (MQ=255)
      aGCCATAAGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGACAGTcttcttaa                                                                                                                                          >  7:188138/1‑147 (MQ=255)
             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                                                                                                                            >  6:189471/1‑132 (MQ=255)
                                cGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGGGTTCTGAAGAACACTGacgctt                                                                                                                            <  5:189471/137‑6 (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                                                                                  <  3:228229/149‑1 (MQ=255)
                                                              ggCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGc                                                                                  >  7:225188/1‑149 (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                >  3:111696/1‑149 (MQ=255)
                                                                                                                                ttCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCt                >  2:250523/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)
                                                                                                                                                  |                                                                                                                                                
AGAGATAGCCATAAGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCAC  >  NC_002947/6016570‑6016860

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