Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,584,322 C→T 5.6% L985L (CTG→CTA) ‡ mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,584,3220CT 5.6% 360.9 / 7.1 107L985L (CTG→CTA) ‡mdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base C (39/62);  new base T (3/3);  total (42/65)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 110 total aligned reads displayed.

CCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGA  >  minE/2584257‑2584392
                                                                 |                                                                      
cccATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCAc                                                                   >  1:1215541/1‑71 (MQ=255)
 ccATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGgcag                                                                      <  1:696123/67‑1 (MQ=255)
      aCGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                                                >  1:1760819/1‑68 (MQ=255)
      aCGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                                                >  1:3895554/1‑68 (MQ=255)
      aCGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGcc                                                                >  1:4143680/1‑68 (MQ=255)
      aCGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCgag                                                             <  1:2300597/71‑1 (MQ=255)
      aCGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCgag                                                             <  1:3409624/71‑1 (MQ=255)
        gCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGaa                                                           >  1:2773088/1‑71 (MQ=255)
          cGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAAt                                                          <  1:1080967/70‑1 (MQ=255)
          cGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATa                                                         <  1:776450/71‑1 (MQ=255)
           ggTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATaa                                                        <  1:2386074/71‑1 (MQ=255)
               cctaccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAgg                                                     <  1:2376324/70‑1 (MQ=255)
               cctaccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAgg                                                     <  1:1136916/70‑1 (MQ=255)
                ctaccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATa                                                         <  1:1114923/65‑1 (MQ=255)
                  accGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGc                                                    >  1:3288401/1‑68 (MQ=255)
                   ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                                                <  1:3885301/71‑1 (MQ=255)
                   ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                                                <  1:3681596/71‑1 (MQ=255)
                   ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                                                <  1:1428751/71‑1 (MQ=255)
                   ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCa                                                  >  1:2490451/1‑69 (MQ=255)
                   ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCa                                                  <  1:3667066/69‑1 (MQ=255)
                     gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAg                                                      <  1:2171797/63‑1 (MQ=255)
                     gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagag                                               <  1:1834526/70‑1 (MQ=255)
                     gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                                                >  1:2647922/1‑69 (MQ=255)
                     gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                                                >  1:2652649/1‑69 (MQ=255)
                      cgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagagag                                             >  1:4116409/1‑71 (MQ=255)
                      cgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagagag                                             <  1:639438/71‑1 (MQ=255)
                       gTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAgg                                            >  1:3922321/1‑71 (MQ=255)
                       gTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAgg                                            >  1:492515/1‑71 (MQ=255)
                          ttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAgg                                            <  1:2368476/68‑1 (MQ=255)
                          ttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtca                                         >  1:3018587/1‑71 (MQ=255)
                          ttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtc                                          <  1:2590390/70‑1 (MQ=255)
                          ctGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGtc                                          >  1:3042924/2‑70 (MQ=255)
                           tGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtca                                         <  1:3132587/70‑1 (MQ=255)
                           tGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtca                                         <  1:3297886/70‑1 (MQ=255)
                           tGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtca                                         <  1:3691148/70‑1 (MQ=255)
                           tGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAAt                                                          <  1:839299/53‑1 (MQ=255)
                           tGCGCGCCACTGCCGGCACCGTTACTGATACCTAGCGGTAATACGCCGAGAAt                                                          <  1:3768472/53‑1 (MQ=255)
                            gcgcgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtcat                                        <  1:3285327/70‑1 (MQ=255)
                              gcgcCAGAACCGGCACCATGACTTATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAg                                     <  1:73211/71‑1 (MQ=255)
                              gcgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAg                                     <  1:3677551/71‑1 (MQ=255)
                               cgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                                    <  1:2766735/71‑1 (MQ=255)
                                gcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                                    <  1:3376530/70‑1 (MQ=255)
                                  cAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGGCATCAGGAtt                                 <  1:1965840/71‑1 (MQ=255)
                                    gAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTgg                               <  1:828765/71‑1 (MQ=255)
                                      aCCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATa                                                         <  1:2270056/43‑1 (MQ=255)
                                      aCCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGc                             >  1:1064371/1‑71 (MQ=255)
                                       ccGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAGGACg                            <  1:1490653/71‑1 (MQ=9)
                                       ccGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                                    <  1:1672802/63‑1 (MQ=255)
                                       ccGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCg                            >  1:119297/1‑71 (MQ=255)
                                        cGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                                    <  1:3487752/62‑1 (MQ=255)
                                         ggCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGt                           >  1:3151704/1‑70 (MQ=255)
                                         ggCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa                          <  1:3973594/71‑1 (MQ=255)
                                         ggCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa                          <  1:1305174/71‑1 (MQ=255)
                                            aCCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                        <  1:3193860/70‑1 (MQ=255)
                                            aCCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAAc                       <  1:2938334/71‑1 (MQ=255)
                                             ccATGACTGATAACCAGCGGCAGCACGCCGAGAATa                                                         <  1:382435/36‑1 (MQ=255)
                                             ccATGACTGATAACCAGCGGCAGCACGCCGAGAATa                                                         <  1:2306178/36‑1 (MQ=255)
                                             ccATGACTGATAACCAGCGGCAGCACGCCGAGAATa                                                         <  1:138680/36‑1 (MQ=255)
                                             ccATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagag                                               <  1:3230313/46‑1 (MQ=255)
                                             ccATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGt                           <  1:1491381/66‑1 (MQ=255)
                                             ccATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa                          <  1:325839/67‑1 (MQ=255)
                                             ccATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACg                      <  1:479912/71‑1 (MQ=255)
                                              cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGt                                           <  1:4049633/49‑1 (MQ=255)
                                              cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGt                                           <  1:2112082/49‑1 (MQ=255)
                                              cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                     <  1:460122/71‑1 (MQ=255)
                                              cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                     <  1:1497753/71‑1 (MQ=255)
                                               aTGACTGATAACCAGCGGCAGCACGCCGAGAATaaa                                                       >  1:2730724/1‑36 (MQ=255)
                                               aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAt                                  >  1:1367422/1‑57 (MQ=255)
                                               aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTggg                              >  1:4004351/1‑61 (MQ=255)
                                               aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                     <  1:3738684/70‑1 (MQ=255)
                                               aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCa                    >  1:3108791/1‑71 (MQ=255)
                                               aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCa                    >  1:2226244/1‑71 (MQ=255)
                                               aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCa                    <  1:1152603/71‑1 (MQ=255)
                                                tGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtca                                         <  1:1936658/49‑1 (MQ=255)
                                                tGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                     <  1:3213435/69‑1 (MQ=255)
                                                 taCTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAGGACGCAGACGCAta                  >  1:4021534/2‑70 (MQ=2)
                                                  aCTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATAGGATAGGACGCAGACGCATAcg                >  1:1001754/1‑71 (MQ=255)
                                                  aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagagag                                             >  1:3195361/1‑43 (MQ=255)
                                                  aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGa                                   >  1:4127276/1‑53 (MQ=255)
                                                  aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAt                                  >  1:155850/1‑54 (MQ=255)
                                                  aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                        >  1:27258/1‑64 (MQ=255)
                                                  aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                        >  1:2387376/1‑64 (MQ=255)
                                                  aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATcc                 >  1:102082/1‑71 (MQ=255)
                                                    tGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAt                                  <  1:2689535/52‑1 (MQ=255)
                                                    tGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCcg                <  1:3504817/70‑1 (MQ=255)
                                                     gATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                                                <  1:1826979/37‑1 (MQ=255)
                                                     gATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGCGTAAACGCATCcgc               <  1:2455645/69‑1 (MQ=255)
                                                      aTAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAtt                                 >  1:1133546/1‑51 (MQ=255)
                                                      aTAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACg                      >  1:1974473/1‑62 (MQ=255)
                                                      aTAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCcgc               <  1:616185/69‑1 (MQ=255)
                                                       tAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGcgcc            >  1:2377341/1‑71 (MQ=255)
                                                        aaCCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCa                    >  1:2717271/1‑62 (MQ=255)
                                                         aCCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGcgccgc          <  1:342431/71‑1 (MQ=255)
                                                            aGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCtt        >  1:3384175/1‑70 (MQ=255)
                                                            aGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTc       <  1:3313112/71‑1 (MQ=255)
                                                            aGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTc       <  1:2644862/71‑1 (MQ=255)
                                                            aGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTc       <  1:656111/71‑1 (MQ=255)
                                                             gaggCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCg      <  1:487447/69‑1 (MQ=255)
                                                              cGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                     >  1:3607859/1‑55 (MQ=255)
                                                                 cagcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCgatga  <  1:1594923/71‑1 (MQ=255)
                                                                 |                                                                      
CCCATCACGCCGGTACCTACCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGA  >  minE/2584257‑2584392

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: