Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 131,549 G→A 39.1% intergenic (‑115/‑159) panD ← / → yadD aspartate 1‑decarboxylase/predicted transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE131,5490GA39.1% 48.4 / 65.6 97intergenic (‑115/‑159)panD/yadDaspartate 1‑decarboxylase/predicted transposase
Reads supporting (aligned to +/- strand):  ref base G (43/16);  new base A (10/28);  total (53/44)
Fisher's exact test for biased strand distribution p-value = 8.85e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

CTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTG  >  minE/131494‑131614
                                                       |                                                                 
ctcACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCg                                                      >  1:2397208/1‑69 (MQ=255)
ctcACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCg                                                      >  1:783610/1‑69 (MQ=255)
          aTTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCgcgc                                              >  1:1037344/1‑67 (MQ=255)
          aTTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCgcgc                                              >  1:370325/1‑67 (MQ=255)
          aTTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCgcgc                                              >  1:2898884/1‑67 (MQ=255)
           ttGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCgcgcag                                            >  1:780126/1‑66 (MQ=255)
           ttGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGa                                           <  1:778377/69‑1 (MQ=255)
           ttGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGa                                           >  1:3035697/1‑69 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTcagc                                                              >  1:1074793/1‑48 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTcag                                                               >  1:1994895/1‑47 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAg                                         >  1:2759544/1‑69 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAg                                         >  1:2736903/1‑69 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAg                                         >  1:3359716/1‑69 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAg                                         >  1:33618/1‑69 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAg                                         >  1:3783/1‑69 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAg                                         >  1:45107/1‑69 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAg                                         >  1:178199/1‑69 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAg                                         >  1:1777627/1‑69 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAg                                         >  1:1775587/1‑69 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAg                                         >  1:1158453/1‑69 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAg                                         >  1:1323462/1‑69 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAg                                         >  1:1576563/1‑69 (MQ=255)
             gCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAg                                         >  1:1429078/1‑69 (MQ=255)
                    ttttGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCa                                  >  1:2768083/1‑69 (MQ=255)
                    ttttGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCa                                  >  1:444675/1‑69 (MQ=255)
                    ttttGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCa                                  >  1:2047778/1‑69 (MQ=255)
                     tttGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCa                                  >  1:524273/1‑68 (MQ=255)
                     tttGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCa                                  >  1:2465573/1‑68 (MQ=255)
                     tttGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCa                                  >  1:3348021/1‑68 (MQ=255)
                               gtCTACTTATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:275548/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTTCAGCAGCAGTTTGCGTACCgtgc                                              <  1:1656906/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCGGTTTGCGTACCgtgc                                              <  1:866466/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:703044/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:3362438/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:344682/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:3194649/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:3083744/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:686271/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:622698/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:837954/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:3026693/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:2838400/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:445520/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:1064725/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:2661399/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:2084287/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:201612/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:116292/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:1221646/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:1364732/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:2470409/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:1640075/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:250711/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:1784661/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:1959318/46‑1 (MQ=255)
                               gtCTACTCATCTGACGGCATTTGCAGCAGCAGATTGCGTACCgtgc                                              <  1:220163/46‑1 (MQ=255)
                                tCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTg                                     >  1:1833775/1‑54 (MQ=255)
                                tCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGt                      <  1:1859696/69‑1 (MQ=255)
                                tCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGt                      <  1:1670661/69‑1 (MQ=255)
                                tCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGt                      <  1:1568279/69‑1 (MQ=255)
                                tCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGt                      <  1:982743/69‑1 (MQ=255)
                                tCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGt                      <  1:990227/69‑1 (MQ=255)
                                tCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGt                      <  1:2713780/69‑1 (MQ=255)
                                tCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGt                      <  1:3241950/69‑1 (MQ=255)
                                tCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGt                      <  1:22596/69‑1 (MQ=255)
                                tCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGt                      <  1:296461/69‑1 (MQ=255)
                                tCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGt                      <  1:2984348/69‑1 (MQ=255)
                                tCTACTCATCTGACGGCATTTGCAGCAGCAGTTTGCGTACCgtgc                                              <  1:2024857/45‑1 (MQ=255)
                                  tACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTa                                                    >  1:1206401/1‑37 (MQ=255)
                                  tACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTa                     >  1:2508605/1‑68 (MQ=255)
                                  tACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTa                     >  1:864082/1‑68 (MQ=255)
                                  tACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTa                     >  1:1646097/1‑68 (MQ=255)
                                  tACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTa                     >  1:2906359/1‑68 (MQ=255)
                                  tACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTa                     >  1:1821515/1‑68 (MQ=255)
                                  tACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTa                     >  1:2341992/1‑68 (MQ=255)
                                  tACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTa                     >  1:29998/1‑68 (MQ=255)
                                      cATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACg                  <  1:2878152/67‑1 (MQ=255)
                                         cTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGt                      <  1:749328/60‑1 (MQ=255)
                                             cGGCATTTGCAGCAGCAGTTTGCGTACCGTGCCGAAGCGTGCCATGTAACGCCTGTAACCAATTGaaa          >  1:2989319/1‑66 (MQ=255)
                                             cGGCATTTGCAGCAGCAGTTTGCGTACCGTGCCGAAGCGTGCCATGTAACGCCTGTAACCAATTGaaa          >  1:2590849/1‑66 (MQ=255)
                                             cGGCATTTGCAGCAGCAGTTTGCGTACCGTGCCGAAGCGTGCCATGTAACGCCTGTAACCAATTGaaa          >  1:972473/1‑66 (MQ=255)
                                             cGGCATTTGCAGCAGCAGTTTGCGTACCGTGCCGAAGCGTGCCATGTAACGCCTGTAACCAATTGaaa          >  1:841889/1‑66 (MQ=255)
                                             cGGCATTTGCAGCAGCAGTTTGCGTACCGTGCCGAAGCGTGCCATGTAACGCCTGTAACCAATTGaaa          >  1:805763/1‑66 (MQ=255)
                                             cGGCATTTGCAGCAGCAGTTTGCGTACCGTGCCGAAGCGTGCCATGTAACGCCTGTAACCAATTGaaa          >  1:2538318/1‑66 (MQ=255)
                                             cGGCATTTGCAGCAGCAGTTTGCGTACCGTGCCGAAGCGTGCCATGTAACGCCTGTAACCAATTGaaa          >  1:760794/1‑66 (MQ=255)
                                             cGGCATTTGCAGCAGCAGTTTGCGTACCGTGCCGAAGCGTGCCATGTAACGCCTGTAACCAATTGaaa          >  1:636171/1‑66 (MQ=255)
                                             cGGCATTTGCAGCAGCAGTTTGCGTACCGTGCCGAAGCGTGCCATGTAACGCCTGTAACCAATTGaaa          >  1:1865756/1‑66 (MQ=255)
                                             cGGCATTTGCAGCAGCAGTTTGCGTACCGTGCCGAAGCGTGCCATGTAACGCCTGTAACCAATTGaaa          >  1:377912/1‑66 (MQ=255)
                                              ggaaTTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATAtt        >  1:3299804/4‑69 (MQ=255)
                                              ggCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATAtt        >  1:140838/1‑69 (MQ=255)
                                              ggCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATAtt        >  1:733408/1‑69 (MQ=255)
                                              ggCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATAtt        >  1:2546379/1‑69 (MQ=255)
                                              ggCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATAtt        >  1:2574907/1‑69 (MQ=255)
                                              ggCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATAtt        >  1:523636/1‑69 (MQ=255)
                                                    tGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTg  <  1:854984/69‑1 (MQ=255)
                                                    tGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTg  <  1:613334/69‑1 (MQ=255)
                                                    tGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTg  <  1:2549188/69‑1 (MQ=255)
                                                       |                                                                 
CTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTG  >  minE/131494‑131614

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: