Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,603,330 C→T 5.4% P447P (CCG→CCA yhiP ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,603,3300CT 5.4% 403.6 / 8.1 128P447P (CCG→CCAyhiPpredicted transporter
Reads supporting (aligned to +/- strand):  ref base C (52/69);  new base T (5/2);  total (57/71)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 131 total aligned reads displayed.

TACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCT  >  minE/2603261‑2603400
                                                                     |                                                                      
tactGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCg                                                                       >  1:3109756/1‑71 (MQ=255)
tactGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCg                                                                       >  1:2925849/1‑71 (MQ=255)
tactGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCg                                                                       >  1:3049624/1‑71 (MQ=255)
 actGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCg                                                                       >  1:2538737/1‑70 (MQ=255)
  ctGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCgg                                                                      <  1:2361893/70‑1 (MQ=255)
  ctGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCgg                                                                      <  1:213395/70‑1 (MQ=255)
  ctGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCgg                                                                      <  1:574132/70‑1 (MQ=255)
  ctGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCgg                                                                      <  1:3517328/70‑1 (MQ=255)
  ctGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGa                                                                     >  1:4047623/1‑71 (MQ=255)
  ctGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGa                                                                     >  1:4383663/1‑71 (MQ=255)
     cAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGAt                                                                    <  1:1600018/69‑1 (MQ=255)
     cAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATc                                                                   <  1:1963115/70‑1 (MQ=255)
     cAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCg                                                                  >  1:4592499/1‑71 (MQ=255)
     cAACGCCCAGAGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCg                                                                  >  1:4285865/1‑71 (MQ=255)
       aCGTCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCgg                                                                 >  1:3296378/1‑70 (MQ=255)
       aCGCTCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGt                                                                <  1:444085/71‑1 (MQ=255)
       aCGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCg                                                                  >  1:3294781/1‑69 (MQ=255)
       aCGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGt                                                                <  1:3500555/71‑1 (MQ=255)
        cGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGt                                                                <  1:4476946/70‑1 (MQ=255)
        cGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTa                                                               <  1:2631430/71‑1 (MQ=255)
          cccAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAt                                                             >  1:4296075/1‑71 (MQ=255)
          cccAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAt                                                             <  1:1545224/71‑1 (MQ=255)
           ccAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtg                                                            <  1:2364543/71‑1 (MQ=255)
            cAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgt                                                           >  1:422128/1‑71 (MQ=255)
              gCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgttg                                                         <  1:2540708/71‑1 (MQ=255)
              gCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgttg                                                         <  1:2304686/71‑1 (MQ=255)
              gCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgtt                                                          >  1:1878381/1‑70 (MQ=255)
                gTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGAATCGGTAATGTTGTc                                                       <  1:1414867/71‑1 (MQ=255)
                  gaccagaccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCg                                                     >  1:2330875/1‑71 (MQ=255)
                  gaccagaccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCg                                                     >  1:3965520/1‑71 (MQ=255)
                   accagaccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCgg                                                    <  1:349810/71‑1 (MQ=255)
                    ccagaccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGc                                                   <  1:3702244/71‑1 (MQ=255)
                    caagaccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGc                                                   <  1:1186986/69‑1 (MQ=255)
                     cagaccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCa                                                  >  1:4038084/1‑71 (MQ=255)
                        accaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCgg                                                    <  1:968381/66‑1 (MQ=255)
                         ccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgt                                                           >  1:3441425/1‑58 (MQ=255)
                         ccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgt                                                           >  1:2841291/1‑58 (MQ=255)
                          caATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACcgcg                                             >  1:4551243/1‑71 (MQ=255)
                           aaTCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgtt                                                          <  1:3168855/57‑1 (MQ=255)
                              cTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                          >  1:818531/1‑70 (MQ=255)
                               ttACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGAt                                                                    >  1:2298824/1‑43 (MQ=255)
                                 aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaaa                                        <  1:1043519/69‑1 (MQ=255)
                                 aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaa                                         <  1:1648742/68‑1 (MQ=255)
                                 aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaa                                         <  1:1843319/68‑1 (MQ=255)
                                 aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                          <  1:4652112/67‑1 (MQ=255)
                                 aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                          <  1:1861925/67‑1 (MQ=255)
                                 aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                          <  1:1999554/67‑1 (MQ=255)
                                 aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAAt                                       <  1:184411/70‑1 (MQ=255)
                                  ccAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGt                                     >  1:4010549/1‑71 (MQ=255)
                                   cAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGc                                                   <  1:2068025/56‑1 (MQ=255)
                                   cAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                          <  1:3939507/65‑1 (MQ=255)
                                    aaaCACGTTGGTATAGACGGGCAACGTCTCAAGCgg                                                                      <  1:1545328/36‑1 (MQ=255)
                                     aaCACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGc                                  >  1:2288308/1‑71 (MQ=255)
                                      acacGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTaa                                                              >  1:2655052/1‑42 (MQ=255)
                                      acacGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTaa                                                              >  1:3888748/1‑42 (MQ=255)
                                      acacGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTaa                                                              >  1:2917394/1‑42 (MQ=255)
                                        acGTTGGTATAGACGGGCAACGTCTCAAGTGGATCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCcac                               <  1:93216/71‑1 (MQ=255)
                                        acGTTGGTATAGACGGGCAACGTCTCAAGTGGAGCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCcac                               <  1:3073503/71‑1 (MQ=255)
                                          gTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACAt                             <  1:822036/71‑1 (MQ=255)
                                          gTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACAt                             <  1:4420413/71‑1 (MQ=255)
                                          gTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACAt                             <  1:2887222/71‑1 (MQ=255)
                                              gTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                          <  1:2229406/54‑1 (MQ=255)
                                              gTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCca                                <  1:3162340/64‑1 (MQ=255)
                                               tataGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACTGCGGTAAATGTTGCcaca                              >  1:4573113/1‑65 (MQ=255)
                                               tataGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCcaca                              >  1:303159/1‑65 (MQ=255)
                                                 taGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaa                                         <  1:1111688/52‑1 (MQ=255)
                                                   gACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAgccgcc                      >  1:3978275/1‑69 (MQ=255)
                                                    aCGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAg                   >  1:2044725/1‑71 (MQ=255)
                                                        gCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAg               <  1:3267562/71‑1 (MQ=255)
                                                         cAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAg               >  1:4345756/1‑70 (MQ=255)
                                                         cAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGa              <  1:504914/71‑1 (MQ=255)
                                                          aaCGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATa                            >  1:607096/1‑56 (MQ=255)
                                                           aCGTCTCAAGTGGATCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCCACATAGCCGCCCAGCaa                >  1:2974053/1‑67 (MQ=255)
                                                           aCGTCTCAAGTGGATCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCCACATAGCCGCCCAGCaa                >  1:1377851/1‑67 (MQ=255)
                                                           aCGTCTCAAGTGGATCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCCACATAGCCGCCCAGCaa                >  1:1601761/1‑67 (MQ=255)
                                                           aCGTCTCAAGTGGATCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCCACATAGCCGCCCAGCaa                >  1:3062828/1‑67 (MQ=255)
                                                           aCGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                          <  1:4439310/41‑1 (MQ=255)
                                                           aCGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAAc            >  1:3017917/1‑71 (MQ=255)
                                                            cGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAAcg           <  1:3884345/71‑1 (MQ=255)
                                                             gTCTCAAGCGGATCGGTAATGTTGTCCGGCACCTCGGTAAATGTTGCCACATAGCCGCCCa                    >  1:365496/1‑61 (MQ=255)
                                                             gTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAAcgc          >  1:2094868/1‑71 (MQ=255)
                                                              tctcAAGTGGATCGGTAATGTTGTCTGGTACTGCGGTAAATGTTGCCACAt                             >  1:857072/1‑51 (MQ=255)
                                                              tctcAAGCGGATCGGTAATGTTGTCCGGCACCGCTGTa                                          <  1:3964316/38‑1 (MQ=255)
                                                              tctcAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                                          <  1:2470818/38‑1 (MQ=255)
                                                              tctcAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAAcgcg         >  1:1500796/1‑71 (MQ=255)
                                                              tctcAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAAcgcg         <  1:2851158/71‑1 (MQ=255)
                                                                tcAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAAt                                       >  1:4173192/1‑39 (MQ=255)
                                                                tcAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCgg        >  1:1823065/1‑70 (MQ=255)
                                                                tcAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCgg        >  1:1414880/1‑70 (MQ=255)
                                                                  aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCca                                >  1:272680/1‑44 (MQ=255)
                                                                  aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAg               <  1:2573887/61‑1 (MQ=255)
                                                                  aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAAcgcg         >  1:2387303/1‑67 (MQ=255)
                                                                  aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCa      <  1:853357/70‑1 (MQ=255)
                                                                  aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAg     >  1:3295283/1‑71 (MQ=255)
                                                                  aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAg     <  1:524230/71‑1 (MQ=255)
                                                                  aaGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAg     >  1:414878/1‑71 (MQ=255)
                                                                   aGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGc    <  1:582622/71‑1 (MQ=255)
                                                                     cGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCa      >  1:1918217/1‑67 (MQ=255)
                                                                     cGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCt  <  1:3155459/71‑1 (MQ=255)
                                                                     cGGATCGGTAATGTTGTCCGGCACCGCCGTAAATGTTg                                   >  1:114533/1‑38 (MQ=255)
                                                                     |                                                                      
TACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGCCGCCCAGCAAGAACGCGGCAGCCT  >  minE/2603261‑2603400

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: