Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,577,649 A→G 30.2% I96I (ATT→ATC xapA ← purine nucleoside phosphorylase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,577,6490AG30.2% 59.1 / 29.9 66I96I (ATT→ATCxapApurine nucleoside phosphorylase II
Reads supporting (aligned to +/- strand):  ref base A (9/37);  new base G (6/14);  total (15/51)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GCATTGGTGCAGAACAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCC  >  minE/1577587‑1577712
                                                              |                                                               
gCATTGGTGCAGAACAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGt                                                            <  1:323788/68‑1 (MQ=255)
gCATTGGTGCAGAACAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGt                                                            <  1:1499705/68‑1 (MQ=255)
gCATTGGTGCAGAACAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGt                                                            <  1:2211548/68‑1 (MQ=255)
gCATTGGTGCAGAACAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGt                                                            <  1:1690958/68‑1 (MQ=255)
gCATTGGTGCAGAACAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGt                                                            <  1:1763562/68‑1 (MQ=255)
gCATTGGTGCAGAACAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGt                                                            <  1:3151951/68‑1 (MQ=255)
             aCAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCtt                                              <  1:2850299/69‑1 (MQ=255)
             aCAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCtt                                              <  1:2660104/69‑1 (MQ=255)
             aCAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCtt                                              <  1:2649441/69‑1 (MQ=255)
             aCAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCtt                                              <  1:2864327/69‑1 (MQ=255)
             aCAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCtt                                              <  1:2500457/69‑1 (MQ=255)
             aCAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCtt                                              <  1:1537908/69‑1 (MQ=255)
             aCAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCtt                                              <  1:2942397/69‑1 (MQ=255)
                 tAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGt                                                            <  1:408162/51‑1 (MQ=255)
                   tcTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGa                                        <  1:2058758/68‑1 (MQ=255)
                   gCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGa                                        <  1:500370/69‑1 (MQ=255)
                   acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                        <  1:3136085/68‑1 (MQ=255)
                   acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                        <  1:3129046/68‑1 (MQ=255)
                   acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                        <  1:1968899/68‑1 (MQ=255)
                   acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                        <  1:2756534/68‑1 (MQ=255)
                   acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                        <  1:470521/68‑1 (MQ=255)
                   acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                        <  1:1672377/68‑1 (MQ=255)
                   acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                        <  1:533281/68‑1 (MQ=255)
                   acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                        <  1:1483195/68‑1 (MQ=255)
                   acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                        <  1:1434304/68‑1 (MQ=255)
                   acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                        <  1:1286382/68‑1 (MQ=255)
                   acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                        <  1:956046/68‑1 (MQ=255)
                   acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                        <  1:1114913/68‑1 (MQ=255)
                   acTCGCAGCCCACCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                        <  1:3119312/68‑1 (MQ=255)
                    cTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                        <  1:346590/68‑1 (MQ=255)
                     tCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGa                                        <  1:801404/67‑1 (MQ=255)
                               gcagcTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGAc                            <  1:985368/69‑1 (MQ=255)
                               gcagcTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGAc                            <  1:182667/69‑1 (MQ=255)
                               gcagcTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGAc                            <  1:1201954/69‑1 (MQ=255)
                               gcagcTTAAAGGTACGAATTGCGTCAGTCATAATGGGCATTCCACGACCTTCGTAGAAATGTCCGCGAc                            <  1:411253/69‑1 (MQ=255)
                                   cTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCttt                        <  1:2366548/69‑1 (MQ=255)
                                   cTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCttt                        <  1:1523376/69‑1 (MQ=255)
                                        aGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTcata                    <  1:981849/68‑1 (MQ=255)
                                        aGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTcata                    <  1:195012/68‑1 (MQ=255)
                                        aGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTcata                    <  1:3253432/68‑1 (MQ=255)
                                         ggTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTcata                    <  1:742293/67‑1 (MQ=255)
                                         ggTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTcata                    <  1:2095732/67‑1 (MQ=255)
                                           tACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCttt                        >  1:3008061/1‑61 (MQ=255)
                                           tACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCttt                        >  1:2439525/1‑61 (MQ=255)
                                           tACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCttt                        >  1:221260/1‑61 (MQ=255)
                                             cGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCttt                        <  1:899445/59‑1 (MQ=255)
                                                   gCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGt        >  1:642171/1‑69 (MQ=255)
                                                      tcagtcaTAATGGTCATTCCACGACCTTCGTAGAAATGt                                   >  1:967670/1‑39 (MQ=255)
                                                      tcagtcaTAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGCCCTTTCATACATACCACCGGTACg     >  1:1587490/1‑69 (MQ=255)
                                                      tcagtcaTAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACg     >  1:826590/1‑69 (MQ=255)
                                                      tcagtcaTAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACg     >  1:2551322/1‑69 (MQ=255)
                                                      tcagtcaTAATGGTCATTCCACGACCTTCGTAGAAATGTCGCGACCTTTCATACATACCACCGGTACGc    <  1:275711/69‑1 (MQ=255)
                                                       cagtcaTAATGTTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGc    <  1:2460676/69‑1 (MQ=255)
                                                       cagtcaTAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGc    <  1:787502/69‑1 (MQ=255)
                                                       cagtcaTAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGc    <  1:1001092/69‑1 (MQ=255)
                                                       cagtcaTAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGc    <  1:3264091/69‑1 (MQ=255)
                                                       cagtcaTAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGc    <  1:2312264/69‑1 (MQ=255)
                                                        cgtcaTGATGGTCATCCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACACGCCACTGGCAcccc   >  1:365099/2‑69 (MQ=255)
                                                        cgtcaTGATGGTCATCCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACACGCCACTGGCAcccc   >  1:3379435/2‑69 (MQ=255)
                                                        cgtcaTGATGGTCATCCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACACGCCACTGGCAcccc   >  1:962384/2‑69 (MQ=255)
                                                        cgtcaTGATGGTCATCCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACACGCCACTGGCAcccc   >  1:3031274/2‑69 (MQ=255)
                                                        cgtcaTGATGGTCATCCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACACGCCACTGGCAcccc   >  1:2760857/2‑69 (MQ=255)
                                                        cgtcaTGATGGTCATCCCAAGACCTTCGTAGAAATGTCCGCGACCTTTCATACACGCCACTGGCAcccc   >  1:273404/2‑69 (MQ=255)
                                                        agtcaTAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGc    <  1:1638485/68‑1 (MQ=255)
                                                        agtcaTAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGc    <  1:2120887/68‑1 (MQ=255)
                                                            aTAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGccc  >  1:635384/1‑66 (MQ=255)
                                                              |                                                               
GCATTGGTGCAGAACAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCC  >  minE/1577587‑1577712

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: