Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 384,439 A→G 84.5% V81V (GTA→GTG ybdB → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE384,4390AG84.5% 154.0 / 13.5 89V81V (GTA→GTGybdBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (5/9);  new base G (24/51);  total (29/60)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGC  >  minE/384388‑384506
                                                   |                                                                   
gCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGc                                                    <  1:3299916/69‑1 (MQ=255)
gCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGc                                                    <  1:263966/69‑1 (MQ=255)
gCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGc                                                    <  1:194406/69‑1 (MQ=255)
gCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGc                                                    <  1:2821617/69‑1 (MQ=255)
gCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGc                                                    <  1:405691/69‑1 (MQ=255)
gCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGc                                                    <  1:1501636/69‑1 (MQ=255)
gCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGc                                                    <  1:1285579/69‑1 (MQ=255)
gCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGc                                                    <  1:1344206/69‑1 (MQ=255)
gCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGc                                                    <  1:1965212/69‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:2003088/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:3249401/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:1865518/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:3079098/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:3022970/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:249682/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:2004928/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:2014183/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:2287219/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:239285/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:2412138/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:2522646/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:1680450/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:1669163/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:664510/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:1135219/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:114562/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:1159318/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:585186/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:1179104/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:1202833/67‑1 (MQ=255)
     gATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTCGTGGGCACGGAACTTAACGCCac                                                 <  1:831766/67‑1 (MQ=255)
      aTCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCac                                                 <  1:1700330/66‑1 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:661855/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:2968798/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:2895240/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:2840201/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:471473/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:461854/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:3401080/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:2636496/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:3217247/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:2540571/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:697290/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:3382143/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:999389/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:2408716/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:1956140/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:1833745/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:2373356/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:1972618/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:2001050/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:2414106/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:1108790/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacaca                                              >  1:1402492/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTGGGCACGGAACTTAACGCCccaca                                              >  1:1016997/1‑68 (MQ=255)
       tCGATGGCCGGATTTATGATGACCCGCGAAGGACAGTGTGTGGTGGGCACGGAACTTAACGCcacacat                                             >  1:2561162/1‑68 (MQ=255)
                 gATTTATGCTGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACCCCATCGCCCg                                    >  1:3333874/1‑68 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTTCgcggt            <  1:436672/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:1504631/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:3316989/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:448566/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:372544/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:461375/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:1492707/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:627524/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:659135/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:1333072/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:1099540/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:1086490/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:927374/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:930038/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:1568018/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:1722071/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:3237502/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:280819/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:3141555/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:3090875/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:2243212/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:3034591/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:2064248/69‑2 (MQ=255)
                                        cAGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:2190099/69‑2 (MQ=255)
                                        cAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACgcggc            >  1:805236/1‑69 (MQ=255)
                                        cAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTATCTGAGGGAAAGGTACgcggc            >  1:3068177/1‑69 (MQ=255)
                                        cAGTGTGTGGGGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:1399713/69‑2 (MQ=255)
                                        cAGTGTGTGGGGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:1930869/69‑2 (MQ=255)
                                         aGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:2326791/68‑2 (MQ=255)
                                         aGTGTGTGGTGGGCACGGAACTTAACGCCACACATCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggt            <  1:2608101/68‑2 (MQ=255)
                                                  tAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAg   >  1:1702730/1‑68 (MQ=255)
                                                  tAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGc  >  1:2536953/1‑69 (MQ=255)
                                                   |                                                                   
GCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGC  >  minE/384388‑384506

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: