Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,376,834 G→A 10.6% V518V (GTC→GTT ubiB ← 2‑octaprenylphenol hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,376,8340GA10.6% 262.6 / 17.1 94V518V (GTC→GTTubiB2‑octaprenylphenol hydroxylase
Reads supporting (aligned to +/- strand):  ref base G (44/40);  new base A (4/6);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

ACCGACAAACCAGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATAACGCGATTGTCCCT  >  minE/2376765‑2376901
                                                                     |                                                                   
aCCGACAAACCAGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGa                                                                    >  1:1076328/1‑71 (MQ=255)
 ccGACAAACCAGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGAc                                                                   >  1:660685/1‑71 (MQ=255)
 ccGACAAACCAGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTAAc                                                                   <  1:2912042/71‑1 (MQ=255)
  cGACAAACCAGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGAc                                                                   >  1:1209584/1‑70 (MQ=255)
  cGACAAACCAGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGAc                                                                   >  1:2355985/1‑70 (MQ=255)
  cGACAAACCAGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGAc                                                                   >  1:3190425/1‑70 (MQ=255)
      aaaCCAGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTAACcaac                                                               >  1:1483493/1‑70 (MQ=255)
        aCCAGGCGATCAGACCACCTGCCATTAACCAGCCTGGCATCAGCCCCCATTCAGGTCGGCTGAcc                                                                  <  1:646073/65‑1 (MQ=255)
        aCCAGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACcaac                                                               >  1:142498/1‑68 (MQ=255)
         ccAGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAg                                                            >  1:3224215/1‑70 (MQ=255)
         ccAGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGa                                                           <  1:820120/71‑1 (MQ=255)
          cAGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGaa                                                          >  1:1598724/1‑71 (MQ=255)
           aGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACcaa                                                                >  1:1296192/1‑64 (MQ=255)
           aGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAAt                                                         <  1:1357748/71‑1 (MQ=255)
           aGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAAt                                                         <  1:2583495/71‑1 (MQ=255)
           aGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAAt                                                         <  1:566805/71‑1 (MQ=255)
            ggCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTAACCAACAAGAAt                                                         <  1:695354/70‑1 (MQ=255)
             gCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAAtgt                                                       <  1:2991693/71‑1 (MQ=255)
             gCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAAtgt                                                       >  1:3330334/1‑71 (MQ=255)
              ggATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACAAACAAGAAtgtg                                                      >  1:3621317/2‑71 (MQ=255)
              cGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAAtgtg                                                      >  1:2302760/1‑71 (MQ=255)
               gATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGc                                                     >  1:782933/1‑71 (MQ=255)
                 tCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCa                                                   <  1:3387727/71‑1 (MQ=255)
                 tCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCa                                                   >  1:548547/1‑71 (MQ=255)
                 tCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCa                                                   <  1:1090802/71‑1 (MQ=255)
                  cAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGcc                                                    >  1:2337749/1‑69 (MQ=255)
                   aGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCAc                                                  <  1:3081840/70‑1 (MQ=255)
                      ccaccTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTaa                                              >  1:2492698/1‑71 (MQ=255)
                          cTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATAc                                           >  1:2997705/1‑70 (MQ=255)
                           tGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTaa                                              <  1:3406727/66‑1 (MQ=255)
                           tGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATa                                            <  1:1777576/68‑1 (MQ=255)
                           tGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACt                                          >  1:3248619/1‑70 (MQ=255)
                            gCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTg                                                                     <  1:1714773/42‑1 (MQ=255)
                            gCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGc                                                     <  1:268733/58‑1 (MQ=255)
                             ccATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAAc                                       >  1:744029/1‑71 (MQ=255)
                              cATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACg                                      <  1:1831394/71‑1 (MQ=255)
                               aTTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGt                                     <  1:220094/71‑1 (MQ=255)
                               aTTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTAACCAACAAGAATGTGCCACTTAATACTAACg                                      <  1:39706/70‑1 (MQ=255)
                                ttAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGcc                                                    <  1:2815010/55‑1 (MQ=255)
                                ttAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGcc                                                    <  1:3206477/55‑1 (MQ=255)
                                ttAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACg                                      >  1:22290/1‑69 (MQ=255)
                                 tAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAg                                   >  1:442796/1‑71 (MQ=255)
                                  aaCCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCAc                                                  >  1:3315866/1‑55 (MQ=255)
                                  aaCCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAgc                                  >  1:40883/1‑71 (MQ=255)
                                  aaCCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAgc                                  >  1:633655/1‑71 (MQ=255)
                                   aCCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAgc                                  >  1:2719519/1‑70 (MQ=255)
                                   aCCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTAACcaaca                                                              >  1:2756807/1‑42 (MQ=255)
                                    ccAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGt                                     <  1:1097428/66‑1 (MQ=255)
                                     cAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAAtgt                                                       <  1:2601259/47‑1 (MQ=255)
                                      aGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACtt                                                <  1:834709/53‑1 (MQ=255)
                                        ccGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAAt                            <  1:3879842/71‑1 (MQ=255)
                                        ccGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAAt                            >  1:1333437/1‑71 (MQ=255)
                                        ccGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAAt                            >  1:1885501/1‑71 (MQ=255)
                                         cGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAAt                            <  1:3724555/70‑1 (MQ=255)
                                           ggCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTcc                         <  1:2900870/71‑1 (MQ=255)
                                           ggCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTcc                         >  1:2072088/1‑71 (MQ=255)
                                           ggCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTc                          <  1:2993359/70‑1 (MQ=255)
                                           ggCATCAGCCCCCATTCAGGGCGGCTGACCAACAAGAATGTGCCACTTAATa                                            <  1:2817367/52‑1 (MQ=255)
                                            gCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCg                        <  1:2257459/71‑1 (MQ=255)
                                              aTCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACt                                                 <  1:2398414/44‑1 (MQ=255)
                                              aTCAGCCCCCATTCAGGTCGGCTAACCAACAAGAATGTGCCACTTAATACTAACg                                      <  1:1684554/55‑1 (MQ=255)
                                               tCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCAc                                                  >  1:3795629/1‑42 (MQ=255)
                                                cAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCga                       >  1:1752356/1‑68 (MQ=255)
                                                cAGCCCCCATTCAGGTCGGCTAACCAACAAGAATGTGCCAc                                                  >  1:3311697/1‑41 (MQ=255)
                                                  gCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTa                                         <  1:1648902/48‑1 (MQ=255)
                                                  gCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTa                                         <  1:1063357/48‑1 (MQ=255)
                                                  gCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTa                                         <  1:3365737/48‑1 (MQ=255)
                                                   cccccATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAAt                 <  1:2192071/71‑1 (MQ=255)
                                                   cccccATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAAATCCGAGAAAAt                 <  1:1931429/71‑1 (MQ=255)
                                                    ccccATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAAt                 >  1:2180242/1‑70 (MQ=255)
                                                    ccccATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAAt                 >  1:632777/1‑70 (MQ=255)
                                                     cccATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGcc                               <  1:3250435/55‑1 (MQ=255)
                                                      ccATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAAt                 <  1:3879401/68‑1 (MQ=255)
                                                       cATTCAGGTCGGCTGCCCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATaa               >  1:2206865/1‑69 (MQ=255)
                                                       cATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTc                          >  1:1118364/1‑58 (MQ=255)
                                                       cATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTc                          >  1:1419566/1‑58 (MQ=255)
                                                       cATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATaa               >  1:2761074/1‑69 (MQ=255)
                                                       cATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATAAcg             <  1:3128775/71‑1 (MQ=255)
                                                        aTTCAGGTCGGCTGACCAACAAGAATGTGCCACTTaa                                              <  1:2850804/37‑1 (MQ=255)
                                                        aTTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATAAcgc            <  1:3607018/71‑1 (MQ=255)
                                                        aTTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATAAcgc            >  1:2782224/1‑71 (MQ=255)
                                                         ttCAGGTCGGCTAACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTcc                         <  1:2197869/57‑1 (MQ=255)
                                                           cAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCga                       >  1:2000905/1‑57 (MQ=255)
                                                           cAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCga                       >  1:3742082/1‑57 (MQ=255)
                                                           cAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATAACGCGa          <  1:2229776/70‑1 (MQ=255)
                                                           cAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATAACGCGa          <  1:2160009/70‑1 (MQ=255)
                                                           cAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATAACGCGAt         >  1:569773/1‑71 (MQ=255)
                                                              gTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATAACGCGATTGt      <  1:1219613/71‑1 (MQ=255)
                                                               tCGGCTAACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATAACGCGATTGTc     <  1:3757994/71‑1 (MQ=255)
                                                                cGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCga                       >  1:2556840/1‑52 (MQ=255)
                                                                cGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATAACGCGATTGTcc    >  1:2398275/1‑71 (MQ=255)
                                                                 ggCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATAACGCGATTGTccc   >  1:3005672/1‑71 (MQ=255)
                                                                  gCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATAACGCGATTGTCCCt  >  1:1565699/1‑71 (MQ=255)
                                                                  gCTAACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATAACGCGATTGTCCCt  >  1:339775/1‑71 (MQ=255)
                                                                   cTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGaaaa                  <  1:399163/54‑1 (MQ=255)
                                                                     |                                                                   
ACCGACAAACCAGGCGATCAGACCACCTGCCATTAACCAGCCGGGCATCAGCCCCCATTCAGGTCGGCTGACCAACAAGAATGTGCCACTTAATACTAACGTAGCGCCAATTCCGAGAAAATAACGCGATTGTCCCT  >  minE/2376765‑2376901

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: