Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,682,573 G→C 7.5% M246I (ATG→ATC)  PP_3248 → Dyp‑type peroxidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,682,5730GC 7.5% 163.4 / 3.8 53M246I (ATG→ATC) PP_3248Dyp‑type peroxidase family protein
Reads supporting (aligned to +/- strand):  ref base G (24/25);  new base C (2/2);  total (26/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 = 8.95e-01

AAGCGCACCGCCCAGGAAAGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGCGAGACGGGTGTCGATCTAAGCCCCTTGCTGCGGGCCTGATCTCTAAACATCCAGATA  >  NC_002947/3682430‑3682717
                                                                                                                                                |                                                                                                                                                
aaGCGCACCGCCCAGGAAAGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCCACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTgg                                                                                                                                             >  4:277434/1‑149 (MQ=255)
   cgcACCGCCCAGGAAAGCTTTGAACCCGAGGCGTTCATGATCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCt                                                                                                                                          >  4:71514/1‑149 (MQ=255)
      aCCGCCCAGGAAAGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTgg                                                                                                                                        >  5:138465/1‑148 (MQ=255)
       ccGCCCAGGAAAGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGaa                                                                                                                                      <  5:140701/149‑1 (MQ=255)
          cccAGGAAAGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGAc                                                                                                                                   >  1:106678/1‑149 (MQ=255)
          cccAGGAAAGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGAc                                                                                                                                   >  3:158967/1‑149 (MQ=255)
             aGGAAAGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGc                                                                                                                                <  3:71514/149‑1 (MQ=255)
               gAAAGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGcat                                                                                                                              <  5:83877/149‑1 (MQ=255)
                aaaGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATCAGTGGCCTGGAAGa                                                                                                                                    <  2:91771/142‑1 (MQ=255)
                aaaGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATCAGTGGCCTGGAAGa                                                                                                                                    <  6:287079/142‑1 (MQ=255)
                aaaGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATCAGTGGCCTGGAAGa                                                                                                                                    >  1:91771/1‑142 (MQ=255)
                aaaGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATCAGTGGCCTGGAAGa                                                                                                                                    >  5:287079/1‑142 (MQ=255)
                 aaGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGcatca                                                                                                                            <  3:277434/149‑1 (MQ=255)
                  aGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGcatca                                                                                                                            >  4:71656/1‑148 (MQ=255)
                    cTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGTCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCatcatc                                                                                                                          >  8:182374/1‑148 (MQ=255)
                              gAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTa                                                                                                               >  6:352147/1‑149 (MQ=255)
                              gAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCAACATCGCCGGGCTCTa                                                                                                               >  6:347016/1‑149 (MQ=255)
                               aGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATTCCCAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTa                                                                                                               <  1:98273/149‑1 (MQ=255)
                                       aTGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCa                                                                                                       <  3:71656/148‑1 (MQ=255)
                                        tGGTCCGTCGCTCAGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTgg                                                                                                                                        <  6:22545/114‑1 (MQ=255)
                                        tGGTCCGTCGCTCAGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTgg                                                                                                                                        >  5:22545/1‑114 (MQ=255)
                                              gTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGccgcccg                                                                                               <  7:153387/149‑1 (MQ=255)
                                                 gCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTGGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTg                                                                                            <  2:106678/149‑1 (MQ=255)
                                                     ggTATCCTGGGTCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGAc                                                                                          <  7:182374/147‑1 (MQ=255)
                                                       tATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACgggg                                                                                      >  4:342009/1‑149 (MQ=255)
                                                            tGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGa                                                                                                                                       >  7:203042/1‑95 (MQ=255)
                                                            tGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGa                                                                                                                                       <  8:203042/95‑1 (MQ=255)
                                                               gCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGctacta                                                                              <  5:352147/149‑1 (MQ=255)
                                                               gCCGACCAACGCGGCGCCGGCCGGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGctacta                                                                              <  5:347016/149‑1 (MQ=255)
                                                                    ccAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGt                                                                         <  3:342009/149‑1 (MQ=255)
                                                                                 ggCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGa                                                                                                                        <  7:24965/89‑1 (MQ=255)
                                                                                 ggCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGa                                                                                                                        >  8:24965/1‑89 (MQ=255)
                                                                                    cTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGc                                                           >  1:309937/1‑147 (MQ=255)
                                                                                    cTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGCGGCTACTACTGGTGCCCGCCAATGAGc                                                           >  1:315076/1‑147 (MQ=255)
                                                                                    cTGGCTTTCGTCGCCCGGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGCga                                                         <  6:367970/149‑1 (MQ=255)
                                                                                       gCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCtt                                                                                                         <  3:140713/98‑1 (MQ=255)
                                                                                       gCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCtt                                                                                                         >  4:140713/1‑98 (MQ=255)
                                                                                       gCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGCGAGAc                                                      >  5:287370/1‑149 (MQ=255)
                                                                                         tttCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACgggggg                                                                                    <  6:21986/117‑1 (MQ=255)
                                                                                         tttCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACgggggg                                                                                    >  5:21986/1‑117 (MQ=255)
                                                                                          ttCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGCGAGACggg                                                   >  3:18239/1‑149 (MQ=255)
                                                                                           tcgtcgCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGCGAGACgg                                                    <  4:158967/147‑1 (MQ=255)
                                                                                                 cccTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGCGAGACGGGTGTc                                               >  5:99774/1‑146 (MQ=255)
                                                                                                 cccTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGCGAGACGGGTGTc                                               <  6:99774/146‑1 (MQ=255)
                                                                                                                       aTT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGCGAGACGGGTGTCGATCTAAGCCCCTTGCTGCGGGcca                      >  5:115460/1‑148 (MQ=255)
                                                                                                                       aTT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGCGAGACGGGTGTCGATCTAAGCCCCTTGCTGCGGGc                        <  6:115460/147‑1 (MQ=255)
                                                                                                                        tt‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTAAAGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGCGAGACGGGTGTCGATCTAAGCCCCTTGCTGCGGGCCt                      >  6:94913/1‑148 (MQ=255)
                                                                                                                             aGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGa                                                             >  2:4428/1‑105 (MQ=255)
                                                                                                                             aGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGa                                                             <  1:4428/105‑1 (MQ=255)
                                                                                                                              ggTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCACTGAGCGAGACGGGTGTCGATCTAAGCCCCTTGCTGCGGGCCTGAtctc                <  6:287370/149‑1 (MQ=255)
                                                                                                                               gTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGCGAGACGGGTGTCGATCTAAGCCCCTTGCTGCGGGCCTGATctct               <  4:18239/149‑1 (MQ=255)
                                                                                                                                  cAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGACGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGCGAGACGGGTGTCGATCTAAGCCCCTTGCTGCGGGCCTGATCTCTaa             <  2:38183/148‑1 (MQ=255)
                                                                                                                                              aTGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGCGAGACGGGTGTCGATCTAAGCCCCTTGCTGCGGGCCTGATCTCTAAACATCCAGATa  >  8:259028/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AAGCGCACCGCCCAGGAAAGCTTTGAACCCGAGGCGTTCATGGTCCGTCGCTCGGTATCCTGGGCCGACCAACGCGGCGCCGGCCTGGCTTTCGTCGCCCTGGGCAAGAGTTTCGAGGCATT‑CGAGGTGCAGCTGCGGCGCATGAGTGGCCTGGAAGACGGCATCATCGACGGGCTGTACCGCTTCAGCCGCCCGCTGACGGGGGGCTACTACTGGTGCCCGCCAATGAGCGAGACGGGTGTCGATCTAAGCCCCTTGCTGCGGGCCTGATCTCTAAACATCCAGATA  >  NC_002947/3682430‑3682717

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