Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,473 C→T 7.6% A277T (GCA→ACA) ‡ ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,4730CT 7.6% 332.6 / 11.2 105A277T (GCA→ACA) ‡ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base C (56/41);  new base T (3/5);  total (59/46)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 110 total aligned reads displayed.

CAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCAGGTAG  >  minE/1439405‑1439543
                                                                    |                                                                      
cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt                                                                      >  1:1470464/1‑71 (MQ=255)
cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt                                                                      >  1:434528/1‑71 (MQ=255)
cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt                                                                      >  1:465393/1‑71 (MQ=255)
cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt                                                                      >  1:3622126/1‑71 (MQ=255)
  ggtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTg                                                                    <  1:1808799‑M2/66‑1 (MQ=255)
  ggtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTg                                                                    >  1:2432084‑M2/6‑71 (MQ=255)
  ggTAAGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTg                                                                    <  1:671175/71‑1 (MQ=12)
   gtaaGCCATGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGt                                                                   >  1:2553311‑M2/5‑71 (MQ=255)
   gtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGt                                                                   >  1:3472081‑M2/5‑71 (MQ=255)
   gtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGt                                                                   >  1:1647389‑M2/5‑71 (MQ=255)
    taaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGtt                                                                  <  1:3317240‑M2/68‑1 (MQ=255)
     aaGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGt                                                                   >  1:3566872/1‑69 (MQ=12)
     aaGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGt                                                                   >  1:3777063/1‑69 (MQ=12)
     aaGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGt                                                                   >  1:1591140/1‑69 (MQ=12)
      aGCCAGGGCCGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc                                                                 <  1:1460458‑M2/69‑1 (MQ=255)
      aGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc                                                                 <  1:3164980‑M2/69‑1 (MQ=255)
          aGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGccc                                                             >  1:2357432/1‑70 (MQ=255)
          aGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGc                                                               >  1:1691457/1‑68 (MQ=255)
           gggACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCaa                                                           >  1:2678961/1‑71 (MQ=255)
              acggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAAccc                                                        <  1:3771531/71‑1 (MQ=255)
               cggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAAccc                                                        <  1:79359/70‑1 (MQ=255)
                  acgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg                                                     >  1:1629976/1‑70 (MQ=255)
                  acgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg                                                     >  1:3512753/1‑70 (MQ=255)
                  acgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg                                                     >  1:3947181/1‑70 (MQ=255)
                    gCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGcc                                                              >  1:2047331/1‑59 (MQ=255)
                    gCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGaa                                                  >  1:3477011/1‑71 (MQ=255)
                     cAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAAc                                                 <  1:2047919/71‑1 (MQ=255)
                     cAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAAc                                                 <  1:1950751/71‑1 (MQ=255)
                      aGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc                                                                 <  1:3145725/54‑1 (MQ=255)
                      aGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAAcc                                                >  1:2904924/1‑71 (MQ=255)
                       gACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCt                                               <  1:1629197/71‑1 (MQ=255)
                        aCCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTa                                              >  1:3220294/1‑71 (MQ=255)
                        aCCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTa                                              >  1:2061859/1‑71 (MQ=255)
                        aCCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTATTCGCCCAACCCAGGGAACCTa                                              >  1:1071219/1‑71 (MQ=255)
                         ccGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCg                                                                >  1:2651472/1‑52 (MQ=255)
                         ccGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGaa                                                  >  1:1765266/1‑66 (MQ=255)
                          cGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc                                                                 <  1:589999/50‑1 (MQ=255)
                          cGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTa                                              <  1:3089718/69‑1 (MQ=255)
                           gAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTAcgc                                           <  1:1510194/71‑1 (MQ=255)
                           aaaGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGtttt                                                                     <  1:898203/44‑1 (MQ=16)
                            aaGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCttt                                                                     >  1:1918927/1‑44 (MQ=255)
                            aaGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCttt                                                                     >  1:3850996/1‑44 (MQ=255)
                            aaGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGc                                                               >  1:3914051/1‑50 (MQ=255)
                              gTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAg                                        <  1:3332638/71‑1 (MQ=255)
                               tCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTg                                                                    <  1:3077257/42‑1 (MQ=255)
                                cGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGtt                                      >  1:1787420/1‑71 (MQ=255)
                                 gAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGTTAGcga                                                             <  1:3387140/47‑3 (MQ=12)
                                 gAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGTTAGcg                                                              <  1:1436413/46‑2 (MQ=12)
                                 gAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGTTAGcg                                                              <  1:503352/46‑2 (MQ=12)
                                   aCTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTg                                                                    <  1:1482405/38‑1 (MQ=255)
                                   aCTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCg                                   >  1:338337/1‑71 (MQ=255)
                                   aCTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCg                                   >  1:2455338/1‑71 (MQ=255)
                                   aCTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCg                                   >  1:520757/1‑71 (MQ=255)
                                   aCTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCAAGGGa                                                   >  1:2513961/1‑55 (MQ=255)
                                      ggTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGccc                                                             >  1:2243854/1‑42 (MQ=255)
                                      ggTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGccc                                                             >  1:3730668/1‑42 (MQ=255)
                                      ggTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCa                                                            <  1:3922949/43‑1 (MQ=255)
                                       gTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGt                               <  1:2675591/71‑1 (MQ=255)
                                       gTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGt                               <  1:1773293/71‑1 (MQ=255)
                                         aCTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAg                             >  1:3236920/1‑71 (MQ=255)
                                          cTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTAcgcg                                          <  1:2125339/57‑1 (MQ=255)
                                          cTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAgg                            <  1:3461837/71‑1 (MQ=255)
                                          cTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAg                             <  1:646809/70‑1 (MQ=255)
                                          cTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAg                             <  1:2708610/70‑1 (MQ=255)
                                           tGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAgg                            >  1:2575439/1‑70 (MQ=255)
                                             aGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCt                         <  1:2508163/71‑1 (MQ=255)
                                             aGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCt                         <  1:1780378/71‑1 (MQ=255)
                                             aGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCt                         >  1:1667622/1‑71 (MQ=255)
                                             aGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCt                         >  1:3069065/1‑71 (MQ=255)
                                               cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTgg                       <  1:131542/71‑1 (MQ=255)
                                                aaCAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTggg                      >  1:131755/1‑71 (MQ=255)
                                                 aCAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGgt                     <  1:1860450/71‑1 (MQ=255)
                                                   aGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGgtg                    >  1:1643527/1‑70 (MQ=255)
                                                    gCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTa                  >  1:1328370/1‑71 (MQ=255)
                                                     cTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTa                  >  1:425300/1‑70 (MQ=255)
                                                      ttCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTg                                <  1:2224268/55‑1 (MQ=255)
                                                      ttCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACt                >  1:1152085/1‑71 (MQ=255)
                                                      ttCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACt                >  1:1313671/1‑71 (MQ=255)
                                                       tCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACt                >  1:2445861/1‑70 (MQ=255)
                                                       tCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACt                >  1:1131851/1‑70 (MQ=255)
                                                        cGAAGTTCTGTGCTTTGTTCGCCCAACCCAGTGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACt                <  1:1263573/69‑1 (MQ=255)
                                                        cGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACt                >  1:3770913/1‑69 (MQ=255)
                                                        cGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACt                >  1:2298744/1‑69 (MQ=255)
                                                        cGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACt                <  1:459931/69‑1 (MQ=255)
                                                        cGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTg               <  1:2606486/70‑1 (MQ=255)
                                                         gAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGag             <  1:648362/71‑1 (MQ=255)
                                                           aGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGagca           <  1:1035734/71‑1 (MQ=255)
                                                            gTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGagcag          >  1:926620/1‑71 (MQ=255)
                                                             ttCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCg                                   <  1:2503139/45‑1 (MQ=255)
                                                             ttCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCt                         <  1:321322/55‑1 (MQ=255)
                                                             ttCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGagcagc         <  1:1756390/71‑1 (MQ=255)
                                                              tCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGagcagc         >  1:75414/1‑70 (MQ=255)
                                                               cTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGcc        <  1:707176/70‑1 (MQ=255)
                                                                tgtgCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCAg      >  1:2779218/1‑71 (MQ=255)
                                                                agtgCTTTGTTCGCCCAACCCAGGGAACCTACGCGa                                         >  1:2715298/2‑36 (MQ=255)
                                                                  tgCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCAgg     <  1:3789678/70‑1 (MQ=255)
                                                                   gCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCAgg     >  1:1138251/1‑69 (MQ=255)
                                                                    cTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCAGGTa   <  1:1408470/70‑1 (MQ=255)
                                                                    cTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCAGGTAg  <  1:776883/71‑1 (MQ=255)
                                                                    cTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGAACTGAGCAGCCAgg     >  1:3311959/1‑68 (MQ=255)
                                                                    |                                                                      
CAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCAGGTAG  >  minE/1439405‑1439543

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: