Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,421,968 C→T 5.3% P163P (CCG→CCA)  ccmF ← heme lyase, CcmF subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,421,9680CT 5.3% 236.1 / 2.3 75P163P (CCG→CCA) ccmFheme lyase, CcmF subunit
Reads supporting (aligned to +/- strand):  ref base C (31/40);  new base T (2/2);  total (33/42)
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.16e-01

CACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAG  >  minE/1421899‑1422036
                                                                     |                                                                    
cACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCg                                                                     <  1:2617290/71‑1 (MQ=255)
    gagaAACCCACGTACCCCATATAAAGCAGAGGCCGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGt                                                                  >  1:3840769/1‑70 (MQ=255)
      gaAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCgc                                                                    >  1:1174800/1‑65 (MQ=255)
      gaAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTaa                                                               <  1:1556441/71‑1 (MQ=255)
       aaaCCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAg                                                              <  1:1375730/71‑1 (MQ=255)
         aCCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGAt                                                            <  1:3444788/71‑1 (MQ=255)
         aCCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGAt                                                            <  1:1490369/71‑1 (MQ=255)
                tACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAACAGTGGGTTa                                                                <  1:2874820/60‑1 (MQ=255)
                  ccccATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCtt                                                   <  1:3597909/71‑1 (MQ=255)
                    ccATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATcg                                                          >  1:2250262/1‑62 (MQ=255)
                    ccATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGAcc                                                     <  1:1696886/67‑1 (MQ=255)
                    ccATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGAc                                                      <  1:1155149/66‑1 (MQ=255)
                     cATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCaa                                                <  1:521641/71‑1 (MQ=255)
                      atatAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAg                                                              <  1:4200204/56‑1 (MQ=255)
                      atatAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGa                                                             <  1:2271675/57‑1 (MQ=255)
                       tataAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGtt                                                                 <  1:312058/52‑1 (MQ=255)
                       tataAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATc                                              <  1:4175596/71‑1 (MQ=255)
                        ataAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCg                                             >  1:1527851/1‑71 (MQ=255)
                         taAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATc                                              >  1:2197769/1‑69 (MQ=255)
                         taAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCgg                                            >  1:1950210/1‑71 (MQ=255)
                         taAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCg                                             <  1:4470507/70‑1 (MQ=255)
                         taAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCATGCAATAGCGGGTTAAGATCGCGACCTTCAATCgg                                            >  1:1504416/1‑71 (MQ=255)
                           aaGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAACAGTGGGTTAAGATCACGCCCTTCAATCGGGa                                          <  1:110618/71‑1 (MQ=255)
                             gCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGTTTAAGATcgc                                                         <  1:720494/54‑1 (MQ=255)
                             gCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGa                                                             >  1:3630223/1‑50 (MQ=255)
                                   gCGGATGGAAGATCAGCCCCGGATCCTGCAATATCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGc                                  <  1:2323918/71‑1 (MQ=255)
                                   gCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGc                                  <  1:129131/71‑1 (MQ=255)
                                     ggATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCaa                                >  1:1095699/1‑71 (MQ=255)
                                     ggATGGAAGATCAGCCCCGGATCCTGCAACAGTGGGTTAAGATCACGCCCTTCAATCGGGAAATTCGGCag                                >  1:3501241/1‑70 (MQ=255)
                                     ggATGGAAGATCAGCCCCGGATCCTGCAACAGTGGGTTAAGATAACGCCCTTCAATCGGGAAATTCGGCag                                >  1:3777843/1‑70 (MQ=255)
                                      gATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCggg                                           <  1:2710327/59‑1 (MQ=255)
                                       aTGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACg                              <  1:1846898/71‑1 (MQ=255)
                                        tGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGt                             >  1:807456/1‑71 (MQ=255)
                                          gAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGAc                                                      >  1:2163388/1‑44 (MQ=255)
                                          gAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGt                             >  1:1983340/1‑69 (MQ=255)
                                           aaGAGCAGCCCCGGAGCCTGCAATAGCGGGTTAAGATCGCTACCTTCAATCGGGAAGGTCGGCAACGTGCg                          <  1:2979957/71‑1 (MQ=255)
                                            aGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCg                                    >  1:1907142/1‑60 (MQ=255)
                                            aGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCga                         <  1:3651497/71‑1 (MQ=255)
                                             gATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATcg                                                          <  1:2656467/37‑1 (MQ=255)
                                             gATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGc                           <  1:1060759/68‑1 (MQ=255)
                                             gATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGc                           <  1:3364582/68‑1 (MQ=255)
                                             gATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgag                        <  1:3842619/71‑1 (MQ=255)
                                              aTCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTg                            <  1:1763299/66‑1 (MQ=255)
                                              aTCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCga                         >  1:3583813/1‑69 (MQ=255)
                                               tCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgagag                      <  1:73783/71‑1 (MQ=255)
                                               tCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgagag                      <  1:2700669/71‑1 (MQ=255)
                                               tCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgagag                      <  1:3602473/71‑1 (MQ=255)
                                                cAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgaga                       >  1:2388965/1‑69 (MQ=255)
                                                cAGCCCCGGATCATGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGa                                          >  1:739160/1‑50 (MQ=255)
                                                   ccccGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACg                  <  1:2688688/71‑1 (MQ=255)
                                                        gATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGa                                          <  1:2625496/42‑1 (MQ=255)
                                                        gATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgag                        <  1:2926917/60‑1 (MQ=255)
                                                        gATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACggg                >  1:22067/1‑68 (MQ=255)
                                                        gATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGt               >  1:3149117/1‑69 (MQ=255)
                                                        gATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGt               >  1:483191/1‑69 (MQ=255)
                                                         aTCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTag            <  1:933947/71‑1 (MQ=255)
                                                         aTCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTag            <  1:426034/71‑1 (MQ=255)
                                                          tCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTaga           <  1:3127737/71‑1 (MQ=255)
                                                             tGCAATAGCGGGTTAAGATCGCGACCTTCAATCggg                                           <  1:2760493/36‑1 (MQ=255)
                                                              gCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGt        >  1:4574115/1‑70 (MQ=255)
                                                                 aTAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAag    <  1:2348358/71‑1 (MQ=255)
                                                                  tAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGaa                    <  1:294102/54‑1 (MQ=255)
                                                                  tAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTaaa     >  1:2606226/1‑69 (MQ=255)
                                                                   aGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgaga                       >  1:61486/1‑50 (MQ=255)
                                                                   aGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGaa                    >  1:2045871/1‑53 (MQ=255)
                                                                    gCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCg                                    >  1:1626247/1‑36 (MQ=255)
                                                                    gCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGGAACGTGCGAGAGAACGGGTTAGAGGTAAagag  >  1:489536/1‑70 (MQ=255)
                                                                    gCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCaa                                <  1:262897/40‑1 (MQ=255)
                                                                    gCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAGCGTGCGAGAGAACGGGTTAGAGGTAAagag  >  1:3279799/1‑70 (MQ=255)
                                                                    gCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACg                              >  1:4042784/1‑42 (MQ=255)
                                                                    gCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgaga                       >  1:356683/1‑49 (MQ=255)
                                                                    gCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCga                         >  1:2658792/1‑47 (MQ=255)
                                                                    gCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGaa                    >  1:1503591/1‑52 (MQ=255)
                                                                    gCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAagag  <  1:4581011/70‑1 (MQ=255)
                                                                    gCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAagag  >  1:4463460/1‑70 (MQ=255)
                                                                     cGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCgagag                      >  1:2594916/1‑49 (MQ=255)
                                                                     |                                                                    
CACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCGGGAAGTTCGGCAACGTGCGAGAGAACGGGTTAGAGGTAAAGAG  >  minE/1421899‑1422036

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