Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 13,594 A→C 6.1% I478L (ATC→CTC)  dnaK → chaperone Hsp70, co‑chaperone with DnaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported13,5940AC 6.1% 417.2 / 19.8 148I478L (ATC→CTC) dnaKchaperone Hsp70, co‑chaperone with DnaJ
Reads supporting (aligned to +/- strand):  ref base A (60/79);  new base C (7/2);  total (67/81)
Fisher's exact test for biased strand distribution p-value = 7.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 158 total aligned reads displayed.

GGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCA  >  Exported/13524‑13663
                                                                      |                                                                     
gggTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGata                                                                       >  1:6261313/1‑71 (MQ=255)
  gTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                                     <  1:4324541/71‑1 (MQ=255)
  gTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                                     <  1:3321881/71‑1 (MQ=255)
   tCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                                     >  1:357250/1‑70 (MQ=255)
     aGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGa                                                                   <  1:6073067/70‑1 (MQ=255)
     aGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGAt                                                                  >  1:5015644/1‑71 (MQ=255)
      gTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATg                                                                 <  1:5629640/71‑1 (MQ=255)
      gTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATg                                                                 <  1:3606120/71‑1 (MQ=255)
        tCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                                     >  1:855763/1‑65 (MQ=255)
        tCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCt                                                               >  1:5241691/1‑71 (MQ=255)
         cAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATg                                                                 >  1:2522621/1‑68 (MQ=255)
         cAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCt                                                               >  1:6611379/1‑70 (MQ=255)
         cAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCt                                                               >  1:5902923/1‑70 (MQ=255)
            ccTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                                     <  1:2908979/61‑1 (MQ=255)
            ccTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACg                                                           >  1:6614363/1‑71 (MQ=255)
                 aTGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                                     <  1:4503383/56‑1 (MQ=255)
                 aTGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                                     <  1:6783266/56‑1 (MQ=255)
                 aTGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                                     <  1:5776140/56‑1 (MQ=255)
                 aTGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                                     <  1:6938292/56‑1 (MQ=255)
                 aTGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                                     <  1:5430815/56‑1 (MQ=255)
                   ggTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGTAGTTACCTTCGATATCGATg                                                                 <  1:4467372/58‑1 (MQ=255)
                   ggTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGa                                                                   >  1:2500743/1‑56 (MQ=255)
                   ggTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATg                                                                 <  1:1805074/58‑1 (MQ=255)
                   ggTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATg                                                                 <  1:4185801/58‑1 (MQ=255)
                   ggTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATg                                                                 <  1:2033667/58‑1 (MQ=255)
                       tCAACCCGTCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCAc                                                <  1:170479/71‑1 (MQ=255)
                       tCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCAc                                                <  1:597100/71‑1 (MQ=255)
                           cccGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGttt                                            <  1:22488/71‑1 (MQ=255)
                           cccGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGttt                                            >  1:5999024/1‑71 (MQ=255)
                            ccGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTc                                           <  1:1931801/71‑1 (MQ=255)
                            ccGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTc                                           <  1:6316488/71‑1 (MQ=255)
                                  ccGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGaa                                     <  1:4124565/71‑1 (MQ=255)
                                     cgcgGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGa                                      >  1:775265/1‑67 (MQ=255)
                                     cgcgGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGa                                      >  1:6926914/1‑67 (MQ=255)
                                       cgGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATa                                <  1:3690011/71‑1 (MQ=255)
                                        ggCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGa                                                             >  1:4993896/1‑41 (MQ=255)
                                         gCATGCCGCAGATCGAAGTTACCTGCGATATCTATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                              <  1:2897972/71‑1 (MQ=255)
                                          cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGt                                                         <  1:4556288/43‑1 (MQ=255)
                                          cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGt                                                         <  1:6820599/43‑1 (MQ=255)
                                          cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGt                                                         <  1:5059197/43‑1 (MQ=255)
                                          cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGaa                                     <  1:1779865/63‑1 (MQ=255)
                                          cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGaa                                     <  1:3107784/63‑1 (MQ=255)
                                          cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGaa                                     <  1:3073577/63‑1 (MQ=255)
                                          cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGaa                                     <  1:2253284/63‑1 (MQ=255)
                                          cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                              >  1:5561404/1‑70 (MQ=255)
                                          cATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                              >  1:4987412/1‑70 (MQ=255)
                                           aTGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                              <  1:2637055/69‑1 (MQ=255)
                                            tGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCTAAAGATaaaaa                            >  1:2704660/1‑70 (MQ=255)
                                            tGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaaaa                            >  1:1107382/1‑70 (MQ=255)
                                            tGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAAc                           <  1:6595622/71‑1 (MQ=255)
                                            tGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAAc                           <  1:467937/71‑1 (MQ=255)
                                            tGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAAc                           <  1:72073/71‑1 (MQ=255)
                                                  aGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCgg                      >  1:4430515/1‑70 (MQ=255)
                                                  aGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCg                       >  1:956456/1‑69 (MQ=255)
                                                   taTCGAAGGTAACTTCGATCTCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCgg                      >  1:322143/2‑69 (MQ=38)
                                                   taTCGAAGGTAACTTCGATCTCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCgg                      >  1:2266511/2‑69 (MQ=38)
                                                   taTCGAAGGTAACTTCGATCTCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGt                     >  1:3985335/2‑70 (MQ=38)
                                                   taTCGAAGGTAACTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGt                     >  1:1383877/2‑70 (MQ=255)
                                                   gATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTa                    <  1:6774982/71‑1 (MQ=255)
                                                    aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                                   <  1:6171405/55‑1 (MQ=255)
                                                    aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                                   <  1:1834339/55‑1 (MQ=255)
                                                    aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                                   <  1:531131/55‑1 (MQ=255)
                                                    aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                                   <  1:244803/55‑1 (MQ=255)
                                                    aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                                   <  1:258929/55‑1 (MQ=255)
                                                    aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                                   <  1:2607715/55‑1 (MQ=255)
                                                    aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                                   <  1:2622870/55‑1 (MQ=255)
                                                    aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                                   <  1:420839/55‑1 (MQ=255)
                                                    aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                                   <  1:6159224/55‑1 (MQ=255)
                                                    aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                              <  1:5041751/60‑1 (MQ=255)
                                                    aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                              <  1:5552084/60‑1 (MQ=255)
                                                    aTCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                              <  1:1901447/60‑1 (MQ=255)
                                                     tCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGt                     >  1:6875360/1‑68 (MQ=255)
                                                     tCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGt                     >  1:4455031/1‑68 (MQ=255)
                                                     tCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaaa                  <  1:6057743/71‑1 (MQ=255)
                                                      cGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCg                       >  1:1695142/1‑65 (MQ=255)
                                                      cGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTa                    >  1:326927/1‑68 (MQ=255)
                                                       gAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaaa                  >  1:3029722/1‑69 (MQ=255)
                                                       gAAGGTAACTTCGATCTCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAag                 <  1:476033/70‑1 (MQ=38)
                                                       gAAGGTAACTTCGATCTCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAag                 <  1:6641940/70‑1 (MQ=38)
                                                        aaGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaaa                  >  1:901120/1‑68 (MQ=255)
                                                        aaGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaaa                  >  1:804194/1‑68 (MQ=255)
                                                        aaGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaaa                  >  1:2500796/1‑68 (MQ=255)
                                                        aaGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTaa                   >  1:2008821/1‑67 (MQ=255)
                                                         aGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAagag               >  1:6718166/1‑70 (MQ=255)
                                                          gTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTa                    <  1:4283242/64‑1 (MQ=255)
                                                          gTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAagag               >  1:1403995/1‑69 (MQ=255)
                                                           ttACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGATAGATAAAAACAGCGGTa                    <  1:4535994/63‑1 (MQ=255)
                                                           ttACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCag            <  1:6436499/71‑1 (MQ=255)
                                                           ttACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCag            <  1:3980314/71‑1 (MQ=255)
                                                           ttACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCag            <  1:2702548/71‑1 (MQ=255)
                                                           ttACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCag            <  1:1442343/71‑1 (MQ=255)
                                                                ttCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaaaa                            <  1:4481160/50‑1 (MQ=255)
                                                                  cGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAg                                   <  1:2361673/41‑1 (MQ=255)
                                                                  cGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCa     >  1:4784074/1‑71 (MQ=255)
                                                                  cGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCa     >  1:1620370/1‑71 (MQ=255)
                                                                  cGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCa     >  1:5668017/1‑71 (MQ=255)
                                                                      atCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCg                       >  1:122970/1‑49 (MQ=255)
                                                                      atCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCg                       >  1:1663570/1‑49 (MQ=255)
                                                                      atCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCg                       >  1:5220306/1‑49 (MQ=255)
                                                                      atCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCTCCa  <  1:2876383/70‑1 (MQ=255)
                                                                      |                                                                     
GGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCA  >  Exported/13524‑13663

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: