Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,879,631 T→C 37.6% E179E (GAA→GAG mazG ← nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,879,6310TC37.6% 60.8 / 70.6 109E179E (GAA→GAGmazGnucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (43/25);  new base C (20/21);  total (63/46)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

AACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCG  >  minE/1879572‑1879697
                                                           |                                                                  
aagagcagGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATc                                                           >  1:351312/4‑69 (MQ=255)
aaCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCTTGCCGCGCTTCGTACATc                                                           >  1:1982335/1‑69 (MQ=255)
aaCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATc                                                           >  1:2257717/1‑69 (MQ=255)
aaCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATc                                                           >  1:223439/1‑69 (MQ=255)
aaCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATc                                                           >  1:2224375/1‑69 (MQ=255)
aaCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATc                                                           >  1:2219855/1‑69 (MQ=255)
aaCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATc                                                           <  1:1587960/69‑1 (MQ=255)
aaCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATc                                                           >  1:265743/1‑69 (MQ=255)
aaCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATc                                                           >  1:1931530/1‑69 (MQ=255)
aaCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATc                                                           >  1:3023798/1‑69 (MQ=255)
aaCAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATc                                                           >  1:3019809/1‑69 (MQ=255)
       ggTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCtcg                                                     >  1:1534617/1‑68 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:1466496/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:156772/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:3110190/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:2924432/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:1928440/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:2006282/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:2039251/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:2526795/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:2502357/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:2296974/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:1007911/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:3327470/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:1019144/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:3383776/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:828266/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:1203352/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:385261/1‑61 (MQ=255)
               aTTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:810567/1‑61 (MQ=255)
               aTTTCATCCTCCAGTTTAGCCTGGTCGACAACAGACTGCCGCGCCTCGTACATCACCtcgt                                                    >  1:660592/1‑61 (MQ=255)
                  tcctcctccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTc                                                                  >  1:1053850/1‑44 (MQ=255)
                  tcctcctccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTAc                                                              >  1:2644166/1‑48 (MQ=255)
                  tcctcctccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAt                                                >  1:811491/1‑62 (MQ=255)
                  tcctcctccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAt                                                >  1:1318063/1‑62 (MQ=255)
                  tcctcctccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAt                                                >  1:1928949/1‑62 (MQ=255)
                  tcctcctccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAt                                                >  1:1371387/1‑62 (MQ=255)
                  tcctcctccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGAt                                                >  1:2657113/1‑62 (MQ=255)
                      cctccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTACGTACATCACCTCGTCGATCTCTTCGTAGa                                     >  1:1255438/1‑69 (MQ=255)
                       ctccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGAc                                    >  1:1261088/1‑69 (MQ=255)
                         ccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGAc                                    <  1:2597432/67‑1 (MQ=255)
                         ccAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGAc                                    <  1:2182377/67‑1 (MQ=255)
                            gTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACttt                                 <  1:1302232/67‑1 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCtcgtcg                                                  >  1:1428773/1‑47 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCtcgtcg                                                  >  1:1297293/1‑47 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCtcgtcg                                                  >  1:1106758/1‑47 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCtcgtcg                                                  >  1:2946784/1‑47 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCtcgtc                                                   >  1:206575/1‑46 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                              >  1:2899427/1‑67 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                              >  1:3008250/1‑67 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                              <  1:361821/67‑1 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                              >  1:1013505/1‑67 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                              <  1:3286069/67‑1 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                              <  1:3289947/67‑1 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                              >  1:213582/1‑67 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                              >  1:1182135/1‑67 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                              <  1:1976044/67‑1 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                              >  1:1422707/1‑67 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                              <  1:1411230/67‑1 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                              >  1:1860309/1‑67 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                              >  1:1774443/1‑67 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCg                             >  1:389998/1‑68 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCg                             >  1:1510629/1‑68 (MQ=255)
                               tAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCg                             >  1:73486/1‑68 (MQ=255)
                                  ccTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACt                          <  1:3162627/68‑1 (MQ=255)
                                  ccTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACt                          <  1:3052668/68‑1 (MQ=255)
                                  ccTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACt                          <  1:1632967/68‑1 (MQ=255)
                                  ccTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACt                          <  1:2770436/68‑1 (MQ=255)
                                  ccTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACt                          <  1:78396/68‑1 (MQ=255)
                                  ccTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACt                          <  1:2222302/68‑1 (MQ=255)
                                  ccTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACt                          <  1:1098019/68‑1 (MQ=255)
                                   cTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATc                              <  1:1316920/63‑1 (MQ=255)
                                         gACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGa                    >  1:1650753/1‑67 (MQ=255)
                                         gACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGa                    >  1:1100994/1‑67 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:824218/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:51184/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:457312/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:649243/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:438930/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:3207826/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:750027/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:762496/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:1067494/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:1026603/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:854568/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:3203311/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:1560619/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:865698/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:1944750/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:1951763/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:1996213/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:2314004/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:2280306/41‑1 (MQ=255)
                                                     ccGCGCCTCGTACATCACCTCGTCGATCTCTTCGTAGACtt                                  <  1:2333034/41‑1 (MQ=255)
                                                         gcTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCg  <  1:918191/69‑1 (MQ=255)
                                                         gcTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCg  <  1:605721/69‑1 (MQ=255)
                                                         gcTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCg  <  1:1536637/69‑1 (MQ=255)
                                                         gcTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCg  <  1:1721504/69‑1 (MQ=255)
                                                         gcTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCg  <  1:2604443/69‑1 (MQ=255)
                                                          cTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCg  <  1:1427935/68‑1 (MQ=255)
                                                           |                                                                  
AACAGCAGGTCCCCCATTTCCTCCTCCAGTTTAGCCTGGTCGACAACAGCCTGCCGCGCTTCGTACATCACCTCGTCGATCTCTTCGTAGACTTTATCGACTACCGGACCAAGCGTCGTCCAATCG  >  minE/1879572‑1879697

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: