Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,382,118 C→T 8.9% V42M (GTG→ATG)  cobA ← uroporphyrinogen‑III C‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,382,1180CT 8.9% 124.9 / 3.0 45V42M (GTG→ATG) cobAuroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (19/22);  new base T (2/2);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

CAAGGGTCGCCGCCCTTGAGCCGCACCACGCAACGGCCCTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCACCAGGGCCGGCACCGACCAGCCATACTTTTGCATTCATCAGCGTTTCCTCATCGGCAAGG  >  NC_002947/2381970‑2382263
                                                                                                                                                    |                                                                                                                                                 
cAAGGGTCGCCGCCCTTGAGCCGCACCACGCAACGGCCCTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACt                                                                                                                                                  <  2:61449/150‑1 (MQ=255)
 aaGGGTCGCCGCCCTTGAGCCGCACCACGCAACGGCCCTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTg                                                                                                                                                 >  2:505621/1‑150 (MQ=255)
    ggTCGCCGCCCTTGAGCCGCACCACGCAACGGCCCTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCg                                                                                                                                              >  1:815140/1‑150 (MQ=255)
     gTCGCCGCCCTTGAGCCGCACCACGCAACGGCCCTGTCGCGCGTGGCGCAGCATCAGGCTCTGAATGAACGCTTGCGGTGTATAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCgg                                                                                                                                             <  1:945769/150‑1 (MQ=255)
         ccgccCTTGAGCCGCACCACGCAACGGCCCTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTg                                                                                                                                         >  2:363689/1‑150 (MQ=255)
              cTTGAGCCGCACCACGCAACGGCCCTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAg                                                                                                                                    >  2:869654/1‑150 (MQ=255)
              cTTGAGCCGCACCACGCAACGGCCCTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAg                                                                                                                                    <  1:753941/150‑1 (MQ=255)
               ttGAGCCGCACCACGCAACGGCCCTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAgg                                                                                                                                   >  1:2477/1‑150 (MQ=255)
                    ccGCACCACGCAACGGCCCTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGgtcgtc                                                                                                                              <  2:523326/150‑1 (MQ=255)
                        accacGCAACGGCCCTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTAcc                                                                                                                          <  2:1309880/150‑1 (MQ=255)
                          cacGCAACGGCCCTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCAt                                                                                                                        >  2:384365/1‑150 (MQ=255)
                                     ccTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCAcggcgg                                                                                                             >  2:914728/1‑150 (MQ=255)
                                     ccTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCAcggcgg                                                                                                             <  1:303536/150‑1 (MQ=255)
                                      cTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGAAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACggcggc                                                                                                            >  1:447328/1‑150 (MQ=255)
                                          cgcgcgTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTtgc                                                                                                        <  1:202638/150‑1 (MQ=255)
                                                    gcagcaTCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAacgcac                                                                                              >  1:1364649/1‑150 (MQ=255)
                                                            aGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAGTGCTCGAGCATTGccc                                                                                                                                              >  2:948050/1‑93 (MQ=255)
                                                            aGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAGTGCTCGAGCATTGccc                                                                                                                                              <  1:948050/94‑2 (MQ=255)
                                                                gcTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGtt                                                                                                                                          >  2:501380/1‑94 (MQ=255)
                                                                gcTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGtt                                                                                                                                          <  1:501380/94‑1 (MQ=255)
                                                                gcTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGa                                                                                  >  1:511634/1‑150 (MQ=255)
                                                                  tgaatgaaCGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAgg                                                                                <  2:511634/150‑1 (MQ=255)
                                                                      tgaaCGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGa                                                                            <  1:947586/150‑1 (MQ=255)
                                                                      tgaaCGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGa                                                                            <  2:2477/150‑1 (MQ=255)
                                                                      ggaaCGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGa                                                                            <  2:195643/149‑1 (MQ=255)
                                                                                cGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAgg                                                                  <  1:363689/150‑1 (MQ=255)
                                                                                           ggCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCAcccggtt                                                       <  2:815140/150‑6 (MQ=255)
                                                                                                       gcTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCACCAGGGCCGGCACCGACCa                                           <  2:447328/150‑1 (MQ=255)
                                                                                                       gcTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAAGGTGAGCAGTTCAGGGTCACCAGGGCCGGCCCCGACCa                                           >  2:184387/1‑150 (MQ=255)
                                                                                                         ttGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTGCCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCACCAGGGCCGGCACCGACCAGc                                         <  2:674680/150‑1 (MQ=255)
                                                                                                         ttGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACAGCGGCTGGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCACCAGGGCCGGCACCGACCAGc                                         >  2:549601/1‑150 (MQ=255)
                                                                                                             ccACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGa                                                                            >  1:734420/1‑111 (MQ=255)
                                                                                                             ccACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGa                                                                            <  2:734420/111‑1 (MQ=255)
                                                                                                                  gTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACGGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGAAGTTCAGGGTCACCAGGGCCGGCACCGACCAGCCATACtttt                                <  2:591940/150‑1 (MQ=255)
                                                                                                                      ttATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCACCAGGGCCGGCACCGACCAGCCATACTTTTGCAt                            >  1:716725/1‑150 (MQ=255)
                                                                                                                         tCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCACCAGGGCCGGCACCGACCAGCCATACTTTTGCATtca                         >  2:1136611/1‑150 (MQ=255)
                                                                                                                           cGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCACCAGGGCCGGCACCGACCAGCCATACTTTTGCATtcatc                       <  1:384365/150‑1 (MQ=255)
                                                                                                                                 tGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCACCAGGGCCGGCACCGACCAGCCATACTTTTGCATTCATCAGCGtt                 <  1:869654/150‑1 (MQ=255)
                                                                                                                                 acccggCAGTGCTCGAGCATTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCACCAGGGCCGGCACCGACCAGCCATACtttt                                <  1:711651/131‑1 (MQ=255)
                                                                                                                                 acccggCAGTGCTCGAGCATTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCACCAGGGCCGGCACCGACCAGCCATACtttt                                >  2:711651/5‑135 (MQ=255)
                                                                                                                                   cGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCACCAGGGCCGGCACCGACCAGCCATACTTTTGCATTCATCAg                     <  2:1314024/144‑1 (MQ=255)
                                                                                                                                   cGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCACCAGGGCCGGCACCGACCAGCCATACTTTTGCATTCATCAg                     >  1:1314024/1‑144 (MQ=255)
                                                                                                                                               gAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAgg                                                                  <  2:242904/87‑1 (MQ=255)
                                                                                                                                               gAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAgg                                                                  >  1:242904/1‑87 (MQ=255)
                                                                                                                                                aGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCCCGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCACCAGGGCCGGCACCGACCAGCCATACTTTTGCATTCATCAGCGTTTCCTCATCGGCAAgg  >  1:249383/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CAAGGGTCGCCGCCCTTGAGCCGCACCACGCAACGGCCCTGTCGCGCGTGGCGCAGCATCAGGCGCTGAATGAACGCTTGCGGTGTAGAGCGGCAGCCGCCGCGCTTGCCCACGGTCATTATCCGGGCCTGCGGGCAATGCTCGAGCACTGCCGGGTTGACGAGGTCGTCTACCATGACCACGGCGGCTTGCTGCAACGCACGCACTGCCTTGAGGGTGAGCAGTTCAGGGTCACCAGGGCCGGCACCGACCAGCCATACTTTTGCATTCATCAGCGTTTCCTCATCGGCAAGG  >  NC_002947/2381970‑2382263

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