Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 139,441 C→T 5.6% S242S (TCG→TCA htrE ← predicted outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE139,4410CT 5.6% 331.6 / 6.3 106S242S (TCG→TCAhtrEpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base C (56/44);  new base T (3/3);  total (59/47)
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.60e-01

Only 100 of 112 total aligned reads displayed.

TCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGT  >  minE/139372‑139511
                                                                     |                                                                      
tCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGc                                                                        >  1:1237731/1‑70 (MQ=255)
tCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGc                                                                        >  1:586030/1‑70 (MQ=255)
 cGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGc                                                                        >  1:4009284/1‑69 (MQ=255)
  ggATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCg                                                                       <  1:356709/69‑1 (MQ=255)
  ggATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCg                                                                       <  1:1133574/69‑1 (MQ=255)
  ggATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAg                                                                     >  1:4461418/1‑71 (MQ=255)
   gATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAgg                                                                    >  1:4123446/1‑71 (MQ=255)
     tACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAgg                                                                    >  1:619415/1‑69 (MQ=255)
     tACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGc                                                                   <  1:4586373/70‑1 (MQ=255)
       cTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAg                                                                     >  1:2498632/1‑66 (MQ=255)
       cTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGa                                                                 <  1:3140726/70‑1 (MQ=255)
       cTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGa                                                                 >  1:445450/1‑70 (MQ=255)
         gACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGa                                                                 <  1:119442/68‑1 (MQ=255)
             gAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCa                                                            >  1:3166327/1‑69 (MQ=255)
             gAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACg                                                          <  1:2577496/71‑1 (MQ=255)
                tCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGc                                                                   <  1:2286835/59‑1 (MQ=255)
                 cAAAGGTTTCGCCCGTCGTATAAGTCTAACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTg                                                       >  1:3542944/1‑70 (MQ=255)
                  aaaGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGc                                                                   >  1:1538598/1‑57 (MQ=255)
                  aaaGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGaa                                                     <  1:1160348/71‑1 (MQ=255)
                  aaaGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGaa                                                     >  1:3697075/1‑71 (MQ=255)
                  aaaGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGaa                                                     <  1:2124840/71‑1 (MQ=255)
                    aGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGTGAGGCGATATCACGCTGAACa                                                   <  1:1219375/71‑1 (MQ=255)
                    aGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACa                                                   <  1:2988587/71‑1 (MQ=255)
                      gTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATa                                                 >  1:3195582/1‑71 (MQ=255)
                       tttCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATa                                                 <  1:4104631/70‑1 (MQ=255)
                       tttCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATAc                                                <  1:2744829/71‑1 (MQ=255)
                       tttCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATAc                                                <  1:4673426/71‑1 (MQ=255)
                        ttCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGat                                                                >  1:3415044/1‑54 (MQ=255)
                        ttCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATAcc                                               >  1:1650237/1‑71 (MQ=255)
                        ttCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATAcc                                               >  1:1679904/1‑71 (MQ=255)
                        ttCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATAcc                                               >  1:3351192/1‑71 (MQ=255)
                         tCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACg                                                          <  1:259991/59‑1 (MQ=255)
                         tCGCCCCTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACa                                                   >  1:604516/1‑66 (MQ=255)
                              cgtcgtATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCt                                         >  1:1788105/1‑71 (MQ=255)
                                tcgtATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCt                                         >  1:793084/1‑69 (MQ=255)
                                  gtATAAGACTCACCAAGAATGAGTTGAGAACGCAGTGAGGCGa                                                                 >  1:4018680/1‑43 (MQ=255)
                                   tataAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCAc                                                           <  1:1178595/48‑1 (MQ=255)
                                    ataAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGaa                                                     <  1:181433/53‑1 (MQ=255)
                                     taAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTg                                                       <  1:2017695/50‑1 (MQ=255)
                                     taAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAAt                                   >  1:2748206/1‑70 (MQ=255)
                                     taAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAAt                                   >  1:160002/1‑70 (MQ=255)
                                     taAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAAt                                   <  1:3714228/70‑1 (MQ=255)
                                       aGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGaa                                                     >  1:264332/1‑50 (MQ=255)
                                       aGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATa                                                 >  1:2957269/1‑54 (MQ=255)
                                       aGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATAcc                                               >  1:1328563/1‑56 (MQ=255)
                                       aGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCa                                 >  1:3360763/1‑70 (MQ=255)
                                       aGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCa                                 >  1:4007521/1‑70 (MQ=255)
                                       aGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCa                                 >  1:3978921/1‑70 (MQ=255)
                                       aGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCa                                 >  1:1448117/1‑70 (MQ=255)
                                       aGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCa                                 >  1:3573198/1‑70 (MQ=255)
                                          cTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGaa                                                     >  1:3730581/1‑47 (MQ=255)
                                          cTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATc                                  <  1:1690020/66‑1 (MQ=255)
                                          cTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAAt                             <  1:1111016/71‑1 (MQ=255)
                                           tCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCt                                                        <  1:3487763/43‑1 (MQ=255)
                                           tCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCAt                                <  1:3879262/67‑1 (MQ=255)
                                             aCCAAGAATGAGTTGAGAACGCAGTGAGGCGATATCa                                                            >  1:437793/1‑37 (MQ=255)
                                             aCCAAGAATGAGTTGAGAACGCAGTGAGGCGATATCa                                                            >  1:2161301/1‑37 (MQ=255)
                                             aCCAAGAATGAGTTGAGAACGCAGTGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTAc                          <  1:2324815/71‑1 (MQ=255)
                                             aCCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATaa                              >  1:717911/1‑67 (MQ=255)
                                             aCCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATa                               <  1:3057123/66‑1 (MQ=255)
                                              ccAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACt                         >  1:2477554/1‑71 (MQ=255)
                                               cAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATa                               <  1:4515200/64‑1 (MQ=255)
                                                aaGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATaa                              >  1:1602160/1‑64 (MQ=255)
                                                aaGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGc                       >  1:4406402/1‑71 (MQ=255)
                                                aaGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGc                       >  1:1539883/1‑71 (MQ=255)
                                                aaGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGc                       >  1:87516/1‑71 (MQ=255)
                                                 aGAATGAGTTGAGAACGCAGCGAGGCGATATCCCGCTGAACATACCGATTCTTaaaa                                    >  1:2784219/1‑57 (MQ=255)
                                                 aGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGaa                                                     <  1:3074589/40‑1 (MQ=255)
                                                 aGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATaa                              >  1:4294300/1‑63 (MQ=255)
                                                 aGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATaa                              >  1:2883953/1‑63 (MQ=255)
                                                   aaTGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGc                       >  1:4263819/1‑68 (MQ=255)
                                                    aTGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATc                                  <  1:1777074/56‑1 (MQ=255)
                                                    aTGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAg                    <  1:3285887/70‑1 (MQ=255)
                                                    aTGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAg                    <  1:1968969/70‑1 (MQ=255)
                                                    aTGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGa                   <  1:2444346/71‑1 (MQ=255)
                                                     tGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCg                                              >  1:3109231/1‑43 (MQ=255)
                                                        gTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCg               <  1:2080839/71‑1 (MQ=255)
                                                        gTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCg               <  1:2460008/71‑1 (MQ=255)
                                                         ttGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATc                                  <  1:1029496/51‑1 (MQ=255)
                                                         ttGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCgg              <  1:1844680/71‑1 (MQ=255)
                                                           gagaACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATaa                              >  1:2459663/1‑53 (MQ=255)
                                                           gagaACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATaa                              >  1:3478955/1‑53 (MQ=255)
                                                           gagaACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGt             >  1:2578394/1‑70 (MQ=255)
                                                           gagaACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGt             >  1:2738068/1‑70 (MQ=255)
                                                           gagaACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGt             >  1:2743490/1‑70 (MQ=255)
                                                           gagaACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGt             <  1:1819094/70‑1 (MQ=255)
                                                            agaACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCa           >  1:3748594/1‑71 (MQ=255)
                                                             gaACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCAt          >  1:1088009/1‑71 (MQ=255)
                                                              aaCGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAg                    <  1:1094865/60‑1 (MQ=255)
                                                              aaCGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAg                    <  1:1010179/60‑1 (MQ=255)
                                                                cGCAGCGAGGCGATATCACGCTGAACATACCGATTCt                                         >  1:2935916/1‑37 (MQ=255)
                                                                cGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAAt                             >  1:2734994/1‑49 (MQ=255)
                                                                cGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAAt                             >  1:3548156/1‑49 (MQ=255)
                                                                 gCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTa                           >  1:2885949/1‑50 (MQ=255)
                                                                   aGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGtt    >  1:2670022/1‑71 (MQ=255)
                                                                    gCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAAt                             <  1:1856230/45‑1 (MQ=255)
                                                                    gCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCTGTCATCCAGt     >  1:4529828/1‑69 (MQ=255)
                                                                     cGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGt  <  1:4499533/71‑1 (MQ=255)
                                                                     cGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGt  <  1:632906/71‑1 (MQ=255)
                                                                     cGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGt  <  1:2270619/71‑1 (MQ=255)
                                                                     |                                                                      
TCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGT  >  minE/139372‑139511

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: