Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,911,785 C→A 12.4% intergenic (+80/+126) 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,7850CA12.4% 209.5 / 11.7 72intergenic (+80/+126)mgtA/yjgFmagnesium transporter/ketoacid‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base A (4/5);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

GGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCAGGCGAT  >  Exported/2911717‑2911855
                                                                    |                                                                      
ggCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCtt                                                                      <  1:2150770/71‑1 (MQ=255)
   aCGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATc                                                                   <  1:1393271/71‑1 (MQ=255)
     gATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATcc                                                                  >  1:5592202/1‑70 (MQ=255)
     gATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATcc                                                                  <  1:5329835/70‑1 (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)
         cttatcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATccggat                                                              >  1:1035075‑M2/49‑68 (MQ=255)
         cttatcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATccggat                                                              >  1:2375641‑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)
                     aTTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAAt                                                  >  1:1179855/1‑70 (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:3166951‑M2/41‑1 (MQ=255)
                                                   ggatgcGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGc                                           <  1:2635218‑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:4256496/71‑1 (MQ=255)
                                                         ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcg             <  1:4635045/71‑1 (MQ=255)
                                                         ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcg             <  1:4971433/71‑1 (MQ=255)
                                                         ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcg             <  1:3069909/71‑1 (MQ=255)
                                                          gcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTa                               >  1:1219835/1‑52 (MQ=255)
                                                          gcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTa                               >  1:5058795/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:3934798/71‑1 (MQ=255)
                                                           cgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCa           >  1:3306600/1‑71 (MQ=255)
                                                           cgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCa           <  1:4497950/71‑1 (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:5577476/1‑60 (MQ=255)
                                                               aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCg                  >  1:5483771/1‑60 (MQ=255)
                                                               aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCa           >  1:4581059/1‑67 (MQ=255)
                                                               aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCAt          >  1:1544242/1‑68 (MQ=255)
                                                               aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCa        <  1:4790980/70‑1 (MQ=255)
                                                               aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCa        <  1:3865948/70‑1 (MQ=255)
                                                               aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCa        <  1:4840835/70‑1 (MQ=255)
                                                               aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCAg       >  1:5323613/1‑71 (MQ=255)
                                                               aaTGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCAg       >  1:2146472/1‑71 (MQ=255)
                                                                 tGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGaa                                      >  1:5377857/1‑38 (MQ=255)
                                                                 tGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGaa                                      >  1:434458/1‑38 (MQ=255)
                                                                 tGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGaa                                      >  1:3930220/1‑38 (MQ=255)
                                                                 tGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCAGGc     <  1:4230306/71‑1 (MQ=255)
                                                                  gCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACg                                 >  1:80636/1‑42 (MQ=255)
                                                                   ccTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCAGGCg    <  1:4190975/70‑1 (MQ=255)
                                                                    cTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCAGGCg    <  1:3237394/69‑1 (MQ=255)
                                                                    cTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCAGGCGAt  <  1:2466022/71‑1 (MQ=255)
                                                                    |                                                                      
GGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCATCAGGCGAT  >  Exported/2911717‑2911855

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: