Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,449,128 A→G 51.2% Q171Q (CAA→CAG 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,1280AG51.2% ‑1.1 / 53.0 62Q171Q (CAA→CAGyifBpredicted bifunctional protein, enzyme and transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (17/13);  new base G (14/17);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

AAAGATAACGAAGATGAAGTGGGGCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCC  >  minE/2449063‑2449194
                                                                 |                                                                  
aaaGATAACGAAGATGAAGTGGGGCTAATTAACGGTGAAGGTTGCCTGATAGCCGATCATCTGCAAGct                                                                 <  1:739405/69‑1 (MQ=255)
aaaGATAACGAAGATGAAGTGGGGCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGct                                                                 <  1:2589161/69‑1 (MQ=255)
aaaGATAACGAAGATGAAGTGGGGCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGct                                                                 <  1:353560/69‑1 (MQ=255)
aaaGATAACGAAGATGAAGTGGGGCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGct                                                                 <  1:13098/69‑1 (MQ=255)
  aGATAACGAAGATGAAGTGGGGCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGct                                                                 <  1:1043847/67‑1 (MQ=255)
                    tgggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                             <  1:1041759/68‑1 (MQ=255)
                    ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                             <  1:2097779/69‑1 (MQ=255)
                    ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                             <  1:2372088/69‑1 (MQ=255)
                    ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                             <  1:2585365/69‑1 (MQ=255)
                    ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                             <  1:1856632/69‑1 (MQ=255)
                    ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                             <  1:2725655/69‑1 (MQ=255)
                    ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                             <  1:1665840/69‑1 (MQ=255)
                    ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                             <  1:3306960/69‑1 (MQ=255)
                    ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                             <  1:1296537/69‑1 (MQ=255)
                    ggggCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAgg                                             <  1:551541/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:282037/1‑69 (MQ=255)
                           aTTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCAc                                      >  1:134519/1‑69 (MQ=255)
                           aTTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCAc                                      >  1:185002/1‑69 (MQ=255)
                           aTTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCAc                                      >  1:533688/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:201694/1‑68 (MQ=255)
                                             cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                     >  1:914199/1‑68 (MQ=255)
                                             cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                     >  1:1231485/1‑68 (MQ=255)
                                             cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                     >  1:1279500/1‑68 (MQ=255)
                                             cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                     >  1:1337024/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:2698910/1‑68 (MQ=255)
                                             cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                     >  1:2035343/1‑68 (MQ=255)
                                             cTGATAGCCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                     >  1:2187993/1‑68 (MQ=255)
                                             cTGATAACCGATCATCTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGaa                     >  1:1991/1‑68 (MQ=255)
                                               gATAGCCGATCATCTGCATGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                  <  1:2002059/69‑1 (MQ=255)
                                               gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                  <  1:2816202/69‑1 (MQ=255)
                                               gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                  <  1:263278/69‑1 (MQ=255)
                                               gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                  <  1:481784/69‑1 (MQ=255)
                                               gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                  <  1:3062869/69‑1 (MQ=255)
                                               gATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAAcc                  <  1:3292829/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:2251165/1‑69 (MQ=255)
                                                     cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa            >  1:2586502/1‑69 (MQ=255)
                                                     cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa            >  1:674143/1‑69 (MQ=255)
                                                     cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa            >  1:2918755/1‑69 (MQ=255)
                                                     cGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGa            >  1:1980186/1‑69 (MQ=255)
                                                          atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg       <  1:2537010/69‑1 (MQ=255)
                                                          atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg       <  1:3057307/69‑1 (MQ=255)
                                                          atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg       <  1:1299031/69‑1 (MQ=255)
                                                          atcTGCAGGCTGTCTGCGCGTTTCTGGAAGGTAAGCACGCTCTCGAACACCCGAAACCAACTAATGCAg       <  1:2928382/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)
                                                                 |                                                                  
AAAGATAACGAAGATGAAGTGGGGCTAATTAACGGTGAAGGATGCCTGATAGCCGATCATCTGCAAGCTGTCTGTGCGTTTCTGGAAGGTAAGCACGCTCTCGAACGCCCGAAACCAACTGATGCAGTATCC  >  minE/2449063‑2449194

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: