Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,774,792 G→A 28.0% R333R (CGG→CGA aceK → isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,774,7920GA28.0% 75.8 / 27.9 71R333R (CGG→CGAaceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base G (18/33);  new base A (8/12);  total (26/45)
Fisher's exact test for biased strand distribution p-value = 7.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

AGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACG  >  minE/2774725‑2774850
                                                                   |                                                          
aGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCggg                                                           <  1:758315/69‑1 (MQ=255)
aGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTCATCggg                                                           <  1:435687/69‑1 (MQ=255)
 gTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCggg                                                           <  1:1399780/68‑1 (MQ=255)
   tCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGGTCGGGTa                                                         >  1:2757653/1‑68 (MQ=255)
   tCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTa                                                         >  1:843140/1‑68 (MQ=255)
    cATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTAt                                                        >  1:805465/1‑68 (MQ=255)
                   tATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCaaa                                          >  1:1077199/1‑67 (MQ=255)
                     ttCGTGGAATTGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCa                                                     <  1:2741893/54‑1 (MQ=255)
                     ttCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCa                                                     <  1:2783398/54‑1 (MQ=255)
                     ttCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCa                                                     <  1:264874/54‑1 (MQ=255)
                     ttCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCa                                                     <  1:2487611/54‑1 (MQ=255)
                       cGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAgg                                    >  1:1508392/1‑69 (MQ=255)
                       cGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAgg                                    >  1:715175/1‑69 (MQ=255)
                                   aTGTTGGTGTTTACGTTGCCGGGTTTTGATCGAGTATTCaaa                                                   >  1:3115139/1‑42 (MQ=37)
                                   aTGTTGGTGTTTACGTTGCCGGGTTTTGATCGAGTATTCaaa                                                   >  1:2608842/1‑42 (MQ=37)
                                   aTGTTGGTGTTTACGTTGCCGGGTTTTGATCGAGTATTCaaa                                                   >  1:2573050/1‑42 (MQ=37)
                                   aTGTTGGTGTTTACGTTGCCGGGTTTTGATCGAGTATTCaaa                                                   >  1:2184792/1‑42 (MQ=37)
                                   aTGTTGGTGTTTACGTTGCCGGGTTTTGATCGAGTATTCaaa                                                   >  1:1374660/1‑42 (MQ=37)
                                   aTGTTGGTGTTTACGTTGCCGGGTTTTGATCGAGTATTCaaa                                                   >  1:623513/1‑42 (MQ=37)
                                   aTGTTGGTGTTTACGTTGCC‑‑GTTTTGATCGAGTATTCaaa                                                   >  1:3428388/1‑40 (MQ=16)
                                      ttGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACa                                      >  1:2148069/1‑52 (MQ=255)
                                      ttGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGa                       >  1:583641/1‑67 (MQ=255)
                                      ttGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGa                       >  1:1144603/1‑67 (MQ=255)
                                        tctgtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAaga                   <  1:2302573/67‑1 (MQ=255)
                                        ggTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGtc                                                >  1:611450/1‑40 (MQ=255)
                                        ggTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGa                       >  1:3074498/1‑65 (MQ=255)
                                        ggTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAaga                   <  1:3165856/69‑1 (MQ=255)
                                        ggTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAaga                   <  1:301741/69‑1 (MQ=255)
                                        ggTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAaga                   <  1:2486489/69‑1 (MQ=255)
                                        ggTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAaga                   <  1:47656/69‑1 (MQ=255)
                                        ggTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAaga                   <  1:3139916/69‑1 (MQ=255)
                                        ggTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAaga                   <  1:1569928/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:2980429/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:3301801/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:3072006/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:1356062/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:2782321/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:2621477/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:2742467/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:848359/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:818209/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:547320/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:184914/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:548447/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:2318951/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:640481/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:567538/69‑1 (MQ=255)
                                           gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:2106254/69‑1 (MQ=255)
                                            tttACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGAt                <  1:2798618/68‑1 (MQ=255)
                                            tttACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATg               <  1:319115/69‑1 (MQ=255)
                                            tttACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATg               <  1:2115352/69‑1 (MQ=255)
                                            tttACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATg               <  1:1279428/69‑1 (MQ=255)
                                              tACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAaga                   <  1:2519885/63‑1 (MQ=255)
                                              tACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTc             >  1:1919972/1‑69 (MQ=255)
                                              tACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTc             >  1:65936/1‑69 (MQ=255)
                                              tACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTc             >  1:773514/1‑69 (MQ=255)
                                              tACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTc             >  1:2074650/1‑69 (MQ=255)
                                              tACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTc             >  1:1666487/1‑69 (MQ=255)
                                              tACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTc             >  1:333965/1‑69 (MQ=255)
                                              tACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTc             >  1:2814027/1‑69 (MQ=255)
                                                  cTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAaga                   <  1:203717/59‑1 (MQ=255)
                                                         gttttGATCGAGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACg  >  1:1587342/3‑69 (MQ=255)
                                                            ttGATCGAGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGa                       <  1:1047356/45‑1 (MQ=255)
                                                            ttGATCGAGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGa                       <  1:1323385/45‑1 (MQ=255)
                                                            ttGATCGAGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGa                       <  1:1176225/45‑1 (MQ=255)
                                                            ttGATCGAGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGa                       <  1:307293/45‑1 (MQ=255)
                                                            ttGATCGAGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGa                       <  1:9698/45‑1 (MQ=255)
                                                            ttGATCGAGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGa                       <  1:195750/45‑1 (MQ=255)
                                                            ttGATCGAGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGa                       <  1:247943/45‑1 (MQ=255)
                                                            ttGATCGAGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGa                       <  1:2602565/45‑1 (MQ=255)
                                                            ttGATCGAGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGa                       <  1:1178309/45‑1 (MQ=255)
                                                            ttGATCGAGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGa                       <  1:3181644/45‑1 (MQ=255)
                                                            ttGATCGAGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGa                       <  1:3037076/45‑1 (MQ=255)
                                                            ttGATCGAGGATTCAAAGTCATCAAAGACAGGTGCGCGCCGCAGa                       <  1:1952370/45‑1 (MQ=255)
                                                                   |                                                          
AGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACG  >  minE/2774725‑2774850

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: