Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,515,674 A→G 7.8% A114A (GCT→GCC cvpA ← membrane protein required for colicin V production

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,515,6740AG 7.8% 302.5 / 11.9 102A114A (GCT→GCCcvpAmembrane protein required for colicin V production
Reads supporting (aligned to +/- strand):  ref base A (45/49);  new base G (2/6);  total (47/55)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 103 total aligned reads displayed.

GTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGCCTGACAACCC  >  minE/1515610‑1515742
                                                                |                                                                    
gtgATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGAAGcaac                                                                 >  1:4181876/1‑70 (MQ=255)
 tgATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGcaac                                                                 <  1:2099471/69‑1 (MQ=255)
   aTTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAAt                                                              >  1:3358004/1‑70 (MQ=255)
   aTTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAAt                                                              >  1:1171598/1‑70 (MQ=255)
     ttGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCGGCAACAATCaa                                                           >  1:2053290/1‑71 (MQ=255)
     ttGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCGGCAACAATCaa                                                           >  1:4165437/1‑71 (MQ=255)
      tGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCaa                                                           <  1:2327828/70‑1 (MQ=255)
      tGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCaa                                                           >  1:232326/1‑70 (MQ=255)
      tGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCaa                                                           <  1:2029814/70‑1 (MQ=255)
      tGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCaa                                                           >  1:866102/1‑70 (MQ=255)
       gCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCaa                                                           <  1:1102742/69‑1 (MQ=255)
          ccAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAacacac                                                      >  1:4396535/1‑71 (MQ=255)
          ccAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAacacac                                                      >  1:4780725/1‑71 (MQ=255)
           cAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACAcc                                                     >  1:246270/1‑71 (MQ=255)
           cAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACAcc                                                     >  1:177600/1‑71 (MQ=255)
            aGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACcg                                                    >  1:157232/1‑71 (MQ=255)
                ttCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACAcc                                                     <  1:4338244/66‑1 (MQ=255)
                ttCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCgcgc                                                 >  1:3324452/1‑70 (MQ=255)
                ttCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCgcgc                                                 >  1:2921515/1‑70 (MQ=255)
                ttCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCgcgc                                                 >  1:1959088/1‑70 (MQ=255)
                ttCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCgcgc                                                 >  1:5195688/1‑70 (MQ=255)
                     ttttCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCGGCAACAATCAACACACcgc                                                   <  1:1965355/63‑1 (MQ=255)
                     ttttCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCGGCAACAATCAACACACcgc                                                   <  1:342046/63‑1 (MQ=255)
                     ttttCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCGGCAACAATCAACACACcgc                                                   <  1:3684561/63‑1 (MQ=255)
                     ttttCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCGGCAACAATCAACACACcgc                                                   <  1:684391/63‑1 (MQ=255)
                     ttttCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCGGCAACAATCAACACACcgc                                                   <  1:2152933/63‑1 (MQ=255)
                      tttCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCGGCAACAATCAACACACcgc                                                   <  1:1487011/62‑1 (MQ=255)
                         cGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCaa                                               >  1:4343462/1‑63 (MQ=255)
                         cGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCaa                                               >  1:3410069/1‑63 (MQ=255)
                         cGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCaa                                               >  1:2489413/1‑63 (MQ=255)
                         cGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCaa                                               >  1:3968912/1‑63 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGaa                                      <  1:3946503/71‑1 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGaa                                      <  1:4940352/71‑1 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGaa                                      <  1:4401569/71‑1 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGaa                                      <  1:4351187/71‑1 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGa                                       <  1:989519/70‑1 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGa                                       <  1:4189676/70‑1 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGa                                       <  1:3016444/70‑1 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGa                                       <  1:1119980/70‑1 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGa                                       <  1:4546336/70‑1 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGa                                       <  1:1853889/70‑1 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGa                                       <  1:4825504/70‑1 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGa                                       <  1:5152027/70‑1 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGa                                       <  1:1179703/70‑1 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGa                                       <  1:1274628/70‑1 (MQ=255)
                          gACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCCCGCAACGCACCGa                                       <  1:1857907/70‑1 (MQ=255)
                           acacCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGTCAGCAACAATCAACACACCGCGCAACGCACCGa                                       <  1:829293/69‑1 (MQ=255)
                                     aaaGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGcc                             >  1:2799266/1‑69 (MQ=255)
                                     aaaGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGcc                             >  1:2097191/1‑69 (MQ=255)
                                        ggAGTCGAGAAAGAAGAGAATGGCAGCAACAATCaa                                                           >  1:3238062/1‑36 (MQ=255)
                                        ggAGTCGAGAAAGAAGAGAATGGCAGCAACAATCaa                                                           >  1:3703158/1‑36 (MQ=255)
                                        ggAGTCGAGAAAGAAGAGAATGGCAGCAACAATCaa                                                           >  1:3899929/1‑36 (MQ=255)
                                        ggAGTCGAGAAAGAAGAGAATGGCAGCAACAATCaa                                                           >  1:4938817/1‑36 (MQ=255)
                                        ggAGTCGAGAAAGAAGAGAATGGCAGCAACAATCaa                                                           >  1:4193534/1‑36 (MQ=255)
                                        ggAGTCGAGAAAGAAGAGAATGGCAGCAACAATCaa                                                           >  1:4463264/1‑36 (MQ=255)
                                        ggAGTCGAGAAAGAAGAGAATGGCAGCAACAATCaa                                                           >  1:140445/1‑36 (MQ=255)
                                        ggAGTCGAGAAAGAAGAGAATGGCAGCAACAATCaa                                                           >  1:319032/1‑36 (MQ=255)
                                        ggAGTCGAGAAAGAAGAGAATGGAAGCAACAATCaa                                                           >  1:1004200/1‑36 (MQ=255)
                                            tCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCgcgc                                                 <  1:3404376/42‑1 (MQ=255)
                                            tCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCgcgc                                                 <  1:106834/42‑1 (MQ=255)
                                            tCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCgcgc                                                 <  1:1411583/42‑1 (MQ=255)
                                            tCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCgcgc                                                 <  1:3217312/42‑1 (MQ=255)
                                            tCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCg                    <  1:4498374/71‑1 (MQ=255)
                                              gagaAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGAt                  >  1:792740/1‑71 (MQ=255)
                                              gagaAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGAt                  >  1:1497034/1‑71 (MQ=255)
                                              gagaAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGAt                  >  1:1592797/1‑71 (MQ=255)
                                               agaaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACa                                   <  1:675245/53‑1 (MQ=255)
                                               agaaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACa                                   <  1:1923791/53‑1 (MQ=255)
                                               agaaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATc                 <  1:3583042/71‑1 (MQ=255)
                                               agaaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATc                 <  1:300430/71‑1 (MQ=255)
                                               agaaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATc                 <  1:2681782/71‑1 (MQ=255)
                                               agaaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACCCCCAGCACCCGATc                 <  1:3893574/71‑1 (MQ=255)
                                                gaaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATc                 <  1:3560685/70‑1 (MQ=255)
                                                 aaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGa                                 >  1:931868/1‑53 (MQ=255)
                                                 aaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGa                                 >  1:200527/1‑53 (MQ=255)
                                                 aaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGa                                 >  1:3354156/1‑53 (MQ=255)
                                                 aaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGa                                 >  1:3021717/1‑53 (MQ=255)
                                                 aaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGa                                 >  1:2748435/1‑53 (MQ=255)
                                                 aaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAg                          <  1:1152638/60‑1 (MQ=255)
                                                 aaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAg                          <  1:948200/60‑1 (MQ=255)
                                                 aaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAg                          <  1:4825943/60‑1 (MQ=255)
                                                 aaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAg                          <  1:262280/60‑1 (MQ=255)
                                                 aaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCgg               >  1:4454878/1‑71 (MQ=255)
                                                 aaagaaGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCgg               >  1:4448652/1‑71 (MQ=255)
                                                  aagaagAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAg                          <  1:2996800/59‑1 (MQ=255)
                                                    gaagaGAATTGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGc            <  1:205427/71‑1 (MQ=255)
                                                    gaagaGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGc            <  1:5023072/71‑1 (MQ=255)
                                                    gaagaGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGc            <  1:3008375/71‑1 (MQ=255)
                                                     aagaGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGc            <  1:4345039/70‑1 (MQ=255)
                                                      agagAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGCCt          <  1:780759/71‑1 (MQ=255)
                                                       gagaATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGCCTg         <  1:3171966/71‑1 (MQ=255)
                                                         gaATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGCCt          <  1:4253633/68‑1 (MQ=255)
                                                         gaATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGCCt          <  1:1076415/68‑1 (MQ=255)
                                                         gaATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGCCt          <  1:1554338/68‑1 (MQ=255)
                                                         gaATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGCCTGa        >  1:1641915/1‑70 (MQ=255)
                                                         gaATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGCCTGa        >  1:2579215/1‑70 (MQ=255)
                                                         gaATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGCCTGAc       <  1:3238440/71‑1 (MQ=255)
                                                         gaATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGCCTGAc       <  1:3804840/71‑1 (MQ=255)
                                                         gaAAGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACa                                   >  1:2477616/1‑43 (MQ=255)
                                                               cagcaACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGCCTGACAAccc  >  1:314271/1‑70 (MQ=255)
                                                                |                                                                    
GTGATTTGCTCCAGTCTTCGCTTTTCGACACCCCGGTAAAGGAGTCGAGAAAGAAGAGAATGGCAGCAACAATCAACACACCGCGCAACGCACCGAAACAGACGCCCAGCACCCGATCGGTGCCTGACAACCC  >  minE/1515610‑1515742

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: