Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,449,137 T→C 54.4% C174C (TGT→TGC yifB → predicted bifunctional protein, enzyme and transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,449,1370TC54.4% 7.9 / 44.6 57C174C (TGT→TGCyifBpredicted bifunctional protein, enzyme and transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (17/9);  new base C (14/17);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGC  >  minE/2449083‑2449201
                                                      |                                                                
tgggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                                    <  1:1041759/68‑1 (MQ=255)
ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                                    <  1:2372088/69‑1 (MQ=255)
ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                                    <  1:551541/69‑1 (MQ=255)
ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                                    <  1:3306960/69‑1 (MQ=255)
ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                                    <  1:2725655/69‑1 (MQ=255)
ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                                    <  1:2585365/69‑1 (MQ=255)
ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                                    <  1:2097779/69‑1 (MQ=255)
ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                                    <  1:1856632/69‑1 (MQ=255)
ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                                    <  1:1665840/69‑1 (MQ=255)
ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                                    <  1:1296537/69‑1 (MQ=255)
ggggCTAATTAACGGTGAAGAATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                                    <  1:1599631/69‑1 (MQ=255)
ggggCTAATAAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                                    <  1:761008/69‑1 (MQ=255)
       aTTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCAc                                             >  1:1356680/1‑69 (MQ=255)
       aTTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCAc                                             >  1:134519/1‑69 (MQ=255)
       aTTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCAc                                             >  1:533688/1‑69 (MQ=255)
       aTTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCAc                                             >  1:282037/1‑69 (MQ=255)
       aTTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCAc                                             >  1:185002/1‑69 (MQ=255)
       aTTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCAc                                             >  1:1887038/1‑69 (MQ=255)
         tAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACg                                            <  1:935896/68‑1 (MQ=255)
              gTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCg                                      <  1:1389444/69‑1 (MQ=255)
                      tGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGccc                               >  1:2997321/1‑68 (MQ=255)
                         cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACAccc                               >  1:1152753/1‑65 (MQ=255)
                         cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                            >  1:914199/1‑68 (MQ=255)
                         cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                            >  1:201694/1‑68 (MQ=255)
                         cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                            >  1:2035343/1‑68 (MQ=255)
                         cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                            >  1:1337024/1‑68 (MQ=255)
                         cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                            >  1:1231485/1‑68 (MQ=255)
                         cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                            >  1:2187993/1‑68 (MQ=255)
                         cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                            >  1:3219560/1‑68 (MQ=255)
                         cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                            >  1:3096494/1‑68 (MQ=255)
                         cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                            >  1:3027641/1‑68 (MQ=255)
                         cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                            >  1:2908961/1‑68 (MQ=255)
                         cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                            >  1:1279500/1‑68 (MQ=255)
                         cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                            >  1:2698910/1‑68 (MQ=255)
                         cTGATAACCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                            >  1:1991/1‑68 (MQ=255)
                           gATAGCCGATCATCTGCATGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                         <  1:2002059/69‑1 (MQ=255)
                           gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                         <  1:263278/69‑1 (MQ=255)
                           gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                         <  1:2816202/69‑1 (MQ=255)
                           gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                         <  1:481784/69‑1 (MQ=255)
                           gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                         <  1:3292829/69‑1 (MQ=255)
                           gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                         <  1:3062869/69‑1 (MQ=255)
                               gCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACt                     >  1:2115206/1‑69 (MQ=255)
                               gCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACt                     >  1:3370659/1‑69 (MQ=255)
                                 cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGctct                                        >  1:1021962/1‑48 (MQ=255)
                                 cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTg                    >  1:17899/1‑68 (MQ=255)
                                 cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                   >  1:1980186/1‑69 (MQ=255)
                                 cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                   >  1:674143/1‑69 (MQ=255)
                                 cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                   >  1:2251165/1‑69 (MQ=255)
                                 cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                   >  1:2918755/1‑69 (MQ=255)
                                 cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa                   >  1:2586502/1‑69 (MQ=255)
                                      atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg              <  1:3057307/69‑1 (MQ=255)
                                      atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg              <  1:2537010/69‑1 (MQ=255)
                                      atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg              <  1:2928382/69‑1 (MQ=255)
                                      atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg              <  1:1299031/69‑1 (MQ=255)
                                      atcTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTATCGAACGCCCGAAACCAACTGATGCAg              <  1:3109224/69‑1 (MQ=255)
                                       tcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg              <  1:1299442/68‑1 (MQ=255)
                                           cAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATcc         >  1:2246788/1‑69 (MQ=255)
                                                    tgtgCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                         <  1:2669019/44‑1 (MQ=255)
                                                    tgtgCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGgcgc  <  1:2165061/67‑1 (MQ=255)
                                                    tgtgCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGgcgc  <  1:1473579/67‑1 (MQ=255)
                                                      |                                                                
GGGGCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCCCGGGCGC  >  minE/2449083‑2449201

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: