Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,373,313 G→T 11.4% R105S (CGC→AGC)  ubiD ← 3‑octaprenyl‑4‑hydroxybenzoate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,373,3130GT11.4% 183.4 / 6.3 61R105S (CGC→AGC) ubiD3‑octaprenyl‑4‑hydroxybenzoate decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base T (3/4);  total (30/31)
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.12e-01

ACCACGCAGCCGCTTTGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACTTCACGCAGCGCCGAAA  >  minE/2373245‑2373381
                                                                    |                                                                    
accacGCAGCCGCTTTGTCGGCATGTTCAATAATTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGa                                                                    <  1:2228231/71‑1 (MQ=255)
 ccacGCAGCCGCTTTGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGa                                                                    <  1:2796405/70‑1 (MQ=255)
 ccacGCAGCCGCTTTGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGa                                                                    >  1:1537850/1‑70 (MQ=255)
  cacGCAGCCGCTTTGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCgg                                                                     <  1:4214923/68‑1 (MQ=255)
   acGCAGCCGCTTTGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAAc                                                                 >  1:2194837/1‑71 (MQ=255)
   acGCAGCCGCTTTGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGATTATCAAACAGGTCGCGGAAAc                                                                 >  1:3332525/1‑71 (MQ=255)
    cGCAGCCGCTTTGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAAcc                                                                >  1:1136171/1‑71 (MQ=255)
    cGCAGCCGCTTTGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAAcc                                                                >  1:2406037/1‑71 (MQ=255)
    cGCAGCCGCTTTGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAAcc                                                                <  1:3490006/71‑1 (MQ=255)
       aGCCGCTTTGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCttt                                                             <  1:2053695/71‑1 (MQ=255)
       aGCCGCTTTGTCGGCATGTTCAATAATTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCttt                                                             <  1:3439768/71‑1 (MQ=255)
        gCCGCTTTGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCtttt                                                            <  1:3989228/71‑1 (MQ=255)
         ccGCTTTGTCGGCATGTTCAATACTTGCTTAAACTGGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTgg                                                          >  1:1759359/1‑71 (MQ=255)
          cGCTTTGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTgg                                                          >  1:3833974/1‑71 (MQ=255)
               tGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTggcggc                                                      <  1:2555951/70‑1 (MQ=255)
               tGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTggcggc                                                      <  1:2550582/70‑1 (MQ=255)
                 tCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCt                                                               <  1:970109/59‑1 (MQ=255)
                      aTGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTGGcggctccggctc                                              >  1:3300556/1‑71 (MQ=255)
                           cAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCttt                                           >  1:3803461/1‑69 (MQ=255)
                           cAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCttt                                           >  1:1475655/1‑69 (MQ=255)
                            aaTACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAg                                        <  1:3868877/71‑1 (MQ=255)
                              tACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGaa                                                                   <  1:2798217/42‑1 (MQ=255)
                              tACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTGGtggctccggc                                                <  1:3415860/61‑1 (MQ=255)
                               ccttgcccggcgttaatgatttgccCAAACAGGTCGCTGaaa                                                                  >  1:1435596‑M2/26‑42 (MQ=255)
                               ccttgcccggcgttaatgatttgccCAAACAGGTCGCTGa                                                                    <  1:2202077‑M2/15‑1 (MQ=255)
                               aCTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGaaa                                     >  1:378812/1‑71 (MQ=255)
                                cTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTgg                                                          >  1:1911739/1‑49 (MQ=255)
                                cTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTg                                                           <  1:459160/48‑1 (MQ=255)
                                cTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTg                                                           <  1:1853750/48‑1 (MQ=255)
                                cTTGCTTAAACTGCGGCAGTTTATCAAACAGGGCGCGGAAACCt                                                               <  1:2714533/44‑1 (MQ=255)
                                 ttgcccggcgttaatgatttgccCAAACAGGTCGCTGaaa                                                                  >  1:282622‑M2/24‑40 (MQ=255)
                                 ttgcccggcgttaatgatttgccCAAACAGGTCGCTGaaa                                                                  <  1:4694086‑M2/17‑1 (MQ=255)
                                 ttgcccggcgttaatgatttgccCAAACAGGTCGCTGa                                                                    <  1:183108‑M2/15‑1 (MQ=255)
                                  tgcccggcgttaatgatttgccCAAACAGGTCGCTGaa                                                                   >  1:4628403‑M2/23‑38 (MQ=255)
                                  tgcccggcgttaatgatttgccCAAACAGGTCGCTGa                                                                    <  1:1307213‑M2/15‑1 (MQ=255)
                                   gCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGaaa                                     >  1:2970999/1‑67 (MQ=255)
                                   gCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGc                                  >  1:225623/1‑70 (MQ=255)
                                       aaaCTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGaaa                                     >  1:1948559/1‑63 (MQ=255)
                                       aaaCTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCaat                              >  1:1329353/1‑70 (MQ=255)
                                         aCTGCGGCAGTTTATCAAACAGGTCGCGTAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCaataat                           <  1:4589126/71‑1 (MQ=255)
                                          cTGCGGCAGTTTATCAAACAGGTCGCGGAAACCtttt                                                            <  1:3510144/37‑1 (MQ=255)
                                            gCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGaaa                                     >  1:2974262/1‑58 (MQ=255)
                                              ggCAGTTTATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTAcc                      <  1:216303/71‑1 (MQ=255)
                                                   tttATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGctct                                             <  1:1930010/43‑1 (MQ=255)
                                                   tttATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACtt                 <  1:1838317/71‑1 (MQ=255)
                                                    ttATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTAcc                      <  1:2915724/65‑1 (MQ=255)
                                                     tataAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACTTCa               >  1:1875606/1‑71 (MQ=255)
                                                     tATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTc                                          >  1:886822/1‑44 (MQ=255)
                                                     tATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTc                                          >  1:4341359/1‑44 (MQ=255)
                                                     tATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCaataa                            >  1:245774/1‑58 (MQ=255)
                                                     tATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACTTc                <  1:1348547/70‑1 (MQ=255)
                                                     tATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACTTCa               >  1:3982837/1‑71 (MQ=255)
                                                       tCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGaaa                                     >  1:4454827/1‑47 (MQ=255)
                                                       tCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACTTCACg             <  1:2403261/71‑1 (MQ=255)
                                                       tCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACTTCACg             >  1:70173/1‑71 (MQ=255)
                                                        cAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACTTCACg             <  1:3504913/70‑1 (MQ=255)
                                                         aaaCAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCa                     <  1:1639995/61‑1 (MQ=255)
                                                            cAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACTTCCCGCAgcg        >  1:889348/1‑71 (MQ=255)
                                                            cAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACTTCACGCAgcg        >  1:2762416/1‑71 (MQ=255)
                                                              ggTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGa                                       <  1:4450637/38‑1 (MQ=255)
                                                               gTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACt                  <  1:2302146/58‑1 (MQ=255)
                                                               gTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACTTCACGCAGCGcc      >  1:93870/1‑70 (MQ=255)
                                                                 cgcgGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCaat                              >  1:1688794/1‑44 (MQ=255)
                                                                  gcgGAAACCTTTTGGCGTCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACTTCACGCAGCGCCGaaa  >  1:2931629/1‑71 (MQ=255)
                                                                  gcgGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACTTCACGCAGCGCCGaa   <  1:2280624/70‑1 (MQ=255)
                                                                    |                                                                    
ACCACGCAGCCGCTTTGTCGGCATGTTCAATACTTGCTTAAACTGCGGCAGTTTATCAAACAGGTCGCGGAAACCTTTTGGCGGCTCCGGCTCTTTCAGAAACGCCAATAATTTACCAACTTCACGCAGCGCCGAAA  >  minE/2373245‑2373381

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: