Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 138,390 A→G 7.9% S593P (TCA→CCA)  htrE ← predicted outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE138,3900AG 7.9% 213.5 / 7.3 76S593P (TCA→CCA) htrEpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base A (34/36);  new base G (3/3);  total (37/39)
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.98e-01

CTGCTAACGTTGAGTTGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTG  >  minE/138320‑138460
                                                                      |                                                                      
ctgctAACGTTGAGTTGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGa                                                                        <  1:1006760/71‑1 (MQ=255)
   ctAACGTTGAGTTGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGGAg                                                                      >  1:1302196/1‑70 (MQ=255)
   ctAACGTTGAGTTGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAg                                                                      <  1:1386946/70‑1 (MQ=255)
    tAACGTTGAGTTGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAg                                                                      <  1:3002416/69‑1 (MQ=255)
      aCGTTGAGTTGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACg                                                                  <  1:1064727/71‑1 (MQ=255)
         ttGAGTTGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACg                                                                  <  1:3333057/68‑1 (MQ=255)
          tGAGTTGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACgttgt                                                              >  1:615209/1‑71 (MQ=255)
            aGTTGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACgttgtt                                                             <  1:340429/70‑1 (MQ=255)
            aGTTGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTc                                                            <  1:2568030/71‑1 (MQ=255)
             gTTGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACgttgtt                                                             >  1:3330556/1‑69 (MQ=255)
             gTTGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCa                                                           <  1:3764596/71‑1 (MQ=255)
              ttGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGt                                                                     <  1:1752647/60‑1 (MQ=255)
              ttGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAg                                                          >  1:2500014/1‑71 (MQ=255)
                ggTTATTACCCTAAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGATCAGTg                                                        >  1:3301964/1‑71 (MQ=255)
                  ttattaCCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACg                                                                  <  1:3465801/59‑1 (MQ=255)
                  ttattaCCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACg                                                                  <  1:3227069/59‑1 (MQ=255)
                  ttattaCCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACg                                                                  <  1:3784008/59‑1 (MQ=255)
                       aCCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAg                                                   >  1:2994735/1‑69 (MQ=255)
                        cccTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGt                                                  <  1:3820370/69‑1 (MQ=255)
                         ccTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAAt                                               <  1:193147/71‑1 (MQ=255)
                         ccTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAAt                                               <  1:1192474/71‑1 (MQ=255)
                           ttAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACgttgtt                                                             <  1:3890438/55‑1 (MQ=255)
                           ttAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAAt                                               >  1:943227/1‑69 (MQ=255)
                            tAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAAtttt                                            >  1:1621286/1‑71 (MQ=255)
                             aaaGTCACTGCTTATTTGAGTATGAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAAttttt                                           >  1:2547254/1‑71 (MQ=255)
                             aaaGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGaa                                                                       >  1:1265277/1‑43 (MQ=255)
                             aaaGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGGAGTACGTTGTTCAGTGCCAAGTAAttt                                             <  1:2578607/69‑1 (MQ=255)
                             aaaGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGGAGTACGTTGTTCAGTGCCAAGTAAttt                                             <  1:514565/69‑1 (MQ=255)
                             aaaGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGGAGTACGTTGTTCAGTGCCAAGTAAttt                                             <  1:1447476/69‑1 (MQ=255)
                             aaaGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAAttttt                                           >  1:3806254/1‑71 (MQ=255)
                             aaaGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAAttttt                                           >  1:3504071/1‑71 (MQ=255)
                              aaGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCa                                                           >  1:3440942/1‑54 (MQ=255)
                               aGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCa                                         <  1:992653/71‑1 (MQ=255)
                               aGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCa                                         <  1:3674709/71‑1 (MQ=255)
                                 tCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCa                                                           <  1:3571074/51‑1 (MQ=255)
                                 tCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCa                                                           <  1:2898386/51‑1 (MQ=255)
                                  cACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCa                                         >  1:642469/1‑68 (MQ=255)
                                   aCTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAAt                                               >  1:779731/1‑61 (MQ=255)
                                    cTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAAt                                               <  1:3125554/60‑1 (MQ=255)
                                    cTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCa                                         >  1:1708467/1‑66 (MQ=255)
                                     tGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAg                                                          >  1:300446/1‑48 (MQ=255)
                                      gCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTc                                          <  1:2594728/63‑1 (MQ=255)
                                       cTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGaa                                  <  1:2887390/70‑1 (MQ=255)
                                       cTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAAt                                 >  1:2756936/1‑71 (MQ=255)
                                       cTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAAt                                 <  1:3740775/71‑1 (MQ=255)
                                        ttATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAAt                                       <  1:3267211/64‑1 (MQ=255)
                                        ttATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATg                                >  1:2175567/1‑71 (MQ=255)
                                          aTTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATg                                <  1:1127118/69‑1 (MQ=255)
                                          aTTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATg                                <  1:1550251/69‑1 (MQ=255)
                                           tttGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCa                                                           <  1:3530143/41‑1 (MQ=255)
                                           tttGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGt                              <  1:774434/70‑1 (MQ=255)
                                             tGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGt                                                         >  1:870520/1‑41 (MQ=255)
                                             tGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGt                                                         >  1:20871/1‑41 (MQ=255)
                                             tGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGa                            >  1:2384632/1‑70 (MQ=255)
                                                gTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAAc                         >  1:1995350/1‑70 (MQ=255)
                                                gTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACt                        <  1:1645564/71‑1 (MQ=255)
                                                gTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACt                        <  1:3855869/71‑1 (MQ=255)
                                                gTATCAATACTCTGGAAACCCGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAAc                         >  1:1987171/1‑70 (MQ=255)
                                                  aTCAATACTCTGGAAACCTGAAGTACGTTGTTCAg                                                          <  1:1799176/35‑1 (MQ=255)
                                                    cAATACTCTGGAAACCTGAAGTACGTTGTTCAGTg                                                        >  1:589498/1‑35 (MQ=255)
                                                     aaTACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGa                    >  1:2080441/1‑70 (MQ=255)
                                                      aTACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCaa                                                    <  1:1481567/37‑1 (MQ=255)
                                                      aTACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCaa                                                    <  1:313916/37‑1 (MQ=255)
                                                       tACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATa                  >  1:3491513/1‑70 (MQ=255)
                                                        aCTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGt                                                  >  1:1780574/1‑37 (MQ=255)
                                                        aCTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGt                                                  >  1:1607216/1‑37 (MQ=255)
                                                        aCTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGa                                   >  1:2169555/1‑52 (MQ=255)
                                                            tGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTc                                          <  1:2099949/41‑1 (MQ=255)
                                                            tGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGc             >  1:460578/1‑70 (MQ=255)
                                                                aaCCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGaaa                          >  1:2395581/1‑53 (MQ=255)
                                                                aaCCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGaaa                          >  1:375659/1‑53 (MQ=255)
                                                                   cTGGAGTACGTTGTTCAGTGCCAAGTAATTTTTCaa                                        >  1:564797/1‑36 (MQ=255)
                                                                   cTGGAGTACGTTGTTCAGTGCCAAGTAATTTTTCaa                                        >  1:1847238/1‑36 (MQ=255)
                                                                   cTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCg        >  1:1356001/1‑68 (MQ=255)
                                                                   cTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGt       >  1:1048384/1‑69 (MQ=255)
                                                                   cTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCa     <  1:948056/71‑1 (MQ=255)
                                                                    tGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTgg                                    <  1:165582/39‑1 (MQ=255)
                                                                     gAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAg    >  1:2296904/1‑70 (MQ=255)
                                                                      aaGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAgtg  >  1:2427264/1‑71 (MQ=255)
                                                                      |                                                                      
CTGCTAACGTTGAGTTGGTTATTACCCTTAAAGTCACTGCTTATTTGAGTATCAATACTCTGGAAACCTGAAGTACGTTGTTCAGTGCCAAGTAATTTTTCAATTGGAATGGTGAAACTAAGATAAACGCTATCGTCAGTG  >  minE/138320‑138460

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: