Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,833,469 T→C 9.6% T261T (ACT→ACC 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,4690TC 9.6% 146.1 / 7.2 52T261T (ACT→ACCgroLCpn60 chaperonin GroEL, large subunit of GroESL
Reads supporting (aligned to +/- strand):  ref base T (25/22);  new base C (2/3);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGT  >  Exported/2833407‑2833538
                                                              |                                                                     
gCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGgttg                                                                <  1:1222525/70‑1 (MQ=255)
 ccAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGgttgt                                                               <  1:1453087/70‑1 (MQ=255)
    aaGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTaa                                                            >  1:323911/1‑70 (MQ=255)
    aaGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTaa                                                            >  1:2281238/1‑70 (MQ=255)
     aGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAaca                                                          >  1:284449/1‑71 (MQ=255)
     aGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAaca                                                          >  1:2138933/1‑71 (MQ=255)
          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)
                     ctgATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAac                                                           <  1:557932/54‑1 (MQ=255)
                       gATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTaa                                                            >  1:106234/1‑51 (MQ=255)
                       gATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTaa                                                            >  1:637397/1‑51 (MQ=255)
                          catcGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGa                                        >  1:2040291/1‑68 (MQ=255)
                              gCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGc                                 >  1:608405/1‑71 (MQ=255)
                                 gAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGt                                          <  1:139691/59‑1 (MQ=255)
                                  aaGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTaa                                                            >  1:2734223/1‑40 (MQ=255)
                                   aGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCgg                            >  1:1806915/1‑71 (MQ=255)
                                    gATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGt                           <  1:2311175/71‑1 (MQ=255)
                                     aTGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCAt                                                      <  1:1715907/43‑1 (MQ=255)
                                     aTGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCAt                                                      <  1:2878379/43‑1 (MQ=255)
                                     aTGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCAt                                                      <  1:810851/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:593785/70‑1 (MQ=255)
                                        tAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTaa                        <  1:345681/70‑1 (MQ=255)
                                        tAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTaa                        <  1:1725063/70‑1 (MQ=255)
                                         aGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAg                      >  1:502148/1‑71 (MQ=255)
                                         aGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAg                      >  1:2846411/1‑71 (MQ=255)
                                          gAAGGCGAAGCGCTGGCGACCCTGGTCGTCAACAACATCCGTGGCAt                                             >  1:1089322/1‑47 (MQ=19)
                                          gAAGGCGAAGCGCTGGCGACCCTGGTCGTCAACAACATCCGTGGCAt                                             >  1:1387551/1‑47 (MQ=19)
                                          gAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATg                                                     <  1:745170/39‑1 (MQ=255)
                                          gAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCAt                                             >  1:1379698/1‑47 (MQ=255)
                                          gAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGa                                        >  1:1731578/1‑52 (MQ=255)
                                          gAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGa                                        >  1:628008/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:796575/55‑1 (MQ=255)
                                                       tGGCCACCCTGGTGGTGAACAACATGCGTGGCATCGTCAAGGTCGCGGCTGTCaa                        <  1:142543/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)
                                                              |                                                                     
GCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGT  >  Exported/2833407‑2833538

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: