Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,631,750 G→T 26.4% A14A (GCG→GCT livM → leucine/isoleucine/valine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,631,7500GT26.4% 46.8 / 13.5 34A14A (GCG→GCTlivMleucine/isoleucine/valine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (11/14);  new base T (7/2);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

CTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTGGATGGCACCAAACTGGTG  >  minE/2631683‑2631816
                                                                   |                                                                  
cTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGa                                                                   >  1:2599055/1‑69 (MQ=255)
cTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGa                                                                   >  1:2676293/1‑69 (MQ=255)
cTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGa                                                                   >  1:2294475/1‑69 (MQ=255)
cTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGa                                                                   >  1:2693919/1‑69 (MQ=255)
cTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGa                                                                   >  1:3006672/1‑69 (MQ=255)
cTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGa                                                                   >  1:3257243/1‑69 (MQ=255)
cTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGa                                                                   >  1:1647182/1‑69 (MQ=255)
            gAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCTTTGTGc                                                       <  1:889389/69‑1 (MQ=255)
             aGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCTTTGTGc                                                       <  1:3078406/68‑1 (MQ=255)
                  gagaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCt                                                             >  1:634277/1‑57 (MQ=255)
                  gagaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCt                                                             >  1:938691/1‑57 (MQ=255)
                  gagaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCt                                                             >  1:2083148/1‑57 (MQ=255)
                  gagaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCt                                                             >  1:1335351/1‑57 (MQ=255)
                  gagaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCt                                                             >  1:1221124/1‑57 (MQ=255)
                  gagaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCt                                                             >  1:2345044/1‑57 (MQ=255)
                  gagaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCt                                                             >  1:2500035/1‑57 (MQ=255)
                    gaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTggcgggc                                               <  1:1161865/69‑1 (MQ=255)
                         tATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCtt                                          <  1:270570/69‑1 (MQ=255)
                         tATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCtt                                          <  1:2440938/69‑1 (MQ=255)
                              aaCCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTggcgggc                                               <  1:1191421/59‑1 (MQ=255)
                                  gATGCATATTGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTggcgggc                                               <  1:1184634/55‑1 (MQ=255)
                                          ttGGAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCtt                                          <  1:888944/52‑1 (MQ=255)
                                          ttGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCtt                                          <  1:2095395/52‑1 (MQ=255)
                                          ttGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCtt                                          <  1:808427/52‑1 (MQ=255)
                                               aTGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTg                    <  1:1917672/69‑1 (MQ=255)
                                                  gcgcTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTg                    <  1:2413408/66‑1 (MQ=255)
                                                  gcgcTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTg                    <  1:2938006/66‑1 (MQ=255)
                                                  gcgcTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTg                    <  1:1784076/66‑1 (MQ=255)
                                                  gcgcTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTg                    <  1:1746936/66‑1 (MQ=255)
                                                  gcgcTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTg                    <  1:1382218/66‑1 (MQ=255)
                                                  gcgcTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTGGAt                 >  1:81463/1‑69 (MQ=255)
                                                          ctctGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTGGATGGCACCaa         >  1:724202/1‑69 (MQ=255)
                                                          ctctGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTGGATGGCACCaa         >  1:1373647/1‑69 (MQ=255)
                                                          ctctGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTGGATGGCACCa          >  1:1538133/1‑68 (MQ=255)
                                                                 gcgATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTGGATGGCACCAAACtggtg  <  1:1957312/69‑1 (MQ=255)
                                                                   |                                                                  
CTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTGGATGGCACCAAACTGGTG  >  minE/2631683‑2631816

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: