Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 980,245 G→A 7.9% F20F (TTC→TTT)  cspI ← cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported980,2450GA 7.9% 252.5 / 8.3 88F20F (TTC→TTT) cspIcold shock protein
Reads supporting (aligned to +/- strand):  ref base G (44/37);  new base A (5/2);  total (49/39)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 108 total aligned reads displayed.

TTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTA  >  Exported/980177‑980309
                                                                    |                                                                
ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                                 <  1:2261872/70‑1 (MQ=255)
ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                                 <  1:4530798/70‑1 (MQ=255)
    aaTCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                                                            <  1:5691882/71‑1 (MQ=255)
     aTCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                                                            >  1:4960457/1‑70 (MQ=255)
     aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                           <  1:4837128/71‑1 (MQ=255)
           gCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCt                                                      <  1:3456392/70‑1 (MQ=255)
            ctctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                                                    >  1:4420352/1‑71 (MQ=255)
              ctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                  >  1:104868/1‑71 (MQ=255)
               tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                  >  1:1620991/1‑70 (MQ=255)
               tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                  >  1:128643/1‑70 (MQ=255)
               tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCATGATGAAACCaaa                                                          <  1:5561120/62‑1 (MQ=255)
                gAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                          >  1:4012266/1‑61 (MQ=255)
                 aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                  <  1:3675894/68‑1 (MQ=255)
                 aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                  <  1:1508089/68‑1 (MQ=255)
                 aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                                 <  1:2949205/69‑1 (MQ=255)
                 aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                               <  1:4446119/71‑1 (MQ=255)
                   ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatg                                                                  >  1:2411473/1‑50 (MQ=255)
                   ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                <  1:2919994/68‑1 (MQ=255)
                   ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                <  1:4550533/68‑1 (MQ=255)
                   ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                <  1:5582488/68‑1 (MQ=255)
                   ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                <  1:4578964/68‑1 (MQ=255)
                     gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAc                                                        >  1:5680592/1‑58 (MQ=255)
                     gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGt                                            >  1:3491141/1‑70 (MQ=255)
                     gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGt                                            >  1:3022718/1‑70 (MQ=255)
                      cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                  <  1:380875/63‑1 (MQ=255)
                      cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                  <  1:5457685/63‑1 (MQ=255)
                      cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa                                          <  1:3922659/71‑1 (MQ=255)
                      cTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                >  1:2389092/1‑65 (MQ=255)
                       tGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                  <  1:3410184/62‑1 (MQ=255)
                        gagaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaaa                                        >  1:1216871/1‑71 (MQ=255)
                          gaAATGGACAAACTCATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAc                                                        >  1:2788336/1‑53 (MQ=37)
                           aaaTGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCa                                     <  1:2167906/71‑1 (MQ=255)
                           aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                       >  1:220450/1‑69 (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)
                                caCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCAtttta                                >  1:1453682/2‑69 (MQ=255)
                                caCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACAAtttta                                >  1:1298922/2‑69 (MQ=255)
                                 aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                          >  1:4190745/1‑44 (MQ=255)
                                 aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                          >  1:3443239/1‑44 (MQ=255)
                                 aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                                 <  1:4048922/53‑1 (MQ=255)
                                 aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                                 <  1:4785105/53‑1 (MQ=255)
                                 aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                <  1:395464/54‑1 (MQ=255)
                                   aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                              >  1:4306035/1‑54 (MQ=255)
                                   aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                             <  1:2315796/71‑1 (MQ=255)
                                   aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                             <  1:720516/71‑1 (MQ=255)
                                   aaaCACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCATTTTACt                              >  1:5153911/1‑70 (MQ=255)
                                   aaaCACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCATTTTACt                              >  1:5127617/1‑70 (MQ=255)
                                     acacATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTAcga                             >  1:4621677/1‑67 (MQ=9)
                                     acacATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTAcga                             >  1:1052620/1‑67 (MQ=9)
                                      cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                                  >  1:1567424/1‑63 (MQ=255)
                                      cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTc                                 >  1:1156252/1‑64 (MQ=255)
                                        caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                               >  1:631642/1‑48 (MQ=255)
                                        caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                               >  1:1659681/1‑48 (MQ=255)
                                        caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                               >  1:5185732/1‑48 (MQ=255)
                                           cTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCa                        >  1:5511648/1‑68 (MQ=255)
                                            tttGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                                                       >  1:1587214/1‑36 (MQ=255)
                                            tttGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                                                       >  1:1263479/1‑36 (MQ=255)
                                            tttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                       <  1:5616151/52‑1 (MQ=255)
                                            tttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCa                    >  1:3707685/1‑71 (MQ=255)
                                             ttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGt                      >  1:3644293/1‑68 (MQ=255)
                                             ttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAt                   >  1:5221626/1‑71 (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:82126/68‑1 (MQ=255)
                                                gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt                 <  1:3892167/68‑1 (MQ=255)
                                                gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt                 <  1:2003423/68‑1 (MQ=255)
                                                ctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAttt                 <  1:1591862/70‑1 (MQ=255)
                                                  gcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                       <  1:3448645/46‑1 (MQ=255)
                                                   ccATCTTTGGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt              >  1:2737445/1‑70 (MQ=255)
                                                   ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt              <  1:4351269/70‑1 (MQ=255)
                                                     aTCTTTCGGTGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                             <  1:3328034/53‑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)
                                                         ttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa        >  1:4028751/1‑70 (MQ=255)
                                                         ttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa        >  1:2488971/1‑70 (MQ=255)
                                                          tCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt       >  1:2333624/1‑70 (MQ=255)
                                                           cGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATaa     >  1:151545/1‑71 (MQ=255)
                                                             gagTGATGAAGCCGAAGCCTTTGTCAGCGTTGAGCCATTTTACGATACCAGTCAtt                  >  1:3514997/3‑56 (MQ=255)
                                                             gagTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtt                  >  1:4568935/3‑56 (MQ=255)
                                                             gagTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtt                  >  1:1857598/3‑56 (MQ=255)
                                                               gTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatta  >  1:2103233/1‑70 (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:8492/55‑1 (MQ=255)
                                                                  atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt              <  1:4753505/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:4645415/1‑52 (MQ=255)
                                                                    gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg               >  1:2865127/1‑52 (MQ=255)
                                                                    gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatt   <  1:1845520/64‑1 (MQ=255)
                                                                    gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatt   <  1:4098546/64‑1 (MQ=255)
                                                                    |                                                                
TTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTA  >  Exported/980177‑980309

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: ‑