Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 720,216 T→C 32.0% intergenic (+14/+221) mviN → / ← rne predicted inner membrane protein/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
*minE720,2160TC32.0% 24.0 / 17.5 25intergenic (+14/+221)mviN/rnepredicted inner membrane protein/fused ribonucleaseE endoribonuclease and scaffold for formation of degradosome
Reads supporting (aligned to +/- strand):  ref base T (9/8);  new base C (5/3);  total (14/11)
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.63e-01

CACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGTCTGCAGTGCAGGCCGGAGATAATCTTTTCGCTTTAAAACATATCATGAAACTGGGTATGTTTTGTC  >  minE/720150‑720282
                                                                  |                                                                  
cACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGTCtgc                                                                <  1:607485/71‑1 (MQ=255)
cACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtg                                                                 <  1:1573393/70‑1 (MQ=255)
     gCGGTACTTGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtg                                                            <  1:1294940/70‑1 (MQ=255)
     gCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtgc                                                           <  1:3230532/71‑1 (MQ=255)
          aCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGTCTGCAGTGCAGGcc                                                      >  1:2435195/1‑71 (MQ=255)
              ggCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGTCtgcagtgc                                                           <  1:175774/62‑1 (MQ=255)
                   cAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATaa                                              >  1:482230/1‑70 (MQ=255)
                   cAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATaa                                              >  1:2142042/1‑70 (MQ=255)
                      aGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGTCTGCAGTGCAGGCCGGAGAt                                                <  1:3749155/65‑1 (MQ=255)
                      aGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCgg                                                    >  1:3344642/1‑61 (MQ=255)
                      aGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCt                                           >  1:721213/1‑70 (MQ=255)
                        ttAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGTCTGCAGTGCAGGCCGGAGATaa                                              >  1:2835238/1‑65 (MQ=255)
                        ttAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGTCTGCAGTGCAGGCCGGAGATAATCtttt                                        <  1:911173/71‑1 (MQ=255)
                        ttAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtgcag                                                         >  1:631701/1‑54 (MQ=255)
                          aaaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGTCtgca                                                               >  1:2959041/1‑46 (MQ=255)
                               aTTTGCCCGCCGGACGGTGTAACAATGCATTCCGGTCTGCAGTGCAGGCCGGAGATAATCTTTTCGCttt                                  >  1:708286/1‑70 (MQ=255)
                                  tGCCCGCCGGACGGTGTAACAATGCATTCCGGTCTGCAGTGCAGGCCGGAGATAATCTTTTCGCTTTaaaa                              >  1:1125406/1‑71 (MQ=255)
                                      cgccgGACGGTGTAACAATGCATTCCGGTCTGCAGTGCAGGCCGGAGATAATCTTTTCGCTTTAAAACata                          <  1:373740/71‑1 (MQ=255)
                                       gccgGACGGTGTAACAATGCATTCCGGTCTGCAGTGCAGGCCGgaga                                                 >  1:1614756/1‑47 (MQ=255)
                                        ccgGACGGTGTAACAATGCATTCCGGTCTGCAGTGCAGGCCGGAGATAATCtt                                          >  1:1104064/1‑53 (MQ=255)
                                         cgGACGGTGTAACAATGCATTCCGGTCTGCAGTGCAGGCCGGAGATAATCTTTTCGCTTTAAAACATATCa                       >  1:1653247/1‑71 (MQ=255)
                                                tgtAACAATGCATTCCGGTCTGCAGTGCAGGCCGGAGATAATCTTTTCGCTTTAAAACATATCATGaaa                  <  1:422232/69‑1 (MQ=255)
                                                        tGCATTCCGGTCTGCAGTGCAGGCCGGAGATAATCTTTTCGCTTTAAAACATATCATGAAACTGGGTATGt        <  1:108330/71‑1 (MQ=255)
                                                          cATTCCGGTCTGCAGTGCAGGCCGGAGATAATCTTTTCGCTTTAAAACATATCATGAAACt                <  1:433860/61‑1 (MQ=255)
                                                              ccGGTCTGCAGTGCAGGCCGGAGATAATCTTTTCGCTTTAAAACATATCATGAAACTGGGTATGTTTTGTc  >  1:1945392/1‑71 (MQ=255)
                                                                  tCTGCAGTGCAGGCCGGAGATAATCTTTTCGCTTTAAAACATATCATGaaa                  >  1:425866/1‑51 (MQ=255)
                                                                  |                                                                  
CACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGTCTGCAGTGCAGGCCGGAGATAATCTTTTCGCTTTAAAACATATCATGAAACTGGGTATGTTTTGTC  >  minE/720150‑720282

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: