Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,823,545 C→T 10.2% intergenic (+290/+668) nrfG → / ← pheU heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,823,5450CT10.2% 394.0 / 27.3 127intergenic (+290/+668)nrfG/pheUheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base C (59/55);  new base T (4/9);  total (63/64)
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.93e-01

Only 100 of 130 total aligned reads displayed.

GATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAA  >  minE/2823475‑2823615
                                                                      |                                                                      
gATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCAc                                                                        >  1:4058435/1‑71 (MQ=255)
  tACTGATTATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGt                                                                      <  1:1916114/71‑1 (MQ=255)
  tACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGt                                                                      <  1:3499193/71‑1 (MQ=255)
  tACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGt                                                                      <  1:3503747/71‑1 (MQ=255)
   aCTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGtt                                                                     >  1:3671552/1‑71 (MQ=255)
    cTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtc                                                                    <  1:2027349/71‑1 (MQ=255)
    cTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtc                                                                    <  1:432796/71‑1 (MQ=255)
      gATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc                                                                  >  1:2354456/1‑71 (MQ=255)
      gATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc                                                                  >  1:2214661/1‑71 (MQ=255)
      gATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc                                                                  >  1:3427788/1‑71 (MQ=255)
      gATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc                                                                  >  1:3529516/1‑71 (MQ=255)
      gATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc                                                                  >  1:3763905/1‑71 (MQ=255)
       aTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcg                                                                 <  1:2538196/71‑1 (MQ=255)
         cATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcgcg                                                               >  1:1644123/1‑71 (MQ=255)
         cATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcgcg                                                               >  1:3308283/1‑71 (MQ=255)
           tAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGtt                                                             <  1:633314/71‑1 (MQ=255)
            aaGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGttt                                                            >  1:3268138/1‑71 (MQ=255)
             aGCGTAAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTc                                                           >  1:3628114/1‑71 (MQ=255)
              gCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtct                                                                   >  1:3774409/1‑62 (MQ=255)
              gCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtct                                                                   >  1:3047985/1‑62 (MQ=255)
               cGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTg                                                         <  1:3752381/71‑1 (MQ=255)
               cGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTg                                                         <  1:3897787/71‑1 (MQ=255)
                  tAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTg                                                         <  1:1544242/68‑1 (MQ=255)
                   aaaaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAAc                                                      <  1:3978730/70‑1 (MQ=255)
                   aaaaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAAc                                                      <  1:2174451/70‑1 (MQ=255)
                   aaaaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACg                                                     >  1:4055124/1‑71 (MQ=255)
                   aaaaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACg                                                     >  1:2100090/1‑71 (MQ=255)
                   aaaaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACg                                                     >  1:3147156/1‑71 (MQ=255)
                     aaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACggg                                                   >  1:210522/1‑71 (MQ=255)
                     aaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACgg                                                    >  1:822238/1‑70 (MQ=255)
                      aaaaTCTACAAGCCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACggg                                                   >  1:266021/1‑70 (MQ=255)
                      aaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCATGTTCTCGCGTTTCTg                                                         <  1:1603401/64‑1 (MQ=255)
                      aaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACgggg                                                  <  1:218033/71‑1 (MQ=255)
                      aaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACgggg                                                  <  1:2921417/71‑1 (MQ=255)
                      aaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACgggg                                                  <  1:2758953/71‑1 (MQ=255)
                       aaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACgggg                                                  >  1:760620/1‑70 (MQ=255)
                       aaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACgggg                                                  <  1:622080/70‑1 (MQ=255)
                        aaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGaa                                                <  1:976160/71‑1 (MQ=255)
                         aTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAAc                                               <  1:2621571/71‑1 (MQ=255)
                            tACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACggg                                                   >  1:855253/1‑64 (MQ=255)
                            tACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACgg                                                    >  1:955528/1‑63 (MQ=255)
                              cAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCATGTTCTCGCGTTTCTGAACGGGGAACGgcgc                                          >  1:738761/1‑71 (MQ=255)
                              cAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGt                                                              >  1:1715858/1‑51 (MQ=255)
                              cAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGgcgc                                          >  1:551641/1‑71 (MQ=255)
                              cAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGgcg                                           >  1:3166126/1‑70 (MQ=255)
                               aaaCCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCt                                         >  1:2612330/1‑71 (MQ=255)
                                aaCCAACGCAACACAATTCATGCCCTGGCAGTATGTCATGATCTCGCGTTTCTGAACGGGGAACGGCGCTc                                        <  1:2510709/71‑1 (MQ=255)
                                aaCCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTc                                        <  1:1490552/71‑1 (MQ=255)
                                aaCCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTc                                        <  1:4016135/71‑1 (MQ=255)
                                 aCCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACgggg                                                  >  1:1018360/1‑60 (MQ=255)
                                 aCCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTc                                        >  1:2470254/1‑70 (MQ=255)
                                 aCCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTc                                        >  1:3501839/1‑70 (MQ=255)
                                 aCCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTc                                        >  1:1986461/1‑70 (MQ=255)
                                 aCCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTc                                        >  1:3591768/1‑70 (MQ=255)
                                 aCCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTc                                        >  1:1081709/1‑70 (MQ=255)
                                 aCCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGCGGAACGGCGCTc                                        >  1:3729594/1‑70 (MQ=255)
                                  ccAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCa                                      <  1:1264433/71‑1 (MQ=255)
                                  ccAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCa                                      <  1:3009886/71‑1 (MQ=255)
                                  ccAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCa                                      <  1:1117315/71‑1 (MQ=255)
                                  ccAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCa                                      <  1:1002968/71‑1 (MQ=255)
                                    aaCGCAACACAATTCATGCCCTGGCAGTATGTCATGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCAtt                                    <  1:4121494/71‑1 (MQ=255)
                                    aaCGCAACACAATTCATGCCCTGGCAGTATGTCATGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCAtt                                    <  1:2559555/71‑1 (MQ=255)
                                    aaCGCAACACAATTCATGCCCTGGCAGTATGTCATGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCAtt                                    <  1:2670042/71‑1 (MQ=255)
                                    aaCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCAtt                                    <  1:2541755/71‑1 (MQ=255)
                                    aaCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCAtt                                    <  1:2584250/71‑1 (MQ=255)
                                        cAACACAATTCATGCCCTGGCAGTATGTCATGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAgg                                <  1:3242920/71‑1 (MQ=255)
                                        cAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAgg                                <  1:1210074/71‑1 (MQ=255)
                                         aaCACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAgg                                <  1:2785468/70‑1 (MQ=255)
                                            acaATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGt                            <  1:782513/71‑1 (MQ=255)
                                             caATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGt                            >  1:311703/1‑70 (MQ=255)
                                                 tCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGtca                          >  1:3845481/1‑68 (MQ=255)
                                                 tCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGtca                          >  1:2107352/1‑68 (MQ=255)
                                                   aTGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGtcattca                      <  1:170640/70‑1 (MQ=255)
                                                    tGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGtcattcat                     >  1:738175/1‑70 (MQ=255)
                                                     gCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGaa                                                       >  1:2845424/1‑35 (MQ=255)
                                                     gCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGaa                                                       >  1:801575/1‑35 (MQ=255)
                                                       ccTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCataata                 <  1:177133/71‑1 (MQ=255)
                                                        cTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCataata                 <  1:3866043/70‑1 (MQ=255)
                                                        cTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCataa                   >  1:1429523/1‑68 (MQ=255)
                                                         tGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCataa                   <  1:3511102/67‑1 (MQ=255)
                                                         tGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAg                >  1:2297445/1‑70 (MQ=255)
                                                         tGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAg                >  1:3861540/1‑70 (MQ=255)
                                                           gCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCataata                 <  1:2293055/67‑1 (MQ=255)
                                                           gCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGa               <  1:1510578/69‑1 (MQ=255)
                                                            cAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAAc            <  1:1621662/71‑1 (MQ=255)
                                                               tATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAg                <  1:2690303/64‑1 (MQ=255)
                                                                aTGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCAt        >  1:453458/1‑71 (MQ=255)
                                                                aTGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCAt        >  1:159237/1‑71 (MQ=255)
                                                                 tGTCATGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCAtt       >  1:1055872/1‑71 (MQ=255)
                                                                 tGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGa               <  1:3895740/63‑1 (MQ=255)
                                                                 tGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGa               <  1:2134919/63‑1 (MQ=255)
                                                                   tCATGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGaaa             >  1:1659133/1‑63 (MQ=255)
                                                                   tCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCa                                      <  1:3669158/38‑1 (MQ=255)
                                                                   tCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCg     <  1:2510198/71‑1 (MQ=255)
                                                                   tCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCg     <  1:2993098/71‑1 (MQ=255)
                                                                   tCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCg     <  1:893090/71‑1 (MQ=255)
                                                                   tCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCg     <  1:3934699/71‑1 (MQ=255)
                                                                     aCGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGaa   <  1:1659651/71‑1 (MQ=255)
                                                                      cGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCataata                 <  1:2818764/56‑1 (MQ=255)
                                                                      cGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGaaa  <  1:1869328/71‑1 (MQ=255)
                                                                      |                                                                      
GATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAA  >  minE/2823475‑2823615

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: