Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 980,233 T→A 7.7% K24N (AAA→AAT) ‡ cspI ← cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported980,2330TA 7.7% 266.7 / 10.7 105K24N (AAA→AAT) ‡cspIcold shock protein
Reads supporting (aligned to +/- strand):  ref base T (41/52);  new base A (5/2);  total (48/57)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 116 total aligned reads displayed.

AGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACA  >  Exported/980167‑980303
                                                                  |                                                                      
aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCttt                                                                      <  1:3185473/69‑1 (MQ=255)
aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCttt                                                                      <  1:785838/69‑1 (MQ=255)
aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCtt                                                                       <  1:4769465/68‑1 (MQ=255)
aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCtt                                                                       <  1:1350852/68‑1 (MQ=255)
aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCtt                                                                       <  1:425608/68‑1 (MQ=255)
  ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                  <  1:1635304/71‑1 (MQ=255)
    aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCt                                                                        >  1:1827069/1‑63 (MQ=255)
    aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCt                                                                        >  1:1162587/1‑63 (MQ=255)
    aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                >  1:2682005/1‑71 (MQ=255)
     aTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                >  1:2279740/1‑70 (MQ=255)
     aTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                <  1:4822840/70‑1 (MQ=255)
        ttttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                                             <  1:5266302/70‑1 (MQ=255)
        ttttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                                             <  1:2931050/70‑1 (MQ=255)
        ttttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                                             <  1:805262/70‑1 (MQ=255)
        ttttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                  <  1:1655403/65‑1 (MQ=255)
          ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                           <  1:4530798/70‑1 (MQ=255)
          ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                           <  1:2261872/70‑1 (MQ=255)
               aTCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                                                      >  1:4960457/1‑70 (MQ=255)
               aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                     <  1:4837128/71‑1 (MQ=255)
                    tGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCttt                                                                      >  1:5331891/1‑49 (MQ=255)
                     gCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCt                                                <  1:3456392/70‑1 (MQ=255)
                      ctctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                                              >  1:4420352/1‑71 (MQ=255)
                       tctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                <  1:2962676/52‑1 (MQ=255)
                        ctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                                              >  1:767543/1‑53 (MQ=255)
                        ctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                                              >  1:2937930/1‑53 (MQ=255)
                        ctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                                              >  1:4610993/1‑53 (MQ=255)
                        ctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                            >  1:104868/1‑71 (MQ=255)
                         tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTc                                                                     <  1:3741904/45‑1 (MQ=255)
                         tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTc                                                                     <  1:3852159/45‑1 (MQ=255)
                         tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                            >  1:128643/1‑70 (MQ=255)
                         tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCATGATGAAACCaaa                                                    <  1:5561120/62‑1 (MQ=255)
                         tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                <  1:5022575/50‑1 (MQ=255)
                         tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                <  1:1798200/50‑1 (MQ=255)
                         tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                <  1:1197165/50‑1 (MQ=255)
                          gAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                    >  1:4012266/1‑61 (MQ=255)
                           aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                            <  1:1508089/68‑1 (MQ=255)
                           aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                            <  1:3675894/68‑1 (MQ=255)
                           aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                           <  1:2949205/69‑1 (MQ=255)
                             ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                          <  1:4578964/68‑1 (MQ=255)
                             ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                          <  1:4550533/68‑1 (MQ=255)
                             ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                          <  1:2919994/68‑1 (MQ=255)
                             ttGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                          <  1:5582488/68‑1 (MQ=255)
                               gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAc                                                  >  1:5680592/1‑58 (MQ=255)
                               gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGt                                      >  1:3022718/1‑70 (MQ=255)
                               gCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGt                                      >  1:3491141/1‑70 (MQ=255)
                                cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                            <  1:380875/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)
                                    gaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                <  1:4193400/39‑1 (MQ=255)
                                     aaaTGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCa                               <  1:2167906/71‑1 (MQ=255)
                                     aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCg                                                                 >  1:4403533/1‑37 (MQ=255)
                                     aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                 >  1:220450/1‑69 (MQ=255)
                                       atggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                             >  1:1039078/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)
                                           aCAAACACATCTTTGCTGCCATCTTTCGGCGtgatg                                                            >  1:286824/1‑36 (MQ=255)
                                           aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                    >  1:3443239/1‑44 (MQ=255)
                                           aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                           <  1:4785105/53‑1 (MQ=255)
                                           aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                           <  1:4048922/53‑1 (MQ=255)
                                           aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                          <  1:395464/54‑1 (MQ=255)
                                           aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                          <  1:4158835/54‑1 (MQ=255)
                                             aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                        >  1:4306035/1‑54 (MQ=255)
                                             aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                       <  1:5510559/71‑1 (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: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: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)
                                                      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)
                                                          gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt           <  1:3892167/68‑1 (MQ=255)
                                                          gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt           <  1:2003423/68‑1 (MQ=255)
                                                          gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAttt           <  1:82126/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:3092487/53‑1 (MQ=255)
                                                               aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                       <  1:160521/53‑1 (MQ=255)
                                                                  tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTa      >  1:5261761/1‑67 (MQ=255)
                                                                  tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa  >  1:1425652/1‑71 (MQ=255)
                                                                  |                                                                      
AGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACA  >  Exported/980167‑980303

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: