Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,396,435 C→G 8.8% E251D (GAG→GAC) ‡ corA ← magnesium/nickel/cobalt transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,396,4350CG 8.8% 269.3 / 13.9 91E251D (GAG→GAC) ‡corAmagnesium/nickel/cobalt transporter
Reads supporting (aligned to +/- strand):  ref base C (31/52);  new base G (3/5);  total (34/57)
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.48e-01

GAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTGGAACAGGGATTCGT  >  Exported/2396367‑2396499
                                                                    |                                                                
gAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTc                                                                <  1:2526692/71‑1 (MQ=255)
gAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGAAGATGCGGTTCTGCTc                                                                <  1:4773092/71‑1 (MQ=255)
  aGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGa                                                              <  1:2167132/71‑1 (MQ=255)
  aGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGCTCTGCTCGa                                                              <  1:3130823/71‑1 (MQ=255)
   gCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGAt                                                             <  1:3354515/71‑1 (MQ=255)
            gTCGGCGGCAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATgaa                                                    <  1:1813713/71‑1 (MQ=255)
            gTCGGCGGCAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATgaa                                                    <  1:2447158/71‑1 (MQ=255)
            gTCGGCGGCAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATgaa                                                    <  1:5057904/71‑1 (MQ=255)
            gTCGGCGGCAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATgaa                                                    <  1:5072763/71‑1 (MQ=255)
            gTCGGCGGCAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATgaa                                                    <  1:549962/71‑1 (MQ=255)
             tCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATaaaa                                                   >  1:258502/1‑71 (MQ=255)
               ggcggcAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATgaagc                                                  >  1:4217374/1‑68 (MQ=255)
               ggcggcAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATGAAGcc                                                 >  1:4184570/1‑71 (MQ=255)
               ggcggcAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATGAAGcc                                                 >  1:1692191/1‑71 (MQ=255)
                   gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGAt                                                       <  1:1047516/61‑1 (MQ=255)
                   gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGAt                                                       <  1:3384193/61‑1 (MQ=255)
                   gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGAt                                                       <  1:3557998/61‑1 (MQ=255)
                   gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGAt                                                       <  1:4141456/61‑1 (MQ=255)
                   gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGAt                                                       <  1:1777772/61‑1 (MQ=255)
                   gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGAt                                                       <  1:1228375/61‑1 (MQ=255)
                   gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAtt                                             <  1:1087850/71‑1 (MQ=255)
                   gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAtt                                             <  1:2477986/71‑1 (MQ=255)
                    cAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt                                              >  1:3800745/1‑69 (MQ=255)
                    cAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt                                              >  1:812618/1‑69 (MQ=255)
                     aGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATaa                                                     >  1:251984/1‑61 (MQ=255)
                       gAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt                                              >  1:1695809/1‑66 (MQ=255)
                       gAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt                                              >  1:1675692/1‑66 (MQ=255)
                       gAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt                                              >  1:2162886/1‑66 (MQ=255)
                         ataccacggataccacggAGAAGATTTTGATGCTGCGGTTCTGCTCGATGTTGATAAAAccc                                                >  1:2891766/1‑62 (MQ=255)
                         aTACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATa                                                      >  1:4469958/1‑56 (MQ=255)
                         aTACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAAccc                                                >  1:4330856/1‑62 (MQ=255)
                          tACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCt                                      >  1:1174252/1‑71 (MQ=255)
                           aCCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGc                                                                  <  1:1085906/42‑1 (MQ=255)
                           aCCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGc                                                                  <  1:3525462/42‑1 (MQ=255)
                           aCCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGc                                                                  <  1:699964/42‑1 (MQ=255)
                           aCCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGc                                                                  <  1:1383812/42‑1 (MQ=255)
                            ccACGGATACCACCGAGAAGATTTTTATGATGCTGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGc                                    <  1:2257654/71‑1 (MQ=255)
                            ccACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCTATGTTGATAAAACCCATTGCCGCCTGc                                    <  1:4529187/71‑1 (MQ=255)
                            ccACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGc                                    <  1:4311684/71‑1 (MQ=255)
                            ccACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGc                                    <  1:824486/71‑1 (MQ=255)
                             cACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATg                                                            <  1:1131899/46‑1 (MQ=255)
                             cACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATg                                                            <  1:1138150/46‑1 (MQ=255)
                             cACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATg                                                            <  1:2826970/46‑1 (MQ=255)
                              aCGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTc                                                                <  1:2991302/41‑1 (MQ=255)
                              aCGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCAt                                  >  1:4366902/1‑71 (MQ=255)
                                ggATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCAt                                  >  1:2249070/1‑69 (MQ=255)
                                ggATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCAt                                  >  1:3429463/1‑69 (MQ=255)
                                  aTACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATaaaa                                                   >  1:773147/1‑50 (MQ=255)
                                  aTACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt                                              >  1:4984013/1‑55 (MQ=255)
                                  aTACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt                                              >  1:2268030/1‑55 (MQ=255)
                                  aTACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt                                              >  1:1645515/1‑55 (MQ=255)
                                  aTACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt                                              >  1:3312330/1‑55 (MQ=255)
                                  aTACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt                                              >  1:1741983/1‑55 (MQ=255)
                                  aTACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt                                              >  1:3814045/1‑55 (MQ=255)
                                  aTACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAgg                              <  1:2766635/71‑1 (MQ=255)
                                  aTACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAgg                              <  1:642244/71‑1 (MQ=255)
                                  aTACCACCGAGAAGATTTAGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAgg                              <  1:698915/71‑1 (MQ=255)
                                   tACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGa                             <  1:4877332/71‑1 (MQ=255)
                                   tACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGa                             <  1:4751950/71‑1 (MQ=255)
                                     ccaccGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAg                           <  1:1530667/71‑1 (MQ=255)
                                     ccaccGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAg                           <  1:3926750/71‑1 (MQ=255)
                                     ccaccGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAACCCATTGCCGCCTGCATCAGGAAg                           <  1:3034435/70‑1 (MQ=255)
                                          gagaAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGa                             <  1:1948/64‑1 (MQ=255)
                                            gaagaTTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAg                           >  1:3563058/1‑64 (MQ=255)
                                            gaagaTTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAg                           >  1:5023836/1‑64 (MQ=255)
                                            gaagaTTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCa                       <  1:822/68‑1 (MQ=255)
                                            gaagaTTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCa                       <  1:1360881/68‑1 (MQ=255)
                                             aagaTTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCtt                   <  1:577200/71‑1 (MQ=255)
                                             aagaTTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCtt                   <  1:742254/71‑1 (MQ=255)
                                             aagaTTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCtt                   <  1:1644396/71‑1 (MQ=255)
                                             aagaTTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCtt                   <  1:1391931/71‑1 (MQ=255)
                                             aagaTTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCtt                   <  1:3319705/71‑1 (MQ=255)
                                              agaTTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCtt                   <  1:4397177/70‑1 (MQ=255)
                                               gaTTTTGATGATGCGGTTCTGCTCGATGTTGATaaaa                                                   >  1:2043189/1‑37 (MQ=255)
                                               gaTTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGAATCAGGAAGTTCACCTTCt                 >  1:4182288/1‑71 (MQ=255)
                                                aTTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTg                <  1:3919714/71‑1 (MQ=255)
                                                 ttttGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTgg               <  1:169153/71‑1 (MQ=255)
                                                 ttttGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTgg               <  1:2935005/71‑1 (MQ=255)
                                                 ttttGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTgg               <  1:1596639/71‑1 (MQ=255)
                                                 ttttGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTgg               <  1:1821522/71‑1 (MQ=255)
                                                 ttttGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTg                <  1:2144787/70‑1 (MQ=255)
                                                 ttttGATGATGCGGTGCTGCGCGATGGTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTgg               <  1:3284914/71‑1 (MQ=255)
                                                  tttGATGATGCGGTTCTGCTCGATGTTGATAAAAccc                                                >  1:654501/1‑37 (MQ=255)
                                                         atgAGGTTCTGCTCGATGTTGATAAAACCCATTgccg                                         >  1:4550234/1‑37 (MQ=255)
                                                          tgCGGTTCTGCTCGATGTTGATAAAACCCATTGCGGCCTGCATCAGGAAGTTCACCTTCTGGAACAGGGAt      >  1:4505098/1‑71 (MQ=255)
                                                          tgCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTGGAACAGGGAt      >  1:3676968/1‑71 (MQ=255)
                                                             ggTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTGGAACAGGGATTCg   <  1:4660083/71‑1 (MQ=255)
                                                             ggTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTGGAACAGGGATTCg   >  1:4380997/1‑71 (MQ=255)
                                                             ggTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTGGAACAGGGATTCg   <  1:1909333/71‑1 (MQ=255)
                                                             ggTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTGGAACAGGGATTCg   >  1:2176394/1‑71 (MQ=255)
                                                             ggTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTGGAACAGGGATTCg   <  1:2305148/71‑1 (MQ=255)
                                                             ggTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTGGAACAGGGATTCg   <  1:2781951/71‑1 (MQ=255)
                                                              gTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTGGAACAGGGATTCg   <  1:2133329/70‑1 (MQ=255)
                                                              gTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTGGAACAGGGATTCGt  <  1:2239164/71‑1 (MQ=255)
                                                                  tGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAgg                              <  1:2354409/39‑1 (MQ=255)
                                                                    |                                                                
GAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCTGCATCAGGAAGTTCACCTTCTGGAACAGGGATTCGT  >  Exported/2396367‑2396499

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