Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,911,779 A→C 15.0% intergenic (+74/+132) 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,7790AC15.0% 120.8 / 20.2 60intergenic (+74/+132)mgtA/yjgFmagnesium transporter/ketoacid‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (20/31);  new base C (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.10e-01

ATAAACTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATC  >  Exported/2911709‑2911848
                                                                      |                                                                     
aTAAACTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                       <  1:3290706/71‑1 (MQ=255)
  aaaCTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaaa                                                                     >  1:100364/1‑71 (MQ=255)
  aaaCTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaaa                                                                     >  1:4000063/1‑71 (MQ=255)
   aaCTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaa                                                                      >  1:1414919/1‑69 (MQ=255)
   aaCTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAAt                                                                    >  1:2340116/1‑71 (MQ=255)
    aCTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATg                                                                   <  1:4279937/71‑1 (MQ=255)
       gggCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                       <  1:1312315/64‑1 (MQ=255)
        ggCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCtt                                                               <  1:2150770/71‑1 (MQ=255)
           aCGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaa                                                                      >  1:532283/1‑61 (MQ=255)
           aCGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATc                                                            <  1:1393271/71‑1 (MQ=255)
            cGATAGCCCCGATTCAGATTGTAGACAAAGTGAGATTTGTTCATGCCGGATGCGGCGTa                                                                       <  1:4294125/59‑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:2375641‑M2/49‑68 (MQ=255)
                 cttatcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATccggat                                                       >  1:1035075‑M2/49‑68 (MQ=255)
                 cttatcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATccggat                                                       >  1:5580942‑M2/49‑68 (MQ=255)
                    ccGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCcaaca                                                   <  1:3010212/71‑1 (MQ=255)
                    ccGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATACGGCGTAAATGCCTTATCCGGCcaac                                                    >  1:664712/1‑70 (MQ=255)
                     tcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCGGAtaacaa                                                  <  1:548195‑M2/27‑1 (MQ=255)
                     tcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCGGAtaacaa                                                  <  1:1471091‑M2/27‑1 (MQ=255)
                     tcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCGGAtaacaa                                                  <  1:591027‑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)
                       aTTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                       <  1:2748961/48‑1 (MQ=255)
                       aTTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                       <  1:3686778/48‑1 (MQ=255)
                        ttCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                       <  1:2204925/47‑1 (MQ=255)
                         tCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                       <  1:396003/46‑1 (MQ=255)
                         tCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGc                                                                  <  1:406676/51‑1 (MQ=255)
                             aTTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAAt                                           >  1:1179855/1‑70 (MQ=255)
                                 tAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                       <  1:816051/38‑1 (MQ=255)
                                 tAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                       <  1:2207174/38‑1 (MQ=255)
                                   gACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaaa                                                                     >  1:5045404/1‑38 (MQ=255)
                                   gACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaaa                                                                     >  1:4444367/1‑38 (MQ=255)
                                   gACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTg                                     <  1:3545159/70‑1 (MQ=255)
                                    aCAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                       <  1:4234065/35‑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:2773408/71‑1 (MQ=255)
                                                      aTGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGa                 <  1:5099439/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:1901803‑M2/68‑1 (MQ=255)
                                                              tgcGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTa         <  1:1182871‑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:4635045/71‑1 (MQ=255)
                                                                 ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcg      <  1:4971433/71‑1 (MQ=255)
                                                                 ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcg      <  1:4256496/71‑1 (MQ=255)
                                                                 ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcg      <  1:3069909/71‑1 (MQ=255)
                                                                  gcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTa                        >  1:5058795/1‑52 (MQ=255)
                                                                  gcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTa                        >  1:1219835/1‑52 (MQ=255)
                                                                  gcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcgc     <  1:2051417/71‑1 (MQ=255)
                                                                  gcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcgc     <  1:3327339/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)
                                                                      |                                                                     
ATAAACTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATC  >  Exported/2911709‑2911848

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: