Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,833,483 C→A 8.3% T266N (ACC→AAC)  groL → Cpn60 chaperonin GroEL, large subunit of GroESL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,833,4830CA 8.3% 203.6 / 4.0 60T266N (ACC→AAC) groLCpn60 chaperonin GroEL, large subunit of GroESL
Reads supporting (aligned to +/- strand):  ref base C (30/25);  new base A (2/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

GCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTG  >  Exported/2833417‑2833553
                                                                  |                                                                      
gCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCAt                                                                     >  1:845421/1‑70 (MQ=255)
    aCCGCTGCTGATCATCGCTGAAGATGTAGAAGGCTAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTg                                                                <  1:1103263/71‑1 (MQ=255)
    aCCGCTGCTGATCATCGCTGAAGATGTAGAAGGCTAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTg                                                                <  1:1106701/71‑1 (MQ=255)
          gctgATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCg                                                          >  1:1291864/1‑71 (MQ=255)
                catcGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGa                                                       >  1:2040291/1‑68 (MQ=255)
                    gCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGc                                                >  1:608405/1‑71 (MQ=255)
                       gAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGt                                                         <  1:139691/59‑1 (MQ=255)
                         aGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCgg                                           >  1:1806915/1‑71 (MQ=255)
                          gATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGt                                          <  1:2311175/71‑1 (MQ=255)
                           aTGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCAt                                                                     <  1:810851/43‑1 (MQ=255)
                           aTGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCAt                                                                     <  1:2878379/43‑1 (MQ=255)
                           aTGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCAt                                                                     <  1:1715907/43‑1 (MQ=255)
                            tGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTa                                        <  1:187616/71‑1 (MQ=255)
                             gTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTg                                                        >  1:411087/1‑54 (MQ=255)
                             gTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTg                                                        >  1:2997129/1‑54 (MQ=255)
                              tAGAAGGCGAAGCGCTGTCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTaa                                       <  1:1683449/70‑1 (MQ=255)
                              tAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTaa                                       <  1:345681/70‑1 (MQ=255)
                              tAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTaa                                       <  1:1725063/70‑1 (MQ=255)
                              tAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTaa                                       <  1:593785/70‑1 (MQ=255)
                               aGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAg                                     >  1:502148/1‑71 (MQ=255)
                               aGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAg                                     >  1:2846411/1‑71 (MQ=255)
                                gAAGGCGAAGCGCTGGCGACCCTGGTCGTCAACAACATCCGTGGCAt                                                            >  1:1387551/1‑47 (MQ=19)
                                gAAGGCGAAGCGCTGGCGACCCTGGTCGTCAACAACATCCGTGGCAt                                                            >  1:1089322/1‑47 (MQ=19)
                                gAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATg                                                                    <  1:745170/39‑1 (MQ=255)
                                gAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCAt                                                            >  1:1379698/1‑47 (MQ=255)
                                gAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGa                                                       >  1:628008/1‑52 (MQ=255)
                                gAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGa                                                       >  1:1731578/1‑52 (MQ=255)
                                   ggCGAAGCGCTGGCGACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCAcc                                 >  1:2033550/1‑71 (MQ=255)
                                     cGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCgg                               <  1:1491823/71‑1 (MQ=255)
                                        aGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCt                            >  1:1038669/1‑71 (MQ=255)
                                        aGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCt                            >  1:1354474/1‑71 (MQ=255)
                                         gcgcTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCtt                           >  1:1673445/1‑71 (MQ=255)
                                         gcgcTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCtt                           <  1:1935551/71‑1 (MQ=255)
                                          cgcTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTg                                                        <  1:2124970/41‑1 (MQ=255)
                                            cTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGaa                                                      >  1:2152445/1‑41 (MQ=255)
                                             tGGCCACCCTGGTGGTGAACAACATGCGTGGCATCGTCAAGGTCGCGGCTGTCaa                                       <  1:1465786/55‑1 (MQ=255)
                                             tGGCCACCCTGGTGGTGAACAACATGCGTGGCATCGTCAAGGTCGCGGCTGTCaa                                       <  1:142543/55‑1 (MQ=255)
                                             tGGCCACCCTGGTGGTGAACAACATGCGTGGCATCGTCAAGGTCGCGGCTGTCaa                                       <  1:796575/55‑1 (MQ=255)
                                               gCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCAc                                  <  1:2051839/58‑1 (MQ=255)
                                               gCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCg                      <  1:2293606/70‑1 (MQ=255)
                                               gCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGa                     <  1:1774300/71‑1 (MQ=255)
                                                 aaCTCTGGTTGTTAACACCATGCGTGGCATCGTGaa                                                      >  1:447672/1‑36 (MQ=255)
                                                   ctctGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCg                                            <  1:1239338/44‑1 (MQ=255)
                                                   ctctGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGAtcgt                 >  1:879183/1‑71 (MQ=255)
                                                      tGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAg                                     >  1:2603096/1‑48 (MQ=255)
                                                      tGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAg                                     >  1:1078707/1‑48 (MQ=255)
                                                      tGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCa                                   >  1:1490926/1‑50 (MQ=255)
                                                      tGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCAc                                  >  1:253955/1‑51 (MQ=255)
                                                      tGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCAc                                  >  1:1567147/1‑51 (MQ=255)
                                                      tGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCAc                                  >  1:1570567/1‑51 (MQ=255)
                                                      tGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGAGGTTAAAg                                     >  1:1426785/1‑48 (MQ=255)
                                                       ggTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTa                                        <  1:1974172/44‑1 (MQ=255)
                                                        gttgttAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTaa            <  1:1530232/71‑1 (MQ=255)
                                                         ttgttAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTaa            <  1:710124/70‑1 (MQ=255)
                                                           gttAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAg          >  1:2446047/1‑70 (MQ=255)
                                                           gttAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAg          >  1:1639359/1‑70 (MQ=255)
                                                           gttAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAg          <  1:2899974/70‑1 (MQ=255)
                                                           gttAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAg          <  1:2389247/70‑1 (MQ=255)
                                                            ttAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGc         >  1:1153287/1‑70 (MQ=255)
                                                             tAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTa       >  1:1126426/1‑71 (MQ=255)
                                                               acacCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGc                       <  1:1704901/53‑1 (MQ=255)
                                                                cacCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCa                                   <  1:1457729/40‑1 (MQ=255)
                                                                cacCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTAtgc    >  1:3049942/1‑71 (MQ=255)
                                                                cacCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTAtgc    >  1:271459/1‑71 (MQ=255)
                                                                 acCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTAtgct   >  1:566288/1‑71 (MQ=255)
                                                                  ccATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTAtgctg  <  1:1180969/71‑1 (MQ=255)
                                                                  |                                                                      
GCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTG  >  Exported/2833417‑2833553

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: