Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,807,377 T→C 25.2% intergenic (+58/+6) emrB → / ← luxS multidrug efflux system protein/S‑ribosylhomocysteinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,807,3770TC25.2% 177.2 / 67.2 99intergenic (+58/+6)emrB/luxSmultidrug efflux system protein/S‑ribosylhomocysteinase
Reads supporting (aligned to +/- strand):  ref base T (32/42);  new base C (7/18);  total (39/60)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

TGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCG  >  minE/1807310‑1807445
                                                                   |                                                                    
tGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACt                                                                   <  1:302199/71‑1 (MQ=255)
 gcgcACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACt                                                                   >  1:3311785/1‑70 (MQ=255)
     aCTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTa                                                              <  1:4074257/71‑1 (MQ=255)
      cTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGAc                                                                >  1:3372775/1‑68 (MQ=255)
      cTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGAc                                                                >  1:3707541/1‑68 (MQ=255)
        aaGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTa                                                                     >  1:4090686/1‑61 (MQ=255)
        aaGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTa                                                                     >  1:3922716/1‑61 (MQ=255)
        aaGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAt                                                           >  1:2152591/1‑71 (MQ=255)
         aGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtg                                                          >  1:2457938/1‑71 (MQ=255)
           tACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTa                                                              <  1:2073249/65‑1 (MQ=255)
            aCAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAt                                                           >  1:928801/1‑67 (MQ=255)
            aCAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAt                                                           >  1:2413393/1‑67 (MQ=255)
            aCAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAt                                                           >  1:2222813/1‑67 (MQ=255)
            aCAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAt                                                           >  1:118882/1‑67 (MQ=255)
            aCAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGCTTACTGACTAGAt                                                           >  1:303842/1‑67 (MQ=255)
             cAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGc                                                       <  1:2450463/70‑1 (MQ=255)
             cAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGc                                                       <  1:2160350/70‑1 (MQ=255)
             cAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGc                                                       <  1:3280643/70‑1 (MQ=255)
               aCTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGCAGt                                                    <  1:2123136/71‑1 (MQ=255)
                     ccAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGCAGTTCCTGc                                              <  1:2720713/71‑1 (MQ=255)
                     ccAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGCAGTTCCTGc                                              <  1:1510228/71‑1 (MQ=255)
                      cAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGCAGTTCCTGc                                              <  1:3546465/70‑1 (MQ=255)
                       aGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGCAGTTCCTGCaa                                            >  1:1529015/1‑71 (MQ=255)
                       aGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGCAGTTCCTGCaa                                            >  1:358252/1‑71 (MQ=255)
                       aGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGCAGTTCCTGCaa                                            >  1:807199/1‑71 (MQ=255)
                        gTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtgtg                                                        <  1:2993018/58‑1 (MQ=255)
                        gTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtgtg                                                        <  1:429526/58‑1 (MQ=255)
                         ttCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACt                                          <  1:1517368/71‑1 (MQ=255)
                          tCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTa                                                                     >  1:2852111/1‑43 (MQ=255)
                          tCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTa                                                                     >  1:3731951/1‑43 (MQ=255)
                          tCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTa                                                                     >  1:1117152/1‑43 (MQ=255)
                          tCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTa                                                                     >  1:1797104/1‑43 (MQ=255)
                          tCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCa                                                      <  1:1130188/58‑1 (MQ=255)
                          tCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCa                                                      <  1:1198723/58‑1 (MQ=255)
                          tCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCa                                                      <  1:609597/58‑1 (MQ=255)
                          tCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGCGCa                                                      >  1:1488340/1‑58 (MQ=255)
                           cATTTGAAGTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCaa                                            >  1:2565937/1‑67 (MQ=255)
                           cATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtc                                        >  1:1772486/1‑71 (MQ=255)
                           cATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtc                                        >  1:3387678/1‑71 (MQ=255)
                           cATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtc                                        >  1:291938/1‑71 (MQ=255)
                           cATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtc                                        >  1:1182642/1‑71 (MQ=255)
                           cATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtc                                        >  1:2331305/1‑71 (MQ=255)
                           cATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtc                                        >  1:2926109/1‑71 (MQ=255)
                           cATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtc                                        >  1:783695/1‑71 (MQ=255)
                           cATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGCAGTTCCTGCAACTtc                                        >  1:712745/1‑71 (MQ=255)
                             tttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCa                                                      <  1:523508/55‑1 (MQ=255)
                             tttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCa                                                      <  1:310935/55‑1 (MQ=255)
                             tttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCa                                                      <  1:276699/55‑1 (MQ=255)
                             tttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtctc                                      <  1:1188760/71‑1 (MQ=255)
                             tttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtctc                                      <  1:129797/71‑1 (MQ=255)
                               tGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCtt                                    >  1:2229973/1‑71 (MQ=255)
                                gAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGCAGTTCCTGCAACTTCTCtt                                    <  1:1091760/70‑1 (MQ=255)
                                  aCTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCg                                 >  1:1544082/1‑71 (MQ=255)
                                   cTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAg                                                             <  1:163367/42‑1 (MQ=255)
                                   cTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAg                                                             <  1:3735591/42‑1 (MQ=255)
                                   cTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAg                                                             <  1:1952386/42‑1 (MQ=255)
                                   cTGGCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAg                                                             <  1:3000402/42‑1 (MQ=255)
                                     ggCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCgg                                >  1:1806525/1‑69 (MQ=255)
                                      gCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                               <  1:1527406/69‑1 (MQ=255)
                                      gCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                               <  1:4055965/69‑1 (MQ=255)
                                      gCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                               <  1:3096161/69‑1 (MQ=255)
                                      gCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                               <  1:1693412/69‑1 (MQ=255)
                                      gCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAg                             >  1:2137279/1‑71 (MQ=255)
                                      gCTTTTTTCAATTAATTGTGAAGATAGTTCACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAg                             >  1:2096926/1‑71 (MQ=255)
                                           tttCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGcc                         <  1:3756337/70‑1 (MQ=255)
                                           tttCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGcc                         <  1:2822192/70‑1 (MQ=255)
                                           tttCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGcc                         <  1:2869375/70‑1 (MQ=255)
                                            ttCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAg                       >  1:2403418/1‑71 (MQ=255)
                                            ttCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAg                       >  1:656504/1‑71 (MQ=255)
                                              cAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGt                      <  1:646741/70‑1 (MQ=255)
                                               tattaattGTGAAGATAGTTTACTGACTAGATGTGCa                                                      <  1:1833321/36‑1 (MQ=255)
                                                attaattGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTctct                                     <  1:2154667/53‑1 (MQ=255)
                                                attaattGTGAAGATAGTTCACTGACTAGATGTGCAGTTCCTGc                                              <  1:3102505/44‑1 (MQ=255)
                                                 ttaattGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGAACTTCTCTTTCGGCAGTGCCAGttcttc                 <  1:2481084/71‑1 (MQ=255)
                                                   aattGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttcttc                 <  1:3260968/70‑1 (MQ=255)
                                                   aattGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttcttc                 <  1:1396952/70‑1 (MQ=255)
                                                   aattGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCACttcttc                 <  1:3366788/70‑1 (MQ=255)
                                                    attGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttcttc                 <  1:1318322/69‑1 (MQ=255)
                                                    attGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttcttc                 <  1:1659450/69‑1 (MQ=255)
                                                    attGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCACttcttc                 <  1:3371565/69‑1 (MQ=255)
                                                         gAAGATAGTTTACTGACTAGATGTGCAGTTCCTGc                                              <  1:2625602/35‑1 (MQ=255)
                                                         gAAGATAGTTTACTGACTAGATGTGCAGTTCCTGc                                              <  1:636801/35‑1 (MQ=255)
                                                         gAAGATAGTTTACTGACTAGATGTGCAGTTCCTGc                                              <  1:290636/35‑1 (MQ=255)
                                                         gAAGATAGTTTACTGACTAGATGTGCAGTTCCTGc                                              <  1:3846460/35‑1 (MQ=255)
                                                         gAAGATAGTTTACTGACTAGATGTGCAGTTCCTGc                                              <  1:1297490/35‑1 (MQ=255)
                                                         gAAGATAGTTTACTGACTAGATGTGCAGTTCCTGc                                              <  1:1901661/35‑1 (MQ=255)
                                                          aaGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTg          <  1:3822649/70‑1 (MQ=255)
                                                          aaGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTg          <  1:892506/70‑1 (MQ=255)
                                                          aaGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTg          <  1:2555205/70‑1 (MQ=255)
                                                           aGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGtt        <  1:1811646/71‑1 (MQ=255)
                                                           aGATAGTTCACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGtt        <  1:2583622/71‑1 (MQ=255)
                                                           aGATAGTTCACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGtt        <  1:2030352/71‑1 (MQ=255)
                                                              tAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGAt     <  1:2895173/71‑1 (MQ=255)
                                                               aGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATg    >  1:921258/1‑71 (MQ=255)
                                                                gTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGc   >  1:479090/1‑71 (MQ=255)
                                                                gTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGc   >  1:2763003/1‑71 (MQ=255)
                                                                 tttACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCg  <  1:1945608/71‑1 (MQ=255)
                                                                 tttACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCg  <  1:3199350/71‑1 (MQ=255)
                                                                 tttACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCg  <  1:2989850/71‑1 (MQ=255)
                                                                  ttACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCg  >  1:2664645/1‑70 (MQ=255)
                                                                   |                                                                    
TGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCG  >  minE/1807310‑1807445

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: