Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,911,781 A→G 15.0% intergenic (+76/+130) mgtA → / ← yjgF magnesium transporter/ketoacid‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,911,7810AG15.0% 137.1 / 15.6 60intergenic (+76/+130)mgtA/yjgFmagnesium transporter/ketoacid‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (20/31);  new base G (4/5);  total (24/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AAACTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCAG  >  Exported/2911711‑2911850
                                                                      |                                                                     
aaaCTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaaa                                                                       >  1:100364/1‑71 (MQ=255)
aaaCTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaaa                                                                       >  1:4000063/1‑71 (MQ=255)
 aaCTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAAt                                                                      >  1:2340116/1‑71 (MQ=255)
  aCTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATg                                                                     <  1:4279937/71‑1 (MQ=255)
      ggCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCtt                                                                 <  1:2150770/71‑1 (MQ=255)
         aCGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATc                                                              <  1:1393271/71‑1 (MQ=255)
           gATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATcc                                                             <  1:5329835/70‑1 (MQ=255)
           gATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATcc                                                             >  1:5592202/1‑70 (MQ=255)
           gATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATc                                                              <  1:1027106/69‑1 (MQ=255)
           gATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCg                                                            >  1:2046015/1‑71 (MQ=255)
            aTAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCgg                                                           >  1:2608133/1‑71 (MQ=255)
             tAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGc                                                          >  1:5080350/1‑71 (MQ=255)
              aGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaaa                                                                       >  1:1647753/1‑57 (MQ=255)
               cttatcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATccggat                                                         >  1:5580942‑M2/49‑68 (MQ=255)
               cttatcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATccggat                                                         >  1:2375641‑M2/49‑68 (MQ=255)
               cttatcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATccggat                                                         >  1:1035075‑M2/49‑68 (MQ=255)
                  ccGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCcaaca                                                     <  1:3010212/71‑1 (MQ=255)
                  ccGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATACGGCGTAAATGCCTTATCCGGCcaac                                                      >  1:664712/1‑70 (MQ=255)
                   tcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCGGAtaacaa                                                    <  1:591027‑M2/27‑1 (MQ=255)
                   tcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCGGAtaacaa                                                    <  1:548195‑M2/27‑1 (MQ=255)
                   tcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCGGAtaacaa                                                    <  1:1471091‑M2/27‑1 (MQ=255)
                    gATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACaaa                                                   <  1:3773129/71‑1 (MQ=255)
                     gtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCgg                                                           <  1:2777510‑M2/20‑1 (MQ=255)
                     gtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCGgataacaata                                                  >  1:179359‑M2/43‑69 (MQ=255)
                     gtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCGgataacaata                                                  <  1:1359697‑M2/29‑3 (MQ=255)
                       tCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGc                                                                    <  1:406676/51‑1 (MQ=255)
                           aTTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAAt                                             >  1:1179855/1‑70 (MQ=255)
                                 gACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaaa                                                                       >  1:4444367/1‑38 (MQ=255)
                                 gACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaaa                                                                       >  1:5045404/1‑38 (MQ=255)
                                 gACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTg                                       <  1:3545159/70‑1 (MQ=255)
                                   cAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAg                                    >  1:1643894/1‑71 (MQ=255)
                                       gTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCgg                                                           >  1:1000726/1‑44 (MQ=255)
                                       gTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAAt                                <  1:4649130/71‑1 (MQ=255)
                                         gAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATcc                                                             >  1:1803372/1‑40 (MQ=255)
                                          aGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCCGCCAACAAATTCAATAAATTGCAGGAATCAc                             >  1:2271938/1‑71 (MQ=255)
                                           gTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGcc                                                         >  1:3416890/1‑42 (MQ=255)
                                           gTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCa                              >  1:2281783/1‑69 (MQ=255)
                                              ttGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAg                         >  1:812365/1‑71 (MQ=255)
                                                  tCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCt                     >  1:3438460/1‑71 (MQ=255)
                                                   cATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCt                     >  1:3526411/1‑70 (MQ=255)
                                                   cATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCt                     <  1:4446270/70‑1 (MQ=255)
                                                    aTGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGa                   <  1:5099439/71‑1 (MQ=255)
                                                    aTGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGa                   <  1:2773408/71‑1 (MQ=255)
                                                         ggatgcGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGc                                      <  1:2635218‑M2/41‑1 (MQ=255)
                                                         ggatgcGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGc                                      <  1:3166951‑M2/41‑1 (MQ=255)
                                                            tgcGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCg             >  1:5076195‑M2/4‑69 (MQ=255)
                                                            tgcGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTa           <  1:1182871‑M2/68‑1 (MQ=255)
                                                            tgcGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTa           <  1:1901803‑M2/68‑1 (MQ=255)
                                                            tgcGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTa           <  1:4959696‑M2/68‑1 (MQ=255)
                                                              cGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGc         <  1:5184338‑M2/70‑1 (MQ=255)
                                                               ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAAt                                <  1:4642394/47‑1 (MQ=255)
                                                               ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcg        <  1:3069909/71‑1 (MQ=255)
                                                               ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcg        <  1:4635045/71‑1 (MQ=255)
                                                               ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcg        <  1:4256496/71‑1 (MQ=255)
                                                               ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcg        <  1:4971433/71‑1 (MQ=255)
                                                                gcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTa                          >  1:1219835/1‑52 (MQ=255)
                                                                gcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTa                          >  1:5058795/1‑52 (MQ=255)
                                                                gcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcgc       <  1:3327339/71‑1 (MQ=255)
                                                                gcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcgc       <  1:2051417/71‑1 (MQ=255)
                                                                 cgTAAATGCCTTATCCTGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCa      <  1:1377405/71‑1 (MQ=255)
                                                                 cgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCa      <  1:4497950/71‑1 (MQ=255)
                                                                 cgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCa      <  1:3934798/71‑1 (MQ=255)
                                                                 cgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCa      >  1:3306600/1‑71 (MQ=255)
                                                                   tAAATGCCTTATCCGGCCAATAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATc    <  1:3104206/71‑1 (MQ=255)
                                                                   tAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATc    <  1:1343368/71‑1 (MQ=255)
                                                                     aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAg                         >  1:3781377/1‑48 (MQ=255)
                                                                     aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCg             >  1:5483771/1‑60 (MQ=255)
                                                                     aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCg             >  1:5577476/1‑60 (MQ=255)
                                                                     aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCa      >  1:4581059/1‑67 (MQ=255)
                                                                     aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCAt     >  1:1544242/1‑68 (MQ=255)
                                                                     aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCa   <  1:3865948/70‑1 (MQ=255)
                                                                     aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCa   <  1:4840835/70‑1 (MQ=255)
                                                                     aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCa   <  1:4790980/70‑1 (MQ=255)
                                                                     aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCAg  >  1:2146472/1‑71 (MQ=255)
                                                                     aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCAg  >  1:5323613/1‑71 (MQ=255)
                                                                      |                                                                     
AAACTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCAG  >  Exported/2911711‑2911850

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: