Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 584,788 T→A 66.7% I369I (ATA→ATT cydC ← fused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE584,7880TA66.7% 51.2 / 36.8 54I369I (ATA→ATTcydCfused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base T (15/3);  new base A (28/8);  total (43/11)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

GGTCCCATGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTC  >  minE/584721‑584847
                                                                   |                                                           
ggTCCCATGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCta                                                            <  1:109590/69‑1 (MQ=255)
ggTCCCATGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCta                                                            <  1:1874045/69‑1 (MQ=255)
ggTCCCATGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCta                                                            <  1:2369885/69‑1 (MQ=255)
ggTCCCATGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCta                                                            <  1:2389412/69‑1 (MQ=255)
ggTCCCATGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCta                                                            <  1:1545235/69‑1 (MQ=255)
ggTCCCATGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCta                                                            <  1:2659394/69‑1 (MQ=255)
ggTCCCATGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCta                                                            <  1:2750141/69‑1 (MQ=255)
ggTCCCATGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCta                                                            <  1:1251072/69‑1 (MQ=255)
ggTCCCATGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCta                                                            <  1:341830/69‑1 (MQ=255)
 gTCCCATGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTAGCCCGAGAATCGCta                                                            <  1:3394160/68‑1 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGCATCGCTATATGTTcccc                                                  >  1:1633486/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTcccc                                                  >  1:3415589/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTcccc                                                  >  1:2737789/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTcccc                                                  >  1:530842/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTcccc                                                  >  1:243907/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTcccc                                                  >  1:2423796/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTcccc                                                  >  1:650173/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTcccc                                                  >  1:88530/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTcccc                                                  >  1:1972102/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTcccc                                                  >  1:2380198/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTcccc                                                  >  1:1781266/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTcccc                                                  >  1:1766109/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTcccc                                                  >  1:1615290/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTcccc                                                  >  1:1632693/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTcccc                                                  >  1:1636598/1‑69 (MQ=255)
          gcgGGTCAGCAGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCAATATGTTccc                                                   >  1:51493/1‑68 (MQ=255)
                 aGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCta                                                            <  1:616924/52‑1 (MQ=255)
                             aGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGaa                               <  1:2023650/69‑1 (MQ=255)
                             aGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGaa                               <  1:1100124/69‑1 (MQ=255)
                             aGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGaa                               <  1:1286643/69‑1 (MQ=255)
                                   gATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATccc                          >  1:949302/1‑68 (MQ=255)
                                       tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAagag                                 >  1:966876/1‑57 (MQ=255)
                                       tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga                                >  1:962986/1‑58 (MQ=255)
                                       tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga                                >  1:1941431/1‑58 (MQ=255)
                                       tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga                                >  1:1164450/1‑58 (MQ=255)
                                       tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga                                >  1:1235738/1‑58 (MQ=255)
                                       tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga                                >  1:2131580/1‑58 (MQ=255)
                                       tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga                                >  1:3161298/1‑58 (MQ=255)
                                       tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga                                >  1:1489610/1‑58 (MQ=255)
                                       tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAgaga                                >  1:2865872/1‑58 (MQ=255)
                                                cGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTg              >  1:1682363/1‑67 (MQ=255)
                                                cGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTg              >  1:866304/1‑67 (MQ=255)
                                                cGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTg              >  1:2075250/1‑67 (MQ=255)
                                                cGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTg              >  1:2914413/1‑67 (MQ=255)
                                                 ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt                       >  1:2093628/1‑57 (MQ=255)
                                                 ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt                       >  1:1845131/1‑57 (MQ=255)
                                                 ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt                       >  1:1046937/1‑57 (MQ=255)
                                                 ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt                       >  1:2151806/1‑57 (MQ=255)
                                                 ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt                       >  1:688586/1‑57 (MQ=255)
                                                 ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt                       >  1:113975/1‑57 (MQ=255)
                                                 ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt                       >  1:61813/1‑57 (MQ=255)
                                                 ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt                       >  1:2515739/1‑57 (MQ=255)
                                                 ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt                       >  1:1226913/1‑57 (MQ=255)
                                                 ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt                       >  1:2588093/1‑57 (MQ=255)
                                                 ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt                       >  1:2932022/1‑57 (MQ=255)
                                                 ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt                       >  1:2743267/1‑57 (MQ=255)
                                                 ggTTCGCCCGAGAATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCttt                       >  1:2899566/1‑57 (MQ=255)
                                                          gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc  <  1:482215/69‑1 (MQ=255)
                                                          gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc  <  1:1299601/69‑1 (MQ=255)
                                                          gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc  <  1:2464840/69‑1 (MQ=255)
                                                          gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc  <  1:1539882/69‑1 (MQ=255)
                                                          gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc  <  1:324056/69‑1 (MQ=255)
                                                          gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc  <  1:2716144/69‑1 (MQ=255)
                                                          gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc  <  1:2036439/69‑1 (MQ=255)
                                                          gagaATCGCAATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGGGATTGCTc  <  1:220931/69‑1 (MQ=255)
                                                                   |                                                           
GGTCCCATGCGCGGGTCAGCTGTTGTAACAGTGTTGATTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTC  >  minE/584721‑584847

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: