Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 980,260 T→A 15.2% E15D (GAA→GAT cspI ← cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported980,2600TA15.2% 242.0 / 29.7 113E15D (GAA→GATcspIcold shock protein
Reads supporting (aligned to +/- strand):  ref base T (49/41);  new base A (10/4);  total (65/48)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 133 total aligned reads displayed.

CTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAG  >  Exported/980191‑980329
                                                                     |                                                                     
ctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                      >  1:104868/1‑71 (MQ=255)
   aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                      <  1:3675894/68‑1 (MQ=255)
   aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                      <  1:1508089/68‑1 (MQ=255)
   aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                                                     <  1:2949205/69‑1 (MQ=255)
     ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                                    <  1:4550533/68‑1 (MQ=255)
     ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                                    <  1:4578964/68‑1 (MQ=255)
       gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGt                                                                >  1:3022718/1‑70 (MQ=255)
        cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                      <  1:380875/63‑1 (MQ=255)
        cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                      <  1:5457685/63‑1 (MQ=255)
        cTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                                    >  1:2389092/1‑65 (MQ=255)
         tGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                      <  1:3410184/62‑1 (MQ=255)
          gagaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaaa                                                            >  1:1216871/1‑71 (MQ=255)
             aaaTGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCa                                                         <  1:2167906/71‑1 (MQ=255)
               atggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                                       >  1:1039078/4‑71 (MQ=255)
               atggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                                       >  1:1490823/4‑71 (MQ=255)
                tggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                                                      <  1:1379447/69‑1 (MQ=255)
                  gACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                                   >  1:2381801/1‑56 (MQ=255)
                  caCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACAAtttta                                                    >  1:1298922/2‑69 (MQ=255)
                   aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                                                     <  1:4785105/53‑1 (MQ=255)
                   aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                                                     <  1:4048922/53‑1 (MQ=255)
                   aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                                    <  1:4158835/54‑1 (MQ=255)
                   aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                                    <  1:395464/54‑1 (MQ=255)
                     aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                                 <  1:720516/71‑1 (MQ=255)
                     aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                                 <  1:2315796/71‑1 (MQ=255)
                     aaaCACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCATTTTACt                                                  >  1:5153911/1‑70 (MQ=255)
                       acacATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTAcga                                                 >  1:1052620/1‑67 (MQ=9)
                       acacATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTAcga                                                 >  1:4621677/1‑67 (MQ=9)
                        cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                                                      >  1:1567424/1‑63 (MQ=255)
                        cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTc                                                     >  1:1156252/1‑64 (MQ=255)
                        cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                                              <  1:1513145/71‑1 (MQ=255)
                          caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                                   >  1:1659681/1‑48 (MQ=255)
                          caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                                   >  1:5185732/1‑48 (MQ=255)
                             cTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCa                                            >  1:5511648/1‑68 (MQ=255)
                              tttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                                           <  1:5616151/52‑1 (MQ=255)
                               ttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGt                                          >  1:3644293/1‑68 (MQ=255)
                                tgctgcCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAAc                                              >  1:4022803/1‑63 (MQ=9)
                                 gctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtt                                      <  1:4479250/70‑1 (MQ=255)
                                 gctgcCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCAtt                                      <  1:2100745/70‑1 (MQ=12)
                                  gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt                                     <  1:2003423/68‑1 (MQ=255)
                                  gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt                                     <  1:3892167/68‑1 (MQ=255)
                                  ctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAttt                                     <  1:1591862/70‑1 (MQ=255)
                                    gcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                                           <  1:3448645/46‑1 (MQ=255)
                                     ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                  <  1:4351269/70‑1 (MQ=255)
                                       aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                                 <  1:160521/53‑1 (MQ=255)
                                       aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                                 <  1:3092487/53‑1 (MQ=255)
                                          tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTa                                >  1:5261761/1‑67 (MQ=255)
                                          tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa                            >  1:1425652/1‑71 (MQ=255)
                                            tCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt                           >  1:2333624/1‑70 (MQ=255)
                                               gagTGATGAAGCCGAAGCCTTTGTCAGCGTTGAGCCATTTTACGATACCAGTCAtt                                      >  1:3514997/3‑56 (MQ=255)
                                               gagTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtt                                      >  1:1857598/3‑56 (MQ=255)
                                                  tgatgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCTCTAAACCAGTCATTTTGTTAGACATAatta                      <  1:4387484/69‑1 (MQ=255)
                                                  tgataaAGCCAAAACCTTTATCTGCGTTAAACCAt                                                        <  1:5295674/35‑1 (MQ=12)
                                                   gatgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatta                      <  1:3978299/68‑1 (MQ=255)
                                                    atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAttt                                     <  1:3203337/52‑1 (MQ=255)
                                                    atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                  <  1:4986770/55‑1 (MQ=255)
                                                    atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                  <  1:210327/55‑1 (MQ=255)
                                                    atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                  <  1:4753505/55‑1 (MQ=255)
                                                    atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                  <  1:8492/55‑1 (MQ=255)
                                                    atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                  <  1:5633836/55‑1 (MQ=255)
                                                     tgaAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtt                                      >  1:899380/1‑50 (MQ=255)
                                                     tgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat                        >  1:2128480/1‑64 (MQ=255)
                                                      gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg                                   >  1:2865127/1‑52 (MQ=255)
                                                      gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg                                   >  1:4645415/1‑52 (MQ=255)
                                                       aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAATCCAGTCATTTTGTTAGACATAATTATTAc                  >  1:3841039/1‑68 (MQ=255)
                                                       aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCtt               >  1:923052/1‑71 (MQ=255)
                                                        aaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCttt              >  1:4390881/1‑71 (MQ=255)
                                                          ccAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg            <  1:2256864/71‑1 (MQ=255)
                                                          ccAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTc                                         >  1:1237644/1‑42 (MQ=25)
                                                          ccAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCa                                        >  1:1923790/1‑43 (MQ=25)
                                                          ccAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCa                                        >  1:3179028/1‑43 (MQ=25)
                                                           cAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtttt                                    >  1:5024372/1‑46 (MQ=255)
                                                           cAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg            >  1:2594071/1‑70 (MQ=255)
                                                           cAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg            >  1:5017903/1‑70 (MQ=255)
                                                            aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                  >  1:5375720/1‑47 (MQ=255)
                                                            aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg            >  1:5498046/1‑69 (MQ=255)
                                                            aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg            >  1:2074788/1‑69 (MQ=255)
                                                            aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg            >  1:3169106/1‑69 (MQ=255)
                                                            aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACAAGTCATTTTGt                                  >  1:2076230/1‑47 (MQ=255)
                                                            aaaaCCTTTATCTGCGTTAAACCATTTTACTAAACCAGt                                          <  1:4281685/39‑1 (MQ=25)
                                                             aaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTa                                >  1:2629414/1‑48 (MQ=255)
                                                             aaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg            >  1:3407416/1‑68 (MQ=255)
                                                              aaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaaga        <  1:405913/71‑1 (MQ=255)
                                                              aaCCTTTATCAGCGTTAAACCATTTTACTAAACCa                                            >  1:1012348/1‑35 (MQ=25)
                                                               aCCTTTTTCAGGGTTAAACCATTTCACTATACCAGTCATTTTGt                                  <  1:5152605/44‑1 (MQ=255)
                                                               aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAattat                     >  1:4194263/1‑57 (MQ=255)
                                                               aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAattat                     >  1:2418489/1‑57 (MQ=255)
                                                               aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaagaa       <  1:5581304/71‑1 (MQ=255)
                                                               aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaagaa       <  1:2837519/71‑1 (MQ=255)
                                                                ccTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGaaa      >  1:691397/1‑71 (MQ=255)
                                                                 cTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAAt     <  1:2764800/71‑1 (MQ=255)
                                                                  tttgtCAGCGTTGAACCATTTTACGATACCAGTCAttt                                     <  1:4874612/38‑1 (MQ=12)
                                                                  tttgtCAGCGTTGAACCATTTTACGATACCAGTCAttt                                     <  1:5566213/38‑1 (MQ=12)
                                                                  tttatCAGCGTTAAACCATTTTACTAAACGAGTCAttt                                     <  1:447067/38‑1 (MQ=21)
                                                                  tttatCAGCGTTAAACCATTTTACTAAACCAGTCAttt                                     <  1:5506022/38‑1 (MQ=25)
                                                                  tttatCAGCGTTAAACCATTTTACTAAACCAGTCATTTTATTTGa                              <  1:1796779/45‑1 (MQ=255)
                                                                    tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAAtt    <  1:2627508/69‑1 (MQ=255)
                                                                    tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAAtt    <  1:4040715/69‑1 (MQ=255)
                                                                    tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTa   <  1:5659966/70‑1 (MQ=255)
                                                                    tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTa   <  1:165492/70‑1 (MQ=255)
                                                                    tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAg  >  1:697425/1‑71 (MQ=255)
                                                                     |                                                                     
CTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAG  >  Exported/980191‑980329

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: