Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 240,531 G→C 23.1% A276P (GCC→CCC)  PP_0188 → uroporphyrin‑III C‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947240,5310GC23.1% 42.2 / 13.2 26A276P (GCC→CCC) PP_0188uroporphyrin‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (11/9);  new base C (3/3);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.58e-01

CCGAGTTCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATGCA  >  NC_002947/240388‑240674
                                                                                                                                               |                                                                                                                                               
ccGAGTTCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGCCCggcc                                                                                                                                           >  2:232220/1‑148 (MQ=255)
      tCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGCCCGGCCAGCAAc                                                                                                                                     <  1:232220/150‑1 (MQ=255)
                              gcgcCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAACGTGCGTCCGTTGCTGGCCg                                                                                                                                              <  2:85536/116‑1 (MQ=255)
                              gcgcCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAACGTGCGTCCGTTGCTGGCCg                                                                                                                                              >  1:85536/1‑116 (MQ=255)
                                    tGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTTTCCCTGAGCCTGa                                                                                                       <  2:550310/150‑1 (MQ=255)
                                       ccgccgACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGGGCCTgctcg                                                                                                    >  1:31724/1‑146 (MQ=255)
                                                                    tCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAATCGTGCGTCCGGTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGCGGGCGGCGCTCAAt                                                                       >  2:208515/1‑150 (MQ=255)
                                                                    tCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAAt                                                                       >  1:132244/1‑150 (MQ=255)
                                                                       cAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGt                                                                    <  2:132244/150‑1 (MQ=255)
                                                                            gtggGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCTGCGCTCAATGGTGATGc                                                               <  1:26817/150‑1 (MQ=255)
                                                                                gCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGTGCCTGGCCCTGAGCCTGACTAt                                                                                                   <  1:349163/110‑1 (MQ=255)
                                                                                gCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGTGCCTGGCCCTGAGCCTGACTAt                                                                                                   >  2:349163/1‑110 (MQ=255)
                                                                                                             cTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTggg                                                                                  <  1:525767/98‑1 (MQ=255)
                                                                                                             cTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTggg                                                                                  >  2:525767/1‑98 (MQ=255)
                                                                                                                        aTGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAg                                                           <  1:672873/110‑1 (MQ=255)
                                                                                                                        aTGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAg                                                           >  2:672873/1‑110 (MQ=255)
                                                                                                                              aCGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCa             >  2:390988/1‑150 (MQ=255)
                                                                                                                               cGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCaa            <  2:694397/150‑1 (MQ=255)
                                                                                                                                      ccGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAAt     >  2:635302/1‑150 (MQ=255)
                                                                                                                                        gTTGCTGCCCGGCCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATgg                                                                     <  1:573340/84‑1 (MQ=255)
                                                                                                                                        gTTGCTGCCCGGCCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATgg                                                                     >  2:573340/1‑84 (MQ=255)
                                                                                                                                         ttGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATACa  >  1:184960/1‑150 (MQ=255)
                                                                                                                                           gCTGCCCGGCCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCaa            >  2:50855/1‑138 (MQ=255)
                                                                                                                                           gCTGCCCGGCCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCaa            <  1:50855/138‑1 (MQ=255)
                                                                                                                                             tGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATg    <  1:800394/144‑1 (MQ=255)
                                                                                                                                             tGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATg    >  2:800394/1‑144 (MQ=255)
                                                                                                                                               |                                                                                                                                               
CCGAGTTCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATGCA  >  NC_002947/240388‑240674

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: