Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 722,240 G→C 54.4% A461A (GCC→GCG rne ← fused ribonucleaseE endoribonuclease and scaffold for formation of degradosome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE722,2400GC54.4% 14.5 / 89.2 79A461A (GCC→GCGrnefused ribonucleaseE endoribonuclease and scaffold for formation of degradosome
Reads supporting (aligned to +/- strand):  ref base G (10/26);  new base C (10/33);  total (20/59)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GTGCGGGGTTTCCATCTGATCGTTTGGCACAATTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCGATTGGCACAGGAACA  >  minE/722174‑722297
                                                                  |                                                         
gTGCGGGGTTTCCATCTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:972465/69‑1 (MQ=255)
gTGCGGGGTTTCCATCTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:1894990/69‑1 (MQ=255)
gTGCGGGGTTTCCATCTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:1774638/69‑1 (MQ=255)
gTGCGGGGTTTCCATCTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:2405669/69‑1 (MQ=255)
gTGCGGGGTTTCCATCTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:1726251/69‑1 (MQ=255)
gTGCGGGGTTTCCATCTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:2409776/69‑1 (MQ=255)
gTGCGGGGTTTCCATCTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:2561557/69‑1 (MQ=255)
gTGCGGGGTTTCCATCTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:2821992/69‑1 (MQ=255)
gTGCGGGGTTTCCATCTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:1329330/69‑1 (MQ=255)
gTGCGGGGTTTCCATCTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:1305203/69‑1 (MQ=255)
gTGCGGGGTTTCCATCTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:45486/69‑1 (MQ=255)
gTGCGGGGTTTCCATCTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:508168/69‑1 (MQ=255)
gTGCGGGGTTTCCATCTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:511943/69‑1 (MQ=255)
gTGCGGGGTTTCCATCTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:900354/69‑1 (MQ=255)
       gTTTCCATCTGATCGTTTGGCACAATTACACAGCGCACACCGTCCTGACGAGTTTCAATGTCATTTAc                                                   >  1:1011095/1‑68 (MQ=255)
       gTTTCCATCTGATCGTTTGGCACAATTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTAc                                                   >  1:1347115/1‑68 (MQ=255)
         ttCCATCTGATCGTTTGGCACAATTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCAttt                                                     >  1:1846899/1‑64 (MQ=255)
         ttCCATCTGATCGTTTGGCACAATTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTAc                                                   >  1:1781288/1‑66 (MQ=255)
         ttCCATCTGATCGTTTGGCACAATTAAAAATCGCACACCGTCCTGACGAGTTTCAATGGCATTTAc                                                   >  1:3383236/1‑66 (MQ=255)
            cATCTGATCGTTTGGCACAATTACAGAGCGCTCACCGTCCTGACGAGTTTCAATGGCATTTACCGCAg                                              <  1:2279040/68‑1 (MQ=255)
            cATCTGATCGTTTGGCACAATTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAg                                              <  1:2282984/68‑1 (MQ=255)
               cTGATCGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGc                                                         <  1:2514841/54‑1 (MQ=255)
                    cGTTTGTCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:1170325/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:2954475/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:2620300/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:2581874/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:639392/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:547251/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:1634571/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:2719732/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:1475876/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:237748/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:311350/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:3172087/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:221766/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:2165726/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:2161832/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:1221047/65‑1 (MQ=255)
                    cGTTTGGCACAATCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:440671/65‑1 (MQ=255)
                              aaTTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTcagcag                           <  1:389123/69‑1 (MQ=255)
                              aaTTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTcagcag                           <  1:3144497/69‑1 (MQ=255)
                              aaTTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTcagcag                           <  1:1079069/69‑1 (MQ=255)
                              aaTTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTcagcag                           <  1:2755908/69‑1 (MQ=255)
                              aaTTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTcagcag                           <  1:3404927/69‑1 (MQ=255)
                              aaTTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTcagcag                           <  1:667697/69‑1 (MQ=255)
                              aaTTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTcagcag                           <  1:740117/69‑1 (MQ=255)
                              aaTTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTcagcag                           <  1:90120/69‑1 (MQ=255)
                              aaTTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTcagcag                           <  1:96829/69‑1 (MQ=255)
                              aaTTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTcagcag                           <  1:1740668/69‑1 (MQ=255)
                              aaTTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTcagcag                           <  1:18053/69‑1 (MQ=255)
                              aaTTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTcagcag                           <  1:2249626/69‑1 (MQ=255)
                              aaTTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTcagcag                           <  1:2234293/69‑1 (MQ=255)
                              aaTCACACAGCGTACGCCGTCCTGACGAGTTTCAATCGCATTTACCGCAGAACGt                                         <  1:1075336/55‑1 (MQ=255)
                                       gcgcACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGc                  <  1:3281205/69‑1 (MQ=255)
                                       gcgcACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGc                  <  1:2080034/69‑1 (MQ=255)
                                              ccGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCGATTgg            >  1:2702380/1‑68 (MQ=255)
                                              ccGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCGATTgg            >  1:3423573/1‑68 (MQ=255)
                                              ccGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCGATTgg            >  1:2046239/1‑68 (MQ=255)
                                              ccGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCGATTgg            >  1:1837388/1‑68 (MQ=255)
                                              ccGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCGATTgg            >  1:2388916/1‑68 (MQ=255)
                                                   cTGACGAGTTTCAATCGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCAATTGGcaca        >  1:1040218/1‑67 (MQ=255)
                                                   cTGACGAGTTTCAATCGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCAATTGGcaca        >  1:2951886/1‑67 (MQ=255)
                                                   cTGACGAGTTTCAATCGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCAATTGGcaca        >  1:1933456/1‑67 (MQ=255)
                                                   cTGACGAGTTTCAATCGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCAATTGGcaca        >  1:3348646/1‑67 (MQ=255)
                                                   cTGACGAGTTTCAATCGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCAATTGGcaca        >  1:3326801/1‑67 (MQ=255)
                                                   cTGACGAGTTTCAATCGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCAATTGGcaca        >  1:1733691/1‑67 (MQ=255)
                                                   cTGACGAGTTTCAATCGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCAATTGGcaca        >  1:1659288/1‑67 (MQ=255)
                                                   cTGACGAGTTTCAATCGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCAATTGGcaca        >  1:1516256/1‑67 (MQ=255)
                                                   cTGACGAGTTTCAATCGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCAATTGGcaca        >  1:2555809/1‑67 (MQ=255)
                                                   cTGACGAGTTTCAATCGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCAATTGGcac         >  1:2905721/1‑66 (MQ=255)
                                                        gAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCGATTGGCACAGGAACa  <  1:2328721/68‑1 (MQ=255)
                                                                aTGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCg                 <  1:2625737/45‑1 (MQ=255)
                                                                aTGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCg                 <  1:9493/45‑1 (MQ=255)
                                                                aTGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCg                 <  1:25147/45‑1 (MQ=255)
                                                                aTGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCg                 <  1:432990/45‑1 (MQ=255)
                                                                aTGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCg                 <  1:539871/45‑1 (MQ=255)
                                                                aTGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCg                 <  1:536036/45‑1 (MQ=255)
                                                                aTGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCg                 <  1:522085/45‑1 (MQ=255)
                                                                aTGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCg                 <  1:1311136/45‑1 (MQ=255)
                                                                  |                                                         
GTGCGGGGTTTCCATCTGATCGTTTGGCACAATTACACAGCGCACACCGTCCTGACGAGTTTCAATGGCATTTACCGCAGAACGTTTTTCATTCAGCAGGTAAGAAGCGATTGGCACAGGAACA  >  minE/722174‑722297

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: