Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,545,256 C→T 5.9% A81T (GCG→ACG)  rfaC ← ADP‑heptose:LPS heptosyl transferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,545,2560CT 5.9% 374.2 / 9.1 118A81T (GCG→ACG) rfaCADP‑heptose:LPS heptosyl transferase I
Reads supporting (aligned to +/- strand):  ref base C (33/78);  new base T (2/5);  total (35/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 123 total aligned reads displayed.

CAGCGCCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGT  >  Exported/2545191‑2545320
                                                                 |                                                                
cAGCGCCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCa                                                              >  1:4785181/1‑70 (MQ=255)
cAGCGCCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCAc                                                             >  1:5093466/1‑71 (MQ=255)
cAGCGCCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCAc                                                             >  1:2947761/1‑71 (MQ=255)
cAGCGCCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCAc                                                             <  1:3038880/71‑1 (MQ=255)
  gcgcCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACg                                                            >  1:4525033/1‑70 (MQ=255)
  gcgcCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACg                                                            >  1:455609/1‑70 (MQ=255)
  gcgcCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGa                                                           <  1:102965/71‑1 (MQ=255)
  gcgcCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGa                                                           <  1:2600203/71‑1 (MQ=255)
  gcgcCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGa                                                           <  1:3925577/71‑1 (MQ=255)
  gcgcCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGa                                                           <  1:1774763/71‑1 (MQ=255)
  gcgcCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGa                                                           <  1:1468892/71‑1 (MQ=255)
    gccgcGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCa                                                              >  1:4262616/1‑66 (MQ=255)
     ccgcGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAAc                                                        <  1:3804144/71‑1 (MQ=255)
     ccgcGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAAc                                                        <  1:4344104/71‑1 (MQ=255)
     ccgcGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAAc                                                        <  1:1108369/71‑1 (MQ=255)
     ccgcGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAAc                                                        <  1:3979380/71‑1 (MQ=255)
     ccgcGCTTTTTACAAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAAc                                                        >  1:5088682/1‑71 (MQ=255)
      cgcgCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACg                                                       <  1:2466778/71‑1 (MQ=255)
          cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                   <  1:2783851/71‑1 (MQ=255)
          cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                   <  1:2125700/71‑1 (MQ=255)
          cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                   <  1:3570713/71‑1 (MQ=255)
          cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                   <  1:2036999/71‑1 (MQ=255)
          cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                   <  1:440054/71‑1 (MQ=255)
          cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                   <  1:1198068/71‑1 (MQ=255)
          cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                   <  1:1150426/71‑1 (MQ=255)
           tttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGTTTCGTGAAGCGCt                                                     <  1:1922041/68‑1 (MQ=255)
           tttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGTTTCGTGAAGCGCt                                                     <  1:4336245/68‑1 (MQ=255)
           tttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGTTTCGTGAAGCGCt                                                     <  1:1000565/68‑1 (MQ=255)
           tttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGTTTCGTGAAGCGCt                                                     <  1:4965671/68‑1 (MQ=255)
            ttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGa                                                           <  1:4085135/61‑1 (MQ=255)
            ttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGa                                                           <  1:5110716/61‑1 (MQ=255)
            ttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGa                                                           <  1:2791315/61‑1 (MQ=255)
            ttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGc                                                 >  1:4245418/1‑71 (MQ=255)
            ttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGc                                                 >  1:4708153/1‑71 (MQ=255)
            ttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGc                                                 >  1:3044812/1‑71 (MQ=255)
            ttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGc                                                 >  1:3616911/1‑71 (MQ=255)
             tttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGTTTCGTGAAGCGCtac                                                   >  1:3412417/1‑66 (MQ=37)
             tttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGTTTCGTGAAGCGCtac                                                   >  1:1482223/1‑66 (MQ=37)
             tttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCg                                                >  1:3973205/1‑71 (MQ=255)
                 ccAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTc                                             <  1:4963426/70‑1 (MQ=255)
                 ccAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTc                                             <  1:550479/70‑1 (MQ=255)
                 ccAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTc                                             <  1:3365253/70‑1 (MQ=255)
                 ccAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTc                                             <  1:521074/70‑1 (MQ=255)
                    gCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCTCTTCACGa                                                           <  1:1094338/53‑1 (MQ=255)
                       ccTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCtttt                                      <  1:4574140/71‑1 (MQ=255)
                       ccTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCtttt                                      <  1:3628929/71‑1 (MQ=255)
                         tGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                   <  1:4081037/56‑1 (MQ=255)
                         tGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                   <  1:2176097/56‑1 (MQ=255)
                         tGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                   <  1:4771655/56‑1 (MQ=255)
                         tGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                   <  1:4996978/56‑1 (MQ=255)
                          gAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATg                                   >  1:2790271/1‑71 (MQ=255)
                          gAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATg                                   >  1:3432281/1‑71 (MQ=255)
                           aGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTAt                                    <  1:3329540/69‑1 (MQ=255)
                           aGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTAt                                    <  1:2647311/69‑1 (MQ=255)
                           aGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATg                                   <  1:542549/70‑1 (MQ=255)
                                 gATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGcc                             >  1:2596579/1‑70 (MQ=255)
                                   tAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCgag                          <  1:2421958/71‑1 (MQ=255)
                                   tAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCgag                          <  1:3017815/71‑1 (MQ=255)
                                   tAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCgag                          <  1:4275101/71‑1 (MQ=255)
                                    aaCTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTAt                                    >  1:3659309/1‑60 (MQ=255)
                                    aaCTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATg                                   >  1:879637/1‑61 (MQ=255)
                                    aaCTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATg                                   >  1:1747926/1‑61 (MQ=255)
                                    aaCTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTAGCGTTCCGCTTTTAt                                    >  1:4556396/1‑60 (MQ=255)
                                     aCTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGTTTTTATGGGGGCCGAGaa                        <  1:3869313/71‑1 (MQ=255)
                                     aCTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGaa                        >  1:4544385/1‑71 (MQ=255)
                                           tCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCttt                                       <  1:2172100/50‑1 (MQ=255)
                                           tCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCttt                                       <  1:3859431/50‑1 (MQ=255)
                                             aTAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCtt                                        >  1:1462039/1‑47 (MQ=255)
                                             aTAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCtt                                        >  1:922801/1‑47 (MQ=255)
                                               aGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCttt               >  1:3613310/1‑70 (MQ=255)
                                               aGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa              >  1:1843934/1‑71 (MQ=255)
                                                 ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCttt                                       <  1:2736096/44‑1 (MQ=255)
                                                 ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCttt                                       <  1:1540566/44‑1 (MQ=255)
                                                 ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCttt                                       <  1:4622100/44‑1 (MQ=255)
                                                 ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTa                                     <  1:1630178/46‑1 (MQ=255)
                                                 ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTa                                     <  1:2884280/46‑1 (MQ=255)
                                                 ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTa                                     <  1:2546798/46‑1 (MQ=255)
                                                 ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTa                                     <  1:2368846/46‑1 (MQ=255)
                                                 ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa              <  1:4216265/69‑1 (MQ=255)
                                                 ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa              <  1:1980312/69‑1 (MQ=255)
                                                 ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa              <  1:129987/69‑1 (MQ=255)
                                                 ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa              <  1:2244955/69‑1 (MQ=255)
                                                 ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa              <  1:4158796/69‑1 (MQ=255)
                                                   ctctGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGc           >  1:5000827/1‑70 (MQ=255)
                                                   ctctGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGc           >  1:1103572/1‑70 (MQ=255)
                                                       gCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGc       >  1:863079/1‑70 (MQ=255)
                                                        cTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCg      <  1:4155942/70‑1 (MQ=255)
                                                        cTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCg      <  1:644275/70‑1 (MQ=255)
                                                        cTTGTAGCGCTTCACGAAACGATTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCg      <  1:4631338/70‑1 (MQ=255)
                                                           gTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACg   <  1:889786/70‑1 (MQ=255)
                                                           gTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACg   >  1:2249930/1‑70 (MQ=255)
                                                            tAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGt  <  1:1980800/70‑1 (MQ=255)
                                                            tAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGt  <  1:4423790/70‑1 (MQ=255)
                                                            tAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGt  <  1:1387420/70‑1 (MQ=255)
                                                             aGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGt  <  1:1412717/69‑1 (MQ=255)
                                                               cgcTTCACGAAACGCTTTGCGTTCCTCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCTGCGACg   <  1:4195220/66‑1 (MQ=255)
                                                                 cTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCtt                <  1:2479647/51‑1 (MQ=255)
                                                                 cTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCtt                <  1:4659195/51‑1 (MQ=255)
                                                                 cTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCtt                <  1:4403758/51‑1 (MQ=255)
                                                                 cTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGCCCGAGAACCAGGCtt                <  1:4432958/51‑1 (MQ=255)
                                                                 |                                                                
CAGCGCCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGT  >  Exported/2545191‑2545320

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: