Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 239,636 T→C 10.6% S239S (AGT→AGC) ‡ hemD → uroporphyrinogen‑III synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947239,6360TC10.6% 118.1 / 5.9 47S239S (AGT→AGC) ‡hemDuroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand):  ref base T (21/21);  new base C (3/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

GTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTCTTGTCCAATAACGATAAACCGTCGGCGCAGGC  >  NC_002947/239488‑239779
                                                                                                                                                        |                                                                                                                                               
gTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGt                                                                                                                                                 <  8:24309/149‑1 (MQ=255)
    gCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCa                                                                                                                                             >  6:353742/1‑149 (MQ=255)
    gCAGGGTTTTCAACATTTGCAGCAGATGGCCCGTCCCGATTGGCTGCGCTTGGCCCAGATGACGCAGTTCCTGCCCAGCCCGCGGCTCGCCGAGCAGGCCAAGGCCCCCGGGCC‑‑CCAACAGGAT‑GT‑GGATTGACGTGCCGCCAGTGcc                                                                                                                                              >  8:263639/1‑148 (MQ=255)
        ggTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGGCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCACTGGCGCCACGGc                                                                                                                                         <  1:255158/149‑1 (MQ=255)
        ggTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCACTGGCGCCACGGc                                                                                                                                         <  1:259881/149‑1 (MQ=255)
         gTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGcc                                                                                                                                        <  7:263639/149‑1 (MQ=255)
             tGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGc                                                                                                                                    >  2:356317/1‑149 (MQ=255)
             tGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGc                                                                                                                                    >  2:257473/1‑149 (MQ=255)
               aaCATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTg                                                                                                                                  <  1:356317/149‑1 (MQ=255)
                  aTTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACCGCCTTGCTgg                                                                                                                                 >  4:305018/1‑147 (MQ=255)
                       cagcagATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCACCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGt                                                                                                                          <  4:7883/149‑1 (MQ=255)
                         gcagATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGc                                                                                                                        <  4:173542/149‑1 (MQ=255)
                          cagATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCa                                                                                                                       <  5:137176/149‑1 (MQ=255)
                           agATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAg                                                                                                                      <  7:284436/149‑1 (MQ=255)
                            gATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGc                                                                                                                     >  2:72324/1‑149 (MQ=255)
                                gCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTa                                                                                                                  <  8:128458/148‑1 (MQ=255)
                                gCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTa                                                                                                                  <  1:255959/148‑1 (MQ=255)
                                gCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTa                                                                                                                  <  5:34649/148‑1 (MQ=255)
                                 ccGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAgc                                                                                                                >  4:230426/1‑149 (MQ=255)
                                               cTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCa                                                                                                                                             >  2:88074/1‑106 (MQ=255)
                                                         cccAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGgcgc                                                                                         >  4:94566/1‑148 (MQ=255)
                                                               tGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCt                                                                                   <  3:230426/148‑1 (MQ=255)
                                                                gCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGATCAGGCCAAGGCCGCCGGGCCTACTATCATCTTAATCCAATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTaa                                                                                      <  1:294181/148‑1 (MQ=255)
                                                                    cTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGACAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGc                                                                             <  2:356689/149‑1 (MQ=255)
                                                                    cTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGACAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGc                                                                             >  1:356689/1‑149 (MQ=255)
                                                                     tGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCg                                                                            <  7:249431/149‑1 (MQ=255)
                                                                       ttCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCg                                                                            <  5:353742/147‑1 (MQ=255)
                                                                           tGCCCAGCCCGCGGGTCGCCGAGCAGGCGAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTgcgc                                                                      >  1:201998/1‑149 (MQ=255)
                                                                             cccAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGccc                                                                    >  4:263358/1‑149 (MQ=255)
                                                                                    cgcgGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGc                                                              >  3:81430/1‑148 (MQ=255)
                                                                                          tCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAggat                                                       <  3:305018/149‑1 (MQ=255)
                                                                                              cGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAggatggat                                                   <  3:94566/149‑1 (MQ=255)
                                                                                                                        ccAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTCTTGTcc                            >  7:178288/1‑148 (MQ=255)
                                                                                                                                       gATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTCTTGTc                             >  5:157802/1‑134 (MQ=255)
                                                                                                                                       gATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTCTTGTc                             <  6:157802/134‑1 (MQ=255)
                                                                                                                                       gATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTCTTGTCCAATAACGATAAAcc              >  4:308483/1‑149 (MQ=255)
                                                                                                                                         ttGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCCCCGGCCGCTTCAGGCGCGAAGCGGCCTGCGGGGCGCCCCCAGGCCAAGGATGGCTACGGGAGCGAGCCTGGCGTGTCCAATAACGATCACCCGt            >  5:47190/1‑149 (MQ=255)
                                                                                                                                         ttGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGGCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTCTTGTCGAATAACGATAAACCGt            >  5:52416/1‑149 (MQ=255)
                                                                                                                                          ggCCGTGGCACTGGCGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTCTTGTCCAATAACGATAAACCGTc           >  5:337804/2‑149 (MQ=255)
                                                                                                                                           gCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTCTTGTCCAATAACGATAAACCGTc           <  2:201998/148‑1 (MQ=255)
                                                                                                                                            ccGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTCTTGTCCAATAACGATAAACCGTCg          <  1:257473/148‑1 (MQ=255)
                                                                                                                                            ccGTGGCACTGGCGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTCTTGTCCAATAACGATAAACCGTCgg         >  3:190261/1‑149 (MQ=255)
                                                                                                                                            ccGTGGCACTGGCGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTCTTGTCCAATAACGATAAACCGTCgg         >  1:42720/1‑149 (MQ=255)
                                                                                                                                              gTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTc                                  >  2:109769/1‑122 (MQ=255)
                                                                                                                                              gTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTc                                  <  1:109769/122‑1 (MQ=255)
                                                                                                                                                   gcCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTTAGCTGAATGCGGTGCGccccc                                                                  <  8:285259/85‑1 (MQ=255)
                                                                                                                                                   gcCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGccccc                                                                  >  7:285259/1‑85 (MQ=255)
                                                                                                                                                   gcCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTCTTGTCCAATAACGATAAACCGTCGGCGCAGgc  >  3:347195/1‑149 (MQ=255)
                                                                                                                                                        |                                                                                                                                               
GTGGGCAGGGTTTTGAACATTTGCAGCAGATGGCCGGTGCCGATTGGCTGCGGTTGGCCCAGTTGCCGCTGTTCGTGCCCAGCCCGCGGGTCGCCGAGCAGGCCAAGGCCGCCGGGGC‑‑CCAACAGGTT‑GT‑GGATTGCCGTGGCGCCAGTGCCACGGCCTTGCTGGCAGCCGTGCAGCGTAGCGCTGCACCTGCCTCTTAAGGCGCTAAGCTGAATGCGGTGCGCCCCCATGCAAAGGATGGATACGTGAGCGAGACTGTCTTGTCCAATAACGATAAACCGTCGGCGCAGGC  >  NC_002947/239488‑239779

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