Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,779,576 C→T 55.7% S302S (TCG→TCA ffh ← Signal Recognition Particle (SRP) component with 4.5S RNA (ffs)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,779,5760CT55.7% 15.4 / 62.8 79S302S (TCG→TCAffhSignal Recognition Particle (SRP) component with 4.5S RNA (ffs)
Reads supporting (aligned to +/- strand):  ref base C (10/25);  new base T (27/17);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 6.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

AGCTTGCTGGCTAATTTCTCTGCCTGCGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGC  >  minE/1779514‑1779633
                                                              |                                                         
aGCTTGCTGGCTAATTTCTCTGCCTGCGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAg                                                      <  1:1084503/68‑1 (MQ=255)
  cTTGCTGGCTAATTTCTCTGCCTGCGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTAc                                                   >  1:1771532/1‑69 (MQ=255)
        ggCTAATTTCTCTGCCTGCGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGcc                                             <  1:718467/69‑1 (MQ=255)
        ggCTAATTTCTCTGCCTGCACGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGcc                                             <  1:568842/69‑1 (MQ=255)
        ggCTAATTTCTCTGCCTGCACGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGcc                                             <  1:2921366/69‑1 (MQ=255)
        ggCTAATTTCTCTGCCTGCACGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGcc                                             <  1:1186813/69‑1 (MQ=255)
        ggCTAATTTCTCTGCCTGCACGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGcc                                             <  1:3319283/69‑1 (MQ=255)
                tctcTGCCTGCGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGCCCATGCCg                                      >  1:1119392/1‑68 (MQ=255)
                tctcTGCCTGCGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGCCCATGCCg                                      >  1:606736/1‑68 (MQ=255)
                tctcTGCCTGCGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGCCCATGCCg                                      >  1:1088139/1‑68 (MQ=255)
                tctcTGCCTGCGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGCCCATGCCg                                      >  1:2879908/1‑68 (MQ=255)
                tctcTGCCTGCGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGCCCATGCCg                                      >  1:1052680/1‑68 (MQ=255)
                tctcTGCCTGCGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGCCCATGCCg                                      >  1:1838912/1‑68 (MQ=255)
                tctcTGCCTGCGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGCCCATGCCg                                      >  1:1065896/1‑68 (MQ=255)
                      cctgcGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGaa                                  <  1:2620374/66‑1 (MQ=255)
                      cctgcGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGaa                                  <  1:948028/66‑1 (MQ=255)
                      cctgcGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGaa                                  <  1:2100838/66‑1 (MQ=255)
                        tgcGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATacgc                             <  1:399430/69‑1 (MQ=255)
                        tgcGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATacgc                             <  1:3260916/69‑1 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:212494/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:275705/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:2816514/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:2670090/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:2902263/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:2298664/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:2960011/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:3063775/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:2234551/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:838424/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:2157221/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:140780/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:1898010/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:1860394/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:1525873/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:875172/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:116071/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggga                           >  1:1106342/1‑68 (MQ=255)
                           gcgcgGTCAACCTTGCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATacggg                            >  1:755325/1‑65 (MQ=255)
                                         gCTTTCGATATCTTCGATCAGTGACAGTACGTCGc                                              >  1:176525/1‑35 (MQ=255)
                                         gCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATACGGGAGGCAATACGGTCCgg            >  1:697971/1‑69 (MQ=255)
                                         gCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATACGGGAGGCAATACGGTCCgg            >  1:735357/1‑69 (MQ=255)
                                         gCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATACGGGAGGCAATACGGTCCgg            >  1:2922752/1‑69 (MQ=255)
                                         gCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATACGGGAGGCAATACGGTCCgg            >  1:1259791/1‑69 (MQ=255)
                                         gCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATACGGGAGGCAATACGGTCCgg            >  1:1303950/1‑69 (MQ=255)
                                         gCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATACGGGAGGCAATACGGTCCgg            >  1:2486217/1‑69 (MQ=255)
                                         gCTTTCGATATCTTCGATCAGTGACAGTACGTCGCCCATGCCGAGAATACGGGAGGCAATACGGTCCg             >  1:1570805/1‑68 (MQ=255)
                                               gATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCTATGCGGTCCGGATgg        <  1:1072393/67‑1 (MQ=255)
                                               gATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATgg        <  1:1179813/67‑1 (MQ=255)
                                               gATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATgg        <  1:1954055/67‑1 (MQ=255)
                                               gATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATgg        <  1:205862/67‑1 (MQ=255)
                                               gATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATgg        <  1:3074017/67‑1 (MQ=255)
                                               gATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATgg        <  1:31589/67‑1 (MQ=255)
                                               gATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATgg        <  1:3161034/67‑1 (MQ=255)
                                               gATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATgg        <  1:3205580/67‑1 (MQ=255)
                                               gATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGa       <  1:809821/68‑1 (MQ=255)
                                               gATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGa       <  1:3398208/68‑1 (MQ=255)
                                               gATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGa       <  1:2291146/68‑1 (MQ=255)
                                               gATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGa       <  1:1290014/68‑1 (MQ=255)
                                               gATATCTTCGATCAGCGACAGTACGGCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATgg        <  1:464721/67‑1 (MQ=255)
                                                atatCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATgg        <  1:173907/66‑1 (MQ=255)
                                                     ttCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGc  >  1:2568135/1‑67 (MQ=255)
                                                     ttCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGc  >  1:2829810/1‑67 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:1571855/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:963945/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:1089552/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:820480/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:774846/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:666778/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:624387/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:463847/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:328889/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:2828252/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:2678477/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:2258710/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:2198576/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:1715908/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:1638458/37‑1 (MQ=255)
                                                      tCGATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:1597262/37‑1 (MQ=255)
                                                        gATCAGTGACAGTACGTCGCCCATGCCGAGAATac                               <  1:3236278/35‑1 (MQ=255)
                                                              |                                                         
AGCTTGCTGGCTAATTTCTCTGCCTGCGCGCGGTCAACTTTGCTTTCGATATCTTCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGC  >  minE/1779514‑1779633

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: