Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 994,533 T→C 24.6% P126P (CCA→CCG ynfK ← predicted dethiobiotin synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE994,5330TC24.6% 96.0 / 27.6 81P126P (CCA→CCGynfKpredicted dethiobiotin synthetase
Reads supporting (aligned to +/- strand):  ref base T (23/38);  new base C (8/12);  total (31/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

GACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTT  >  minE/994485‑994599
                                                |                                                                  
gACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATtcatca                                                 <  1:856344/68‑1 (MQ=255)
gACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATtcatca                                                 <  1:2299297/68‑1 (MQ=255)
gACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATtcatca                                                 <  1:3029391/68‑1 (MQ=255)
gACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATtcatca                                                 <  1:53222/68‑1 (MQ=255)
gACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATtcatca                                                 <  1:420709/68‑1 (MQ=255)
          aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCa                                       >  1:2066386/1‑68 (MQ=255)
          aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCa                                       >  1:2568230/1‑68 (MQ=255)
          aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAg                                      >  1:2061862/1‑69 (MQ=255)
          aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAg                                      >  1:2168533/1‑69 (MQ=255)
          aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAg                                      >  1:2309264/1‑69 (MQ=255)
          aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAg                                      >  1:292822/1‑69 (MQ=255)
          aaCACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAg                                      >  1:3110280/1‑69 (MQ=255)
              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                                  <  1:2813521/69‑1 (MQ=255)
              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                                  <  1:885758/69‑1 (MQ=255)
              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                                  <  1:430748/69‑1 (MQ=255)
              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                                  <  1:3116580/69‑1 (MQ=255)
              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                                  <  1:2184704/69‑1 (MQ=255)
              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                                  <  1:3249810/69‑1 (MQ=255)
              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                                  <  1:3400536/69‑1 (MQ=255)
              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                                  <  1:380011/69‑1 (MQ=255)
              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                                  <  1:1172492/69‑1 (MQ=255)
              ccGGCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAgccgc                                  <  1:1443740/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:1322687/1‑69 (MQ=255)
                 gCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAg                               >  1:969335/1‑69 (MQ=255)
                 gCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAg                               >  1:2735370/1‑69 (MQ=255)
                 gCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAg                               >  1:2211601/1‑69 (MQ=255)
                 gCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAg                               >  1:2796535/1‑69 (MQ=255)
                 gCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAg                               >  1:888375/1‑69 (MQ=255)
                 gCAGTTGTTCCTGCACCACCCATTCAGAGAGCGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAg                               >  1:2431597/1‑69 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          >  1:1236752/1‑68 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          >  1:693465/1‑68 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          <  1:2986783/68‑1 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          >  1:568410/1‑68 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          >  1:1218510/1‑68 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          <  1:3179276/68‑1 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          <  1:3324458/68‑1 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          <  1:3402971/68‑1 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          <  1:643021/68‑1 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          >  1:1479079/1‑68 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          >  1:979975/1‑68 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          <  1:243843/68‑1 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          <  1:1481723/68‑1 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          <  1:166063/68‑1 (MQ=255)
                       gttCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          >  1:2210646/1‑68 (MQ=255)
                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa                <  1:1561544/69‑1 (MQ=255)
                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa                <  1:1804998/69‑1 (MQ=255)
                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa                <  1:822771/69‑1 (MQ=255)
                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa                <  1:3300040/69‑1 (MQ=255)
                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa                <  1:2291082/69‑1 (MQ=255)
                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa                <  1:86436/69‑1 (MQ=255)
                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa                <  1:3017406/69‑1 (MQ=255)
                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa                <  1:3010130/69‑1 (MQ=255)
                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa                <  1:1306719/69‑1 (MQ=255)
                                cTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa                <  1:852809/69‑1 (MQ=255)
                                 tACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGa                <  1:772377/68‑1 (MQ=255)
                                        tCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCGAGTCCCTTCTACCACGACATGATCg        <  1:2819768/69‑1 (MQ=255)
                                        tCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCt                          <  1:1608954/51‑1 (MQ=255)
                                        tCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGTCATGATCg        <  1:3350990/69‑1 (MQ=255)
                                        tCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCg        <  1:768638/69‑1 (MQ=255)
                                        tCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCg        <  1:2457329/69‑1 (MQ=255)
                                        tCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCg        <  1:2885749/69‑1 (MQ=255)
                                        tCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCg        <  1:2485926/69‑1 (MQ=255)
                                        tCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCg        <  1:2447117/69‑1 (MQ=255)
                                        tCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCg        <  1:2020308/69‑1 (MQ=255)
                                        tCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCg        <  1:1319306/69‑1 (MQ=255)
                                        tCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCg        <  1:1286883/69‑1 (MQ=255)
                                        tCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCg        <  1:1274120/69‑1 (MQ=255)
                                        tCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGACGCCAGTCCCTTCTACCACGACATGATCg        <  1:110409/69‑1 (MQ=255)
                                         cAGAGAGTGGACGCTAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGa       >  1:1136023/1‑69 (MQ=255)
                                         cAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGa       >  1:2591225/1‑69 (MQ=255)
                                         cAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGa       >  1:3211115/1‑69 (MQ=255)
                                         cAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGa       >  1:1030537/1‑69 (MQ=255)
                                         cAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGa       >  1:2399133/1‑69 (MQ=255)
                                             gagTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACttt   <  1:508828/69‑1 (MQ=255)
                                              agTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACtttt  >  1:2566498/1‑69 (MQ=255)
                                              agTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACtttt  >  1:2164574/1‑69 (MQ=255)
                                              agTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACtttt  >  1:1697840/1‑69 (MQ=255)
                                              agTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACtttt  >  1:991264/1‑69 (MQ=255)
                                                |                                                                  
GACAACCATCAACACCGGCAGTTGTTCCTGCACTACCCATTCAGAGAGTGGACGCAAATCATTCATCAGACTGCGCCAGCCGCCAGTCCCTTCTACCACGACATGATCGACTTTT  >  minE/994485‑994599

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: