Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,090,843 G→A 45.0% N286N (AAC→AAT sufB ← component of SufBCD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,090,8430GA45.0% 10.8 / 46.4 60N286N (AAC→AATsufBcomponent of SufBCD complex
Reads supporting (aligned to +/- strand):  ref base G (6/27);  new base A (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

TTCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCCACCACTGC  >  minE/1090777‑1090902
                                                                  |                                                           
ttCGGACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:1503515/69‑1 (MQ=255)
ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:1180195/69‑1 (MQ=255)
ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:80149/69‑1 (MQ=255)
ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:661977/69‑1 (MQ=255)
ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:528922/69‑1 (MQ=255)
ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:2997663/69‑1 (MQ=255)
ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:2750429/69‑1 (MQ=255)
ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:2733074/69‑1 (MQ=255)
ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:256180/69‑1 (MQ=255)
ttCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGGAACCAAtt                                                           <  1:2496266/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)
 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:652572/69‑1 (MQ=255)
                                          ggTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:1431013/69‑1 (MQ=255)
                                          ggTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:106760/69‑1 (MQ=255)
                                          ggTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:3185818/69‑1 (MQ=255)
                                          ggTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:159179/69‑1 (MQ=255)
                                          ggTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:857410/69‑1 (MQ=255)
                                          ggTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:2829190/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:2952533/65‑1 (MQ=255)
                                              ttgttATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGgatgat                 <  1:874287/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:3240794/1‑66 (MQ=255)
                                                         aGGGAACCAATTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCcaccaccg   >  1:2529995/1‑66 (MQ=255)
                                                         aGGGAACCAATTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCCACCACCGc  >  1:735165/1‑69 (MQ=255)
                                                         aGGGAACCAATTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCCACCACCGc  >  1:3251393/1‑69 (MQ=255)
                                                         aGGGAACCAATTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCCACCACCGc  >  1:1406453/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:1722107/1‑69 (MQ=255)
                                                                  |                                                           
TTCGCACAAAGCACGCTTGGTGACGAAGTTGAGAATACCGCCGGTGTTGTTATCGCCAGGAAACCAGTTTTGTACCGTGGAATATTTCACCTCGGCGTTTTTATGGATGATGACTTCCACCACTGC  >  minE/1090777‑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: