Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 9,238 G→A 34.7% intergenic (+47/‑68) talB → / → mog transaldolase B/predicted molybdochelatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE9,2380GA34.7% 50.0 / 36.4 69intergenic (+47/‑68)talB/mogtransaldolase B/predicted molybdochelatase
Reads supporting (aligned to +/- strand):  ref base G (42/3);  new base A (9/15);  total (51/18)
Fisher's exact test for biased strand distribution p-value = 1.12e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAA  >  minE/9179‑9292
                                                           |                                                      
cGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctgtctgt                                               >  1:3208337/1‑69 (MQ=255)
cGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctgtctgt                                               >  1:2644501/1‑69 (MQ=255)
cGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctgtctgt                                               >  1:2530757/1‑69 (MQ=255)
cGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctgtctgt                                               >  1:733581/1‑69 (MQ=255)
cGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctgtctgt                                               >  1:2014885/1‑69 (MQ=255)
cGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctgtctgt                                               >  1:1664740/1‑69 (MQ=255)
cGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctgtctg                                                <  1:1230009/68‑1 (MQ=255)
cGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctgtctg                                                <  1:2295663/68‑1 (MQ=255)
cGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctgtctg                                                <  1:2252385/68‑1 (MQ=255)
    ctgctgTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctg                                                    >  1:2083894/1‑60 (MQ=255)
    ctgctgTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctg                                                    >  1:378912/1‑60 (MQ=255)
    ctgctgTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctg                                                    >  1:380926/1‑60 (MQ=255)
    ctgctgTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctg                                                    >  1:2440151/1‑60 (MQ=255)
    ctgctgTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctg                                                    >  1:601761/1‑60 (MQ=255)
    ctgctgTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctg                                                    >  1:1613023/1‑60 (MQ=255)
    ctgctgTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctg                                                    >  1:1315523/1‑60 (MQ=255)
    ctgctgTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctg                                                    >  1:1411014/1‑60 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCTGTCTGTCACTCCCTTCGCAgtg                               <  1:1011727/67‑1 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCTGTCTGTCACTCCCTTCGCAgtg                               <  1:1619589/67‑1 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCTGTCTGTCACTCCCTTCGCAgtg                               <  1:3055933/67‑1 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCTGTCTGTCACTCCCTTCGCAgtg                               <  1:3358692/67‑1 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCTGTCTGTCACTCCCTTCGCAgtg                               <  1:3401741/67‑1 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCTGTCTGTCACTCCCTTCGCAgtg                               <  1:1098695/67‑1 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCTGTCTGTCACTCCCTTCGCAgtg                               <  1:284696/67‑1 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCTGTCTGTCACTCCCTTCGCAgtg                               <  1:2694342/67‑1 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCTGTCTGTCACTCCCTTCGCAgtg                               <  1:1261834/67‑1 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCTGTCTGTCACTCCCTTCGCAgtg                               <  1:1006743/67‑1 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCTGTCTGTCACTCCCTTCGCAgtg                               <  1:1363775/67‑1 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCTGTCTGTCACTCCCTTCGCAgtg                               <  1:64763/67‑1 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCTGTCTGTCACTCCCTTCGCAgtg                               <  1:1635318/67‑1 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCTGTCTGTCACTCCCTTCGCAgtg                               <  1:1527855/67‑1 (MQ=255)
                 acTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAACCCGTCTGTCACTCCCTTCGCAgtg                               <  1:794430/67‑1 (MQ=255)
                     aGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATc                           >  1:2630063/1‑68 (MQ=255)
                            ccGGGAAGTCGGTCACGCTACCTCTTCTGAAACctg                                                    >  1:23116/1‑36 (MQ=255)
                            ccGGGAAGTCGGTCACGCTACCTCTTCTGAAACctg                                                    >  1:17265/1‑36 (MQ=255)
                            ccGGGAAGTCGGTCACGCTACCTCTTCTGAAACctg                                                    >  1:1832899/1‑36 (MQ=255)
                            ccGGGAAGTCGGTCACGCTACCTCTTCTGAAACctg                                                    >  1:1878389/1‑36 (MQ=255)
                            ccGGGAAGTCGGTCACGCTACCTCTTCTGAAACctg                                                    >  1:802261/1‑36 (MQ=255)
                            ccGGGAAGTCGGTCACGCTACCTCTTCTGAAACctg                                                    >  1:2145272/1‑36 (MQ=255)
                            ccGGGAAGTCGGTCACGCTACCTCTTCTGAAACctg                                                    >  1:1539680/1‑36 (MQ=255)
                            ccGGGAAGTCGGTCACGCTACCTCTTCTGAAACctg                                                    >  1:2443145/1‑36 (MQ=255)
                            ccGGGAAGTCGGTCACGCTACCTCTTCTGAAACctg                                                    >  1:869972/1‑36 (MQ=255)
                             cGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:582593/1‑69 (MQ=255)
                             cGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:590218/1‑69 (MQ=255)
                             cGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:883078/1‑69 (MQ=255)
                             cGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:748709/1‑69 (MQ=255)
                             cGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:2578730/1‑69 (MQ=255)
                             cGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:3243501/1‑69 (MQ=255)
                             cGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:1794488/1‑69 (MQ=255)
                             cGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:2693209/1‑69 (MQ=255)
                             cGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:2299487/1‑69 (MQ=255)
                             cGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGtt                   >  1:962286/1‑68 (MQ=255)
                                      ggTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGc                                   >  1:2944328/1‑43 (MQ=255)
                                      ggTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGc                                   >  1:1177203/1‑43 (MQ=255)
                                      ggTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGc                                   >  1:287187/1‑43 (MQ=255)
                                      ggTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGc                                   >  1:1970798/1‑43 (MQ=255)
                                      ggTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGc                                   >  1:2267825/1‑43 (MQ=255)
                                      ggTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGc                                   >  1:2325871/1‑43 (MQ=255)
                                          aCGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:1559130/1‑56 (MQ=255)
                                          aCGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:482686/1‑56 (MQ=255)
                                          aCGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:3389957/1‑56 (MQ=255)
                                             cTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTaaa  >  1:1480639/1‑69 (MQ=255)
                                             cTACCTCTTCTGAAGCCTGTCTGTAACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTaaa  >  1:2691091/1‑69 (MQ=255)
                                                   cttctGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:992273/1‑47 (MQ=255)
                                                   cttctGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:843478/1‑47 (MQ=255)
                                                   cttctGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:2709896/1‑47 (MQ=255)
                                                   cttctGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:2904896/1‑47 (MQ=255)
                                                   cttctGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGttt                  >  1:3276877/1‑47 (MQ=255)
                                                   cttctGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGtt                   >  1:346390/1‑46 (MQ=255)
                                                           |                                                      
CGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAA  >  minE/9179‑9292

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: