Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,879,891 C→A 10.9% intergenic (+75/+61) PP_1686 → / ← PP_1687 glutathione peroxidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,879,8910CA10.9% 136.7 / 7.6 55intergenic (+75/+61)PP_1686/PP_1687glutathione peroxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base A (4/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

CCAGCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCGGCAT  >  NC_002947/1879746‑1880032
                                                                                                                                                 |                                                                                                                                             
ccAGCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGt                                                                                                                                            >  1:354949/1‑149 (MQ=255)
 cAGCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGcc                                                                                                                                              <  3:107397/146‑1 (MQ=255)
 cAGCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCg                                                                                                                                             <  3:115269/147‑1 (MQ=255)
  atccTGACCAAAGCCGAGGATACCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGGCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCTTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTtc                                                                                                                                          <  6:259206/147‑1 (MQ=38)
   gCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATAGCGACTTCAGGAGGCCGTtct                                                                                                                                         >  5:348725/1‑149 (MQ=255)
    ccTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTAACGAAAAGGCCATTGACTCACAGCCGGGCCCCGGCTGCTCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGAGACTTAAGGAGGCCGTtctc                                                                                                                                        >  1:352482/1‑149 (MQ=38)
     cTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTtctac                                                                                                                                       <  7:110723/149‑3 (MQ=255)
       gACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCt                                                                                                                                      >  7:146532/1‑148 (MQ=255)
         ccAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATCGCGACTTCAGGAGGCCGTTCTCCTGaa                                                                                                                                   <  6:348725/149‑1 (MQ=255)
          cAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAg                                                                                                                                  >  8:309645/1‑149 (MQ=255)
            aaCCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGAACGGCCTCCTGaa                                                                                                                                   <  4:194863/146‑1 (MQ=255)
            aaCCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGAACGGCCTCCTGaa                                                                                                                                   >  3:194863/1‑146 (MQ=255)
                    gatCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgat                                                                                                                         >  7:111253/1‑148 (MQ=255)
                      tCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatc                                                                                                                        >  3:136838/1‑147 (MQ=255)
                      tCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatc                                                                                                                        >  6:213448/1‑147 (MQ=255)
                        ccGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatc                                                                                                                        >  6:48662/1‑145 (MQ=255)
                        ccGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGAACGGCCTCCTGAAGt                                                                                                                                 >  3:261619/1‑136 (MQ=255)
                        ccGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGAACGGCCTCCTGAAGt                                                                                                                                 <  4:261619/136‑1 (MQ=255)
                                    ctgctATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGa                                                                                                        <  2:354949/149‑1 (MQ=255)
                                       ctATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCGTGag                                                                                                                                   >  2:305279/1‑118 (MQ=255)
                                         aTCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaagatcaa                                                                                                                <  8:48850/136‑1 (MQ=255)
                                         aTCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaagatcaa                                                                                                                >  7:48850/1‑136 (MQ=255)
                                          tCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtggc                                                                                                  <  1:150209/149‑1 (MQ=255)
                                             aaaaGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtggctgg                                                                                               >  4:137915/1‑149 (MQ=255)
                                               aaGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGAAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTggctggct                                                                                             >  7:278205/1‑149 (MQ=255)
                                                 ggCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGc                                                                                           >  7:337597/1‑149 (MQ=255)
                                                  gCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCg                                                                                          <  5:19883/149‑1 (MQ=255)
                                                    cATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCggg                                                                                        <  5:213448/149‑1 (MQ=255)
                                                              caGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCAcc                                                                              <  7:309645/149‑1 (MQ=255)
                                                                        cAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCgg                                                                                         <  7:199762/128‑1 (MQ=255)
                                                                        cAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCgg                                                                                         >  8:199762/1‑128 (MQ=255)
                                                                        cAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCAc                                                                    <  4:136838/149‑1 (MQ=255)
                                                                           gCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACgg                                                                  >  3:97668/1‑148 (MQ=255)
                                                                           gCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGt                                                                 <  1:4173/149‑1 (MQ=255)
                                                                             tGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTaa                                                               >  5:342649/1‑149 (MQ=255)
                                                                             tGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTaa                                                               >  5:298994/1‑149 (MQ=255)
                                                                               cgcTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAAt                                                              <  6:258223/148‑1 (MQ=255)
                                                                                  ttGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCATGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTTACCTTGCATGCACGGtgaagtgc                                                          <  3:137915/149‑8 (MQ=39)
                                                                                           cgTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGc                                                 <  6:123544/149‑1 (MQ=255)
                                                                                           cgTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGc                                                 <  8:337597/149‑1 (MQ=255)
                                                                                            gTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGca                                                <  5:48662/149‑1 (MQ=255)
                                                                                                 tttAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcaccac                                            >  5:172855/1‑148 (MQ=255)
                                                                                                    aaaGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGcca                                                         >  3:209264/1‑132 (MQ=255)
                                                                                                    aaaGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGcca                                                         <  4:209264/132‑1 (MQ=255)
                                                                                                        cgcgAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCAccggcc                                    <  5:330751/149‑1 (MQ=255)
                                                                                                             ggCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGacc                               <  8:341190/149‑1 (MQ=255)
                                                                                                                 gTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcaccaccac                                         >  3:291306/1‑135 (MQ=255)
                                                                                                                 gTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcaccaccac                                         <  4:291306/135‑1 (MQ=255)
                                                                                                                 gTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGAccacca                           <  8:172234/149‑1 (MQ=255)
                                                                                                                           cGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGcacacacac                 <  8:111253/149‑1 (MQ=255)
                                                                                                                              aGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACaaa              >  1:239632/1‑149 (MQ=255)
                                                                                                                                  cgcaaCTTCAGGAGAACGGCCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGc           >  1:226454/5‑148 (MQ=255)
                                                                                                                                      aCTTCATGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCg      <  8:146532/149‑1 (MQ=255)
                                                                                                                                       cTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACAACCAGCACACACACAAATGCCACCgg     >  8:156668/1‑149 (MQ=255)
                                                                                                                                        ttCAGGAGAACGGCCTCCTGAAGTCGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCgg     >  6:115894/1‑148 (MQ=255)
                                                                                                                                          cAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCCCCGGCAt  >  8:313488/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CCAGCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCGGCAT  >  NC_002947/1879746‑1880032

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