Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,490,422 A→G 6.0% E107E (GAA→GAG recF → gap repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,490,4220AG 6.0% 236.8 / 4.3 83E107E (GAA→GAGrecFgap repair protein
Reads supporting (aligned to +/- strand):  ref base A (38/40);  new base G (3/2);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GCTTAACCAAAGACAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTCAACGGCGGCCCCAAA  >  Exported/2490352‑2490491
                                                                      |                                                                     
gCTTAACCAAAGACAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGaa                                                                       >  1:2962593/1‑71 (MQ=255)
    aaCCAAAGACAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTgg                                                                   >  1:2136791/1‑71 (MQ=255)
    aaCCAAAGACAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTgg                                                                   >  1:2034583/1‑71 (MQ=255)
    aaCCAAAGACAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTg                                                                    <  1:1984870/70‑1 (MQ=255)
    aaCCAAAGACAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTg                                                                    <  1:4770967/70‑1 (MQ=255)
     aCCAAAGACAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTgg                                                                   >  1:4984166/1‑70 (MQ=255)
     aCCAAAGACAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGgc                                                                  >  1:1320164/1‑71 (MQ=255)
       cAAAGACAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGgcgc                                                                <  1:3577684/71‑1 (MQ=255)
        aaaGACAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCa                                                               <  1:785923/71‑1 (MQ=255)
         aaGACAAACAGGGCGTCAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCAc                                                              >  1:1074693/1‑71 (MQ=255)
           gACAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCt                                                            <  1:584575/71‑1 (MQ=255)
           gACAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCt                                                            <  1:2447182/71‑1 (MQ=255)
           gACAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCt                                                            >  1:5429478/1‑71 (MQ=255)
             cAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGa                                                          >  1:1714750/1‑71 (MQ=255)
              aaaCAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGAt                                                         >  1:2447450/1‑71 (MQ=255)
              aaaCAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGAt                                                         >  1:3776579/1‑71 (MQ=255)
                 cAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAAc                                                                      >  1:2731703/1‑55 (MQ=255)
                 cAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGc                                                       >  1:270892/1‑70 (MQ=255)
                 cAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGc                                                       >  1:4465243/1‑70 (MQ=255)
                  aGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGcc                                                      >  1:5213531/1‑70 (MQ=255)
                    ggCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGc                                                       <  1:4109961/67‑1 (MQ=255)
                    ggCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGc                                                       <  1:3397314/67‑1 (MQ=255)
                    ggCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGc                                                       <  1:3258770/67‑1 (MQ=255)
                     gCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATg                                                  <  1:651961/71‑1 (MQ=255)
                     gCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATg                                                  <  1:5055305/71‑1 (MQ=255)
                      cGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGc                                                 <  1:3911793/71‑1 (MQ=255)
                      cGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGc                                                 <  1:1911686/71‑1 (MQ=255)
                      cGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGc                                                 <  1:4631496/71‑1 (MQ=255)
                       gACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATgaa                                                <  1:2589906/71‑3 (MQ=255)
                        aCAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAg                                               <  1:3987775/71‑1 (MQ=255)
                        aCAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAg                                               <  1:4818926/71‑1 (MQ=255)
                          aGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGtt                                             >  1:2223141/1‑71 (MQ=255)
                           gCATAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGgcgc                                                                <  1:3181597/51‑1 (MQ=255)
                           gCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTg                                            <  1:2185773/71‑1 (MQ=255)
                            cAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATg                                                  <  1:2253710/64‑1 (MQ=255)
                             aaaGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCAc                                                              >  1:3560169/1‑51 (MQ=255)
                             aaaGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGa                                           >  1:5543408/1‑70 (MQ=255)
                              aaGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGAt                                          >  1:3222839/1‑70 (MQ=255)
                               agcccGCAGCGACGGTACAGACGGGCATAAGGGCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATaa                                        <  1:2528336/68‑1 (MQ=255)
                                 tCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCt                                                            <  1:1348486/49‑1 (MQ=255)
                                      accGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGAt                                          <  1:918747/60‑1 (MQ=255)
                                       tCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGAt                                                         >  1:1333118/1‑46 (MQ=255)
                                       tCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGAt                                                         >  1:573505/1‑46 (MQ=255)
                                       tCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGAt                                          <  1:3176980/61‑1 (MQ=255)
                                        cGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCa                                                     <  1:5032303/49‑1 (MQ=255)
                                        cGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAg                               >  1:194331/1‑71 (MQ=255)
                                         gACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCAc                                                              >  1:2928928/1‑39 (MQ=255)
                                         gACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCAc                                                              >  1:5489064/1‑39 (MQ=255)
                                         gACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCAc                                                              >  1:1632622/1‑39 (MQ=255)
                                         gACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATaa                                        >  1:1185790/1‑61 (MQ=255)
                                         gACCGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAgg                              >  1:3269094/1‑71 (MQ=255)
                                          aCGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCa                                                     <  1:4876407/47‑1 (MQ=255)
                                          aCGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGa                                           <  1:5416725/57‑1 (MQ=255)
                                          aCGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGa                                           <  1:3590902/57‑1 (MQ=255)
                                           cGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCAACTGATGCCAATGc                                                 >  1:5569068/1‑50 (MQ=255)
                                           cGGTACAGACGGGCATAAAGTCGCGGAGCTGGCGCACCt                                                            >  1:2823387/1‑39 (MQ=38)
                                           cGGTACAGACGGGCATAAAGTCGCGGAGCTGGCGCACCt                                                            >  1:1756925/1‑39 (MQ=38)
                                           cGGTACAGACGGGCATAAAGTCGCGGAGCTGGCGCACCt                                                            >  1:4770837/1‑39 (MQ=38)
                                              tACAGACGGGCATAAGGTCGCGGAACTGGCGCACCt                                                            <  1:2987577/36‑1 (MQ=255)
                                               aCAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAg                                  >  1:3491037/1‑61 (MQ=255)
                                               aCAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAg                                  >  1:4812758/1‑61 (MQ=255)
                                                  gACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGtttacttt                     <  1:3643135/71‑1 (MQ=255)
                                                      ggCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTc                 <  1:1799218/71‑1 (MQ=255)
                                                      ggCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTc                 <  1:1526110/71‑1 (MQ=255)
                                                         aTAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTCAAc              <  1:2409904/71‑1 (MQ=255)
                                                          tAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTCAAcg             >  1:503179/1‑71 (MQ=255)
                                                          tAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTCAAcg             >  1:1326416/1‑71 (MQ=255)
                                                          tAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTCAAcg             <  1:2264923/71‑1 (MQ=255)
                                                           aaGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGa                                 <  1:1672298/50‑1 (MQ=255)
                                                           aaGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTCAAcgg            >  1:2575997/1‑71 (MQ=255)
                                                           aaGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTCAAcgg            >  1:1208613/1‑71 (MQ=255)
                                                            aGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGAt                                          <  1:1789410/40‑1 (MQ=255)
                                                            aGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGAt                                          <  1:4925006/40‑1 (MQ=255)
                                                            aGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGaa                                >  1:5604016/1‑50 (MQ=255)
                                                            aGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTCAAcggc           <  1:3045665/71‑1 (MQ=255)
                                                            aGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTCAAcggc           >  1:1857838/1‑71 (MQ=255)
                                                            aGGTCGCAGAACTGGCGCACCTGATGCCAATGCAGTTGATaa                                        >  1:1440043/1‑42 (MQ=255)
                                                              gTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAg                                  <  1:5356294/46‑1 (MQ=255)
                                                               tCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGttt                          >  1:243093/1‑53 (MQ=255)
                                                                cgcgGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTCAACGGCGGcc       >  1:3493140/1‑71 (MQ=255)
                                                                  cgGAGCTGGCGCACCTGATGCCAATGCAGCTGATAACGCCAGAAGGGttt                          <  1:639481/50‑1 (MQ=255)
                                                                  cgGAGCTGGCGCACCTGATGCCAATGCAGCTGATAACGCCAGAAGGGttt                          <  1:85521/50‑1 (MQ=255)
                                                                  cgGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGt                            <  1:434512/48‑1 (MQ=255)
                                                                    gAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTCAACGGCGGcccc     >  1:5476894/1‑69 (MQ=255)
                                                                     aaCTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTc                 <  1:5434073/56‑1 (MQ=255)
                                                                     aaCTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTCAACGGCGGCCCCaaa  >  1:2106698/1‑71 (MQ=255)
                                                                     aaCTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTCAACGGCGGCCCCaaa  <  1:3210843/71‑1 (MQ=255)
                                                                      |                                                                     
GCTTAACCAAAGACAAACAGGGCGACAGCAAAGTCCGCATCGACGGTACAGACGGGCATAAGGTCGCGGAACTGGCGCACCTGATGCCAATGCAGTTGATAACGCCAGAAGGGTTTACTTTACTCAACGGCGGCCCCAAA  >  Exported/2490352‑2490491

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: