Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,577,658 T→C 32.7% G93G (GGA→GGG xapA ← purine nucleoside phosphorylase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,577,6580TC32.7% 45.6 / 30.4 61G93G (GGA→GGGxapApurine nucleoside phosphorylase II
Reads supporting (aligned to +/- strand):  ref base T (10/31);  new base C (6/14);  total (16/45)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGATGG  >  minE/1577600‑1577721
                                                          |                                                               
aCAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCtt                                                       <  1:2942397/69‑1 (MQ=255)
aCAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCtt                                                       <  1:2500457/69‑1 (MQ=255)
aCAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCtt                                                       <  1:2864327/69‑1 (MQ=255)
aCAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCtt                                                       <  1:2649441/69‑1 (MQ=255)
aCAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCtt                                                       <  1:2850299/69‑1 (MQ=255)
aCAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCtt                                                       <  1:1537908/69‑1 (MQ=255)
aCAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCtt                                                       <  1:2660104/69‑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:533281/68‑1 (MQ=255)
      acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                                 <  1:470521/68‑1 (MQ=255)
      acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                                 <  1:1968899/68‑1 (MQ=255)
      acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                                 <  1:1672377/68‑1 (MQ=255)
      acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                                 <  1:2756534/68‑1 (MQ=255)
      acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                                 <  1:1483195/68‑1 (MQ=255)
      acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                                 <  1:1434304/68‑1 (MQ=255)
      acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                                 <  1:956046/68‑1 (MQ=255)
      acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                                 <  1:1286382/68‑1 (MQ=255)
      acTCGCAGCCCAGCAACTTAAAGGTACGGATTGCATCCGTCATGATGGTCATCCCACGACCTTCGTAGa                                                 <  1:3129046/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:1201954/69‑1 (MQ=255)
                  gcagcTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGAc                                     <  1:182667/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:195012/68‑1 (MQ=255)
                           aGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTcata                             <  1:3253432/68‑1 (MQ=255)
                           aGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTcata                             <  1:981849/68‑1 (MQ=255)
                            ggTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTcata                             <  1:2095732/67‑1 (MQ=255)
                            ggTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTcata                             <  1:742293/67‑1 (MQ=255)
                              tACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCttt                                 >  1:221260/1‑61 (MQ=255)
                              tACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCttt                                 >  1:2439525/1‑61 (MQ=255)
                              tACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCttt                                 >  1:3008061/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:2551322/1‑69 (MQ=255)
                                         tcagtcaTAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACg              >  1:826590/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:3031274/2‑69 (MQ=255)
                                           cgtcaTGATGGTCATCCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACACGCCACTGGCAcccc            >  1:2760857/2‑69 (MQ=255)
                                           cgtcaTGATGGTCATCCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACACGCCACTGGCAcccc            >  1:3379435/2‑69 (MQ=255)
                                           cgtcaTGATGGTCATCCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACACGCCACTGGCAcccc            >  1:365099/2‑69 (MQ=255)
                                           cgtcaTGATGGTCATCCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACACGCCACTGGCAcccc            >  1:962384/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)
                                                      tCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGATgg  >  1:2593684/1‑68 (MQ=255)
                                                        aTTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGc             <  1:2293681/55‑1 (MQ=255)
                                                          tCCACGACCTTCGTAGAAATGTCCGCGACCTTTca                               <  1:1363106/35‑1 (MQ=255)
                                                          tCCACGACCTTCGTAGAAATGTCCGCGACCTTTca                               <  1:1282952/35‑1 (MQ=255)
                                                          tCCACGACCTTCGTAGAAATGTCCGCGACCTTTca                               <  1:1207240/35‑1 (MQ=255)
                                                          tCCACGACCTTCGTAGAAATGTCCGCGACCTTTca                               <  1:91263/35‑1 (MQ=255)
                                                          tCCACGACCTTCGTAGAAATGTCCGCGACCTTTca                               <  1:1052151/35‑1 (MQ=255)
                                                          tCCACGACCTTCGTAGAAATGTCCGCGACCTTTca                               <  1:882873/35‑1 (MQ=255)
                                                          tCCACGACCTTCGTAGAAATGTCCGCGACCTTTca                               <  1:2086762/35‑1 (MQ=255)
                                                          tCCACGACCTTCGTAGAAATGTCCGCGACCTTTca                               <  1:2104535/35‑1 (MQ=255)
                                                          tCCACGACCTTCGTAGAAATGTCCGCGACCTTTca                               <  1:645169/35‑1 (MQ=255)
                                                          tCCACGACCTTCGTAGAAATGTCCGCGACCTTTca                               <  1:616405/35‑1 (MQ=255)
                                                          tCCACGACCTTCGTAGAAATGTCCGCGACCTTTca                               <  1:2664011/35‑1 (MQ=255)
                                                          tCCACGACCTTCGTAGAAATGTCCGCGACCTTTca                               <  1:2935748/35‑1 (MQ=255)
                                                          tCCACGACCTTCGTAGAAATGTCCGCGACCTTTca                               <  1:3018278/35‑1 (MQ=255)
                                                          |                                                               
ACAGTAGCTCGCAGCCCAGCAGCTTAAAGGTACGAATTGCGTCAGTCATAATGGTCATTCCACGACCTTCGTAGAAATGTCCGCGACCTTTCATACATACCACCGGTACGCCCTGGAGATGG  >  minE/1577600‑1577721

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: