Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,146,056 G→C 23.8% P46A (CCG→GCG)  chbF ← cryptic phospho‑beta‑glucosidase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,146,0560GC23.8% 101.3 / 28.0 63P46A (CCG→GCG) chbFcryptic phospho‑beta‑glucosidase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (28/20);  new base C (8/7);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

AAGCTTCATCGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAA  >  minE/1145988‑1146121
                                                                    |                                                                 
aaGCTTCATCGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCgg                                                                   <  1:2755424/69‑1 (MQ=255)
aaGCTTCATCGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCgg                                                                   <  1:1382758/69‑1 (MQ=255)
aaGCTTCATCGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCgg                                                                   <  1:559046/69‑1 (MQ=255)
aaGCTTCATCGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCgg                                                                   <  1:1554428/69‑1 (MQ=255)
     tCATCGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTac                                                              <  1:2284515/69‑1 (MQ=255)
     tCATCGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTac                                                              <  1:3254858/69‑1 (MQ=255)
      cATCGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTac                                                              <  1:1269294/68‑1 (MQ=255)
        tCGGGACGCCATCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTacc                                                             <  1:2534330/67‑1 (MQ=255)
        tCGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTacc                                                             <  1:874969/67‑1 (MQ=255)
        tCGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTacc                                                             <  1:2180966/67‑1 (MQ=255)
        tCGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTacc                                                             <  1:2849793/67‑1 (MQ=255)
        tCGGGACGCCAGCGTTATCAATCATCCGTTAGCAGAGATCGAAAATAATATCCAGTTTCGGTTTacc                                                             <  1:2356756/67‑1 (MQ=255)
         cGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTaccac                                                           <  1:1120217/68‑1 (MQ=255)
         cGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTaccac                                                           <  1:3247617/68‑1 (MQ=255)
         cGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTaccac                                                           <  1:3435157/68‑1 (MQ=255)
             aCGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCgg                                                                   <  1:2015365/56‑1 (MQ=255)
                        atcaatcaTCCGTTGGCAGAGATCAAAAATAATATCCAGTTTCGCtt                                                                 >  1:1343573/1‑47 (MQ=39)
                        atcaatcaTCCGTTGGCAGAGATCAAAAATAATATCCAGTTTCGCtt                                                                 >  1:2796999/1‑47 (MQ=39)
                        atcaatcaTCCGTTGGCAGAGATCAAAAATAATATCCAGTTTCGCtt                                                                 >  1:1339803/1‑47 (MQ=39)
                        atcaatcaTCCGTTGGCAGAGATCAAAAATAATATCCAGTTTCGCtt                                                                 >  1:2539531/1‑47 (MQ=39)
                        atcaatcaTCCGTTGGCAGAGATCAAAAATAATATCCAGTTTCGCtt                                                                 >  1:1271958/1‑47 (MQ=39)
                        atcaatcaTCCGTTGGCAGAGATCAAAAATAATATCCAGTTTCGCtt                                                                 >  1:2138626/1‑47 (MQ=39)
                        atcaatcaTCCGTTGGCAGAGATCAAAAATAATATCCAGTTTCGCtt                                                                 >  1:1178160/1‑47 (MQ=39)
                        atcaatcaTCCGTTGGCAGAGATCAAAAATAATATCCAGTTTCGCtt                                                                 >  1:2258682/1‑47 (MQ=39)
                               atcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGt                                                                  >  1:2301189/1‑39 (MQ=255)
                               atcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATcc                                                >  1:3101218/1‑57 (MQ=255)
                               atcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATcc                                                >  1:2379098/1‑57 (MQ=255)
                               atcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATcc                                                >  1:691004/1‑57 (MQ=255)
                               atcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATcc                                                >  1:2519206/1‑57 (MQ=255)
                               atcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATcc                                                >  1:1779274/1‑57 (MQ=255)
                               atcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATcc                                                >  1:1918672/1‑57 (MQ=255)
                               atcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATcc                                                >  1:1414541/1‑57 (MQ=255)
                               atcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATc                                                 >  1:1935402/1‑56 (MQ=255)
                                            gaTCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                       <  1:1210714/69‑1 (MQ=255)
                                            gaTCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                       <  1:2287865/69‑1 (MQ=255)
                                            gaTCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                       <  1:2968027/69‑1 (MQ=255)
                                            gaTCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                       <  1:2949591/69‑1 (MQ=255)
                                            gaTCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                       <  1:367670/69‑1 (MQ=255)
                                            gaTCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                       <  1:1132533/69‑1 (MQ=255)
                                            gaTCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                       <  1:2398913/69‑1 (MQ=255)
                                            gaTCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                       <  1:994138/69‑1 (MQ=255)
                                                      aataTCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGAt              >  1:2501181/1‑68 (MQ=255)
                                                      aataTCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa             >  1:733005/1‑69 (MQ=255)
                                                      aataTCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa             >  1:1752485/1‑69 (MQ=255)
                                                      aataTCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa             >  1:399642/1‑69 (MQ=255)
                                                      aataTCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa             >  1:910758/1‑69 (MQ=255)
                                                      aataTCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa             >  1:976854/1‑69 (MQ=255)
                                                      aataTCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa             >  1:989515/1‑69 (MQ=255)
                                                      aataTCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa             >  1:99993/1‑69 (MQ=255)
                                                      aataTCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa             >  1:2048190/1‑69 (MQ=255)
                                                      aataTCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa             >  1:1863516/1‑69 (MQ=255)
                                                      aataTCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa             >  1:2656666/1‑69 (MQ=255)
                                                       atatCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                       <  1:2653804/58‑1 (MQ=255)
                                                       atatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATaa            <  1:2243945/69‑1 (MQ=255)
                                                       atatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATaa            <  1:99240/69‑1 (MQ=255)
                                                       atatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATaa            <  1:2046492/69‑1 (MQ=255)
                                                       atatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATaa            <  1:880895/69‑1 (MQ=255)
                                                       atatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATaa            <  1:40110/69‑1 (MQ=255)
                                                       atatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATaa            <  1:2740365/69‑1 (MQ=255)
                                                       atatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATaa            <  1:364980/69‑1 (MQ=255)
                                                                ttCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaata   >  1:1993855/1‑69 (MQ=255)
                                                                 tCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaataa  >  1:2739879/1‑69 (MQ=255)
                                                                 tCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaataa  >  1:625462/1‑69 (MQ=255)
                                                                 tCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaataa  >  1:2610052/1‑69 (MQ=255)
                                                                 tCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaataa  >  1:369298/1‑69 (MQ=255)
                                                                 tCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaataa  >  1:1114736/1‑69 (MQ=255)
                                                                 tCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaataa  >  1:1589724/1‑69 (MQ=255)
                                                                 tCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaataa  >  1:3306848/1‑69 (MQ=255)
                                                                    gTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa                       <  1:1816261/45‑1 (MQ=255)
                                                                    |                                                                 
AAGCTTCATCGGGACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGGTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAA  >  minE/1145988‑1146121

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: