Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,090,837 A→G 45.0% F288F (TTT→TTC sufB ← component of SufBCD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,090,8370AG45.0% 7.4 / 48.2 60F288F (TTT→TTCsufBcomponent of SufBCD complex
Reads supporting (aligned to +/- strand):  ref base A (6/27);  new base G (10/17);  total (16/44)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

TTCGCCTTCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCCACCACTGC  >  minE/1090771‑1090902
                                                                  |                                                                 
ttCGCCTTCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGaa                                                                  >  1:2554292/1‑68 (MQ=255)
ttCGCCTTCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGa                                                                   >  1:2558226/1‑67 (MQ=255)
      ttCGGACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:1503515/69‑1 (MQ=255)
      ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:1180195/69‑1 (MQ=255)
      ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:2398891/69‑1 (MQ=255)
      ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:1956366/69‑1 (MQ=255)
      ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:1782099/69‑1 (MQ=255)
      ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:2496266/69‑1 (MQ=255)
      ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:256180/69‑1 (MQ=255)
      ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:2733074/69‑1 (MQ=255)
      ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:2750429/69‑1 (MQ=255)
      ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:2997663/69‑1 (MQ=255)
      ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:528922/69‑1 (MQ=255)
      ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:661977/69‑1 (MQ=255)
      ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:80149/69‑1 (MQ=255)
       tCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:812721/68‑1 (MQ=255)
       tCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:2182396/68‑1 (MQ=255)
           acaAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCTGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:1479191/64‑1 (MQ=255)
                 cACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:1888437/58‑1 (MQ=255)
                   cGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGaa                                              <  1:1792398/69‑1 (MQ=255)
                   cGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGaa                                              <  1:1811377/69‑1 (MQ=255)
                    gCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGaa                                              <  1:2687044/68‑1 (MQ=255)
                            aCGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGa                                               >  1:1777241/1‑59 (MQ=255)
                            aCGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCAcc                                     >  1:1638752/1‑69 (MQ=255)
                            aCGAAGGTGAGAATACCGCCGGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCAcc                                     >  1:435175/1‑69 (MQ=255)
                                       aaTACCGCCGGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGtttt                           <  1:995354/68‑1 (MQ=255)
                                       aaTACCGCCGGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGtttt                           <  1:407400/68‑1 (MQ=255)
                                         tACCGCCGGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTGTTTa                         >  1:2737171/1‑68 (MQ=255)
                                          aCCGCCGGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATg                       <  1:1228184/69‑1 (MQ=255)
                                           ccgccgGTGTTGTTATCGCCAGGACACCAGTTTTGTACCGTGGAat                                             <  1:3368327/46‑1 (MQ=255)
                                           ccgccgGTGTTGTTATCGCCAGGAAACCAGTTTTg                                                        >  1:780659/1‑35 (MQ=255)
                                           ccgccgGTGTTGTTATCGCCAGGAAACCAGTTTTg                                                        >  1:2858103/1‑35 (MQ=255)
                                           ccgccgGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAat                                             <  1:199535/46‑1 (MQ=255)
                                                ggTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:159179/69‑1 (MQ=255)
                                                ggTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:652572/69‑1 (MQ=255)
                                                ggTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:3185818/69‑1 (MQ=255)
                                                ggTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:857410/69‑1 (MQ=255)
                                                ggTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:2829190/69‑1 (MQ=255)
                                                ggTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:1431013/69‑1 (MQ=255)
                                                ggTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:106760/69‑1 (MQ=255)
                                                 gtgtTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:880721/68‑1 (MQ=255)
                                                  tgttgtTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGAtgatga               <  1:1292803/69‑1 (MQ=255)
                                                    ttgttATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATg                       <  1:1317569/59‑1 (MQ=255)
                                                    ttgttATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:874287/65‑1 (MQ=255)
                                                    ttgttATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:2952533/65‑1 (MQ=255)
                                                    ttgttATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:748279/65‑1 (MQ=255)
                                                    ttgttATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGTATGATGACt             <  1:1184301/69‑1 (MQ=255)
                                                    ttgttATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACt             <  1:2679018/69‑1 (MQ=255)
                                                    ttgttATCGCCAGGAAACAAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:2564934/65‑1 (MQ=255)
                                                     tgttATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:2907149/64‑1 (MQ=255)
                                                     tgttATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:2582539/64‑1 (MQ=255)
                                                      gttATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTc           <  1:1278360/69‑1 (MQ=255)
                                                               aGGGAACCAATTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTccacc       >  1:662360/1‑64 (MQ=255)
                                                               aGGGAACCAATTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCcaccaccg   >  1:2529995/1‑66 (MQ=255)
                                                               aGGGAACCAATTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCcaccaccg   >  1:3240794/1‑66 (MQ=255)
                                                               aGGGAACCAATTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCCACCACCGc  >  1:1722107/1‑69 (MQ=255)
                                                               aGGGAACCAATTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCCACCACCGc  >  1:735165/1‑69 (MQ=255)
                                                               aGGGAACCAATTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCCACCACCGc  >  1:3251393/1‑69 (MQ=255)
                                                               aGGGAACCAATTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCCACCACCGc  >  1:3143555/1‑69 (MQ=255)
                                                               aGGGAACCAATTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCCACCACCGc  >  1:3012062/1‑69 (MQ=255)
                                                               aGGGAACCAATTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCCACCACCGc  >  1:2838428/1‑69 (MQ=255)
                                                               aGGGAACCAATTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCCACCACCGc  >  1:1406453/1‑69 (MQ=255)
                                                                  |                                                                 
TTCGCCTTCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCCACCACTGC  >  minE/1090771‑1090902

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