Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,489,649 A→G 19.5% V216V (GTA→GTG dnaN → DNA polymerase III, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,489,6490AG19.5% 78.7 / 11.6 56V216V (GTA→GTGdnaNDNA polymerase III, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (20/25);  new base G (6/5);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GTACCGCGTAAAGGCGTGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCCAAACTGG  >  minE/2489584‑2489710
                                                                 |                                                             
gTACCGCGTAAAGGCGTGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAg                                                            <  1:2232767/69‑1 (MQ=255)
gTACCGCGTAAAGGCGTGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAg                                                            <  1:2606458/69‑1 (MQ=255)
gTACCGCGTAAAGGCGTGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAg                                                            <  1:2805602/69‑1 (MQ=255)
         aaaGGCGTGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   <  1:108596/69‑1 (MQ=255)
         aaaGGCGTGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   <  1:1994873/69‑1 (MQ=255)
              cGTGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGcaacaa                                              <  1:3153467/69‑1 (MQ=255)
              cGTGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGcaacaa                                              <  1:2828996/69‑1 (MQ=255)
              cGTGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGcaacaa                                              <  1:1067775/69‑1 (MQ=255)
                tGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGCAacaaca                                            <  1:458354/69‑1 (MQ=255)
                tGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGCAacaaca                                            <  1:1518548/69‑1 (MQ=255)
                tGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGCAacaaca                                            <  1:3388830/69‑1 (MQ=255)
                tGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGCAacaaca                                            <  1:3045169/69‑1 (MQ=255)
                tGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGCAacaaca                                            <  1:672418/69‑1 (MQ=255)
                tGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGCAacaaca                                            <  1:1233385/69‑1 (MQ=255)
                tGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGCAacaaca                                            <  1:2496430/69‑1 (MQ=255)
                 gATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGCAacaaca                                            <  1:1471497/68‑1 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACa                                                             <  1:1072963/50‑1 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACa                                                             <  1:3000410/50‑1 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACa                                                             <  1:268377/50‑1 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACa                                                             <  1:2362881/50‑1 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACa                                                             <  1:80891/50‑1 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACa                                                             <  1:1550051/50‑1 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACa                                                             <  1:2028181/50‑1 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACa                                                             <  1:1828437/50‑1 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACa                                                             <  1:1839206/50‑1 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGa                                                           >  1:1218649/1‑52 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   >  1:1620466/1‑60 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   >  1:982975/1‑60 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   >  1:515413/1‑60 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   >  1:460440/1‑60 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   >  1:1488450/1‑60 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   >  1:381760/1‑60 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   >  1:1527405/1‑60 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   >  1:2885546/1‑60 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   >  1:1899123/1‑60 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   >  1:2020122/1‑60 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   >  1:280482/1‑60 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   >  1:2071597/1‑60 (MQ=255)
                  aTTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGgcagc                                                   >  1:2339445/1‑60 (MQ=255)
                    tGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGCAACAACATTcg                                        >  1:3071738/1‑69 (MQ=255)
                                    cTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTgg                                                       >  1:2902899/1‑38 (MQ=255)
                                    cTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGAc                        >  1:1890864/1‑69 (MQ=255)
                                    cTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGAc                        >  1:1419952/1‑69 (MQ=255)
                                                      ccGCTGCGCGTACAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCCaaa      >  1:3119777/1‑69 (MQ=255)
                                                      ccGCTGCGCGTACAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCCaaa      >  1:2533158/1‑69 (MQ=255)
                                                        gCTGCGCGTGCAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCGAAAc     >  1:1910857/1‑68 (MQ=255)
                                                        gCTGCGCGTGCAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCGAAAc     >  1:1919422/1‑68 (MQ=255)
                                                        gCTGCGCGTGCAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCGAAACt    >  1:3114784/1‑69 (MQ=255)
                                                        gCTGCGCGTGCAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCGAAACt    >  1:386519/1‑69 (MQ=255)
                                                        gCTGCGCGTGCAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCGAAACt    >  1:2315356/1‑69 (MQ=255)
                                                        gCTGCGCGTGCAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCGAAACt    >  1:487520/1‑69 (MQ=255)
                                                          tGCGCGTGCAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCGAAACtgg  <  1:2057545/69‑1 (MQ=255)
                                                          tGCGCGTGCAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCGAAACtgg  <  1:3157873/69‑1 (MQ=255)
                                                          tGCGCGTGCAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCGAAACtgg  <  1:3098715/69‑1 (MQ=255)
                                                          tGCGCGTGCAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCGAAACtgg  <  1:1599602/69‑1 (MQ=255)
                                                           gcgcGTGCAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCGAAACtgg  <  1:2234021/68‑1 (MQ=255)
                                                                 |                                                             
GTACCGCGTAAAGGCGTGATTGAACTGATGCGTATGCTCGACGGCGGCGACAATCCGCTGCGCGTACAGATTGGCAGCAACAACATTCGCGCCCACGTTGGCGACTTTATCTTCACCTCCAAACTGG  >  minE/2489584‑2489710

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: