Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 994,518 T→C 22.9% V131V (GTA→GTG ynfK ← predicted dethiobiotin synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE994,5180TC22.9% 112.4 / 26.9 87V131V (GTA→GTGynfKpredicted dethiobiotin synthetase
Reads supporting (aligned to +/- strand):  ref base T (31/36);  new base C (8/12);  total (39/48)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

AGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGA  >  minE/994453‑994585
                                                                 |                                                                   
aGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTAcc                                                                  <  1:2825644/69‑1 (MQ=255)
aGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTAcc                                                                  <  1:997871/69‑1 (MQ=255)
aGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTAcc                                                                  <  1:3407663/69‑1 (MQ=255)
aGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTAcc                                                                  <  1:2633243/69‑1 (MQ=255)
aGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTAcc                                                                  <  1:2515559/69‑1 (MQ=255)
aGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTAcc                                                                  <  1:2490808/69‑1 (MQ=255)
aGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTAcc                                                                  <  1:2483918/69‑1 (MQ=255)
aGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTAcc                                                                  <  1:483206/69‑1 (MQ=255)
aGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTAcc                                                                  <  1:2442205/69‑1 (MQ=255)
aGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTAcc                                                                  <  1:411305/69‑1 (MQ=255)
aGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTAcc                                                                  <  1:2408520/69‑1 (MQ=255)
   gCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCa                                                                >  1:209746/1‑68 (MQ=255)
   gCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCa                                                                >  1:1128018/1‑68 (MQ=255)
   gCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCa                                                                >  1:3213035/1‑68 (MQ=255)
   gCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCa                                                                >  1:2406418/1‑68 (MQ=255)
   gCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCa                                                                >  1:3154263/1‑68 (MQ=255)
   gCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCa                                                                >  1:2980486/1‑68 (MQ=255)
   gCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCa                                                                >  1:290491/1‑68 (MQ=255)
   gCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCa                                                                >  1:29077/1‑68 (MQ=255)
   gCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCa                                                                >  1:2003214/1‑68 (MQ=255)
   gCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCa                                                                >  1:682839/1‑68 (MQ=255)
   gCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCa                                                                >  1:1761244/1‑68 (MQ=255)
   gCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCa                                                                >  1:531521/1‑68 (MQ=255)
   gCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACAGGCAGTTGTTCCTGCACTACCCa                                                                >  1:1835891/1‑68 (MQ=255)
       ggTTAATGCAACCTTCCTGAATACCGACACCCATCAACACCGGCAGTTGTTCCTGCACTACCCAt                                                               >  1:265712/1‑65 (MQ=255)
       ggTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCAtt                                                              >  1:26890/1‑66 (MQ=255)
       ggTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCAtt                                                              >  1:1990178/1‑66 (MQ=255)
       ggTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCAtt                                                              >  1:1205478/1‑66 (MQ=255)
        gTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCaga                                                          <  1:1194651/69‑1 (MQ=255)
                           aTACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTAcc                                                                  <  1:3331623/42‑1 (MQ=255)
                                gACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATtcatca                                   <  1:3029391/68‑1 (MQ=255)
                                gACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATtcatca                                   <  1:856344/68‑1 (MQ=255)
                                gACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATtcatca                                   <  1:53222/68‑1 (MQ=255)
                                gACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATtcatca                                   <  1:420709/68‑1 (MQ=255)
                                gACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATtcatca                                   <  1:2299297/68‑1 (MQ=255)
                                          aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCa                         >  1:2568230/1‑68 (MQ=255)
                                          aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCa                         >  1:2066386/1‑68 (MQ=255)
                                          aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAg                        >  1:292822/1‑69 (MQ=255)
                                          aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAg                        >  1:2309264/1‑69 (MQ=255)
                                          aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAg                        >  1:2168533/1‑69 (MQ=255)
                                          aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAg                        >  1:2061862/1‑69 (MQ=255)
                                          aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAg                        >  1:3110280/1‑69 (MQ=255)
                                              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                    <  1:885758/69‑1 (MQ=255)
                                              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                    <  1:1172492/69‑1 (MQ=255)
                                              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                    <  1:430748/69‑1 (MQ=255)
                                              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                    <  1:3116580/69‑1 (MQ=255)
                                              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                    <  1:380011/69‑1 (MQ=255)
                                              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                    <  1:3249810/69‑1 (MQ=255)
                                              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                    <  1:3400536/69‑1 (MQ=255)
                                              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                    <  1:2813521/69‑1 (MQ=255)
                                              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                    <  1:1443740/69‑1 (MQ=255)
                                              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                    <  1:2184704/69‑1 (MQ=255)
                                              ccGGCAGTGGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                    <  1:1781034/69‑1 (MQ=255)
                                               cGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                    <  1:2720543/68‑1 (MQ=255)
                                                 gCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCa                  >  1:1229926/1‑68 (MQ=255)
                                                 gCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAg                 >  1:2796535/1‑69 (MQ=255)
                                                 gCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAg                 >  1:969335/1‑69 (MQ=255)
                                                 gCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAg                 >  1:2735370/1‑69 (MQ=255)
                                                 gCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAg                 >  1:888375/1‑69 (MQ=255)
                                                 gCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAg                 >  1:2431597/1‑69 (MQ=255)
                                                 gCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAg                 >  1:2211601/1‑69 (MQ=255)
                                                 gCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAg                 >  1:1322687/1‑69 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            <  1:1481723/68‑1 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            >  1:1218510/1‑68 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            >  1:1479079/1‑68 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            >  1:568410/1‑68 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            <  1:643021/68‑1 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            >  1:693465/1‑68 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            >  1:1236752/1‑68 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            >  1:979975/1‑68 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            <  1:166063/68‑1 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            <  1:2986783/68‑1 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            <  1:243843/68‑1 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            <  1:3179276/68‑1 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            <  1:3324458/68‑1 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            >  1:2210646/1‑68 (MQ=255)
                                                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt            <  1:3402971/68‑1 (MQ=255)
                                                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa  <  1:1306719/69‑1 (MQ=255)
                                                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa  <  1:3300040/69‑1 (MQ=255)
                                                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa  <  1:822771/69‑1 (MQ=255)
                                                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa  <  1:852809/69‑1 (MQ=255)
                                                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa  <  1:2291082/69‑1 (MQ=255)
                                                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa  <  1:86436/69‑1 (MQ=255)
                                                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa  <  1:1804998/69‑1 (MQ=255)
                                                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa  <  1:3017406/69‑1 (MQ=255)
                                                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa  <  1:3010130/69‑1 (MQ=255)
                                                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa  <  1:1561544/69‑1 (MQ=255)
                                                                 tACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa  <  1:772377/68‑1 (MQ=255)
                                                                 |                                                                   
AGTGCATGGTTAATGCAACCTTCCTGAATACCGACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGA  >  minE/994453‑994585

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: