Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 980,296 G→A 7.5% N3N (AAC→AAT cspI ← cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported980,2960GA 7.5% 311.8 / 10.0 106N3N (AAC→AATcspIcold shock protein
Reads supporting (aligned to +/- strand):  ref base G (53/45);  new base A (3/5);  total (56/50)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 108 total aligned reads displayed.

CCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCA  >  Exported/980228‑980362
                                                                    |                                                                  
ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                                                   <  1:4351269/70‑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)
            gTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatta                                                       >  1:2103233/1‑70 (MQ=255)
             tgatgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCTCTAAACCAGTCATTTTGTTAGACATAatta                                                       <  1:4387484/69‑1 (MQ=255)
              gatgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatta                                                       <  1:3978299/68‑1 (MQ=255)
               atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                                                   <  1:4986770/55‑1 (MQ=255)
               atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                                                   <  1:5633836/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)
                 gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg                                                                    >  1:4645415/1‑52 (MQ=255)
                 gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg                                                                    >  1:2865127/1‑52 (MQ=255)
                 gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatt                                                        <  1:4098546/64‑1 (MQ=255)
                 gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatt                                                        <  1:1845520/64‑1 (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)
                      cAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                                             >  1:5017903/1‑70 (MQ=255)
                      cAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                                             >  1:2594071/1‑70 (MQ=255)
                       aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                                                   >  1:5375720/1‑47 (MQ=255)
                       aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                                             >  1:2074788/1‑69 (MQ=255)
                       aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                                             >  1:5498046/1‑69 (MQ=255)
                       aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                                             >  1:5448399/1‑69 (MQ=255)
                       aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                                             >  1:3169106/1‑69 (MQ=255)
                       aaaaCCTTTTTCAGGGTTAAACCATTTCACTAAACAAGTCATTTTGt                                                                   >  1:2076230/1‑47 (MQ=255)
                        aaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTa                                                                 >  1:2629414/1‑48 (MQ=255)
                        aaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                                             >  1:3407416/1‑68 (MQ=255)
                        aaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                                             >  1:3166722/1‑68 (MQ=255)
                         aaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaaga                                         <  1:405913/71‑1 (MQ=255)
                          aCCTTTTTCAGGGTTAAACCATTTCACTATACCAGTCATTTTGt                                                                   <  1:5152605/44‑1 (MQ=255)
                          aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAattat                                                      >  1:4194263/1‑57 (MQ=255)
                          aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAattat                                                      >  1:4752905/1‑57 (MQ=255)
                          aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAattat                                                      >  1:2418489/1‑57 (MQ=255)
                          aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaagaa                                        <  1:2837519/71‑1 (MQ=255)
                          aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaagaa                                        <  1:5581304/71‑1 (MQ=255)
                           ccTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGaaa                                       >  1:691397/1‑71 (MQ=255)
                            cTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAAt                                      <  1:2764800/71‑1 (MQ=255)
                             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)
                               tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAg                                   >  1:697425/1‑71 (MQ=255)
                                 tCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGc                                  <  1:2869713/70‑1 (MQ=255)
                                 tCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGc                                  <  1:1376651/70‑1 (MQ=255)
                                 tCAGCGTTAAACCATTTTACTAAACCAGTCATTTTAttt                                                                 <  1:3098439/39‑2 (MQ=25)
                                  cAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGc                                  <  1:778949/69‑1 (MQ=255)
                                   aGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCt                               <  1:1923509/71‑1 (MQ=255)
                                       ttAAACCATTTTACTAAACCAGTCATTTTATTTGa                                                               <  1:5405127/35‑1 (MQ=25)
                                       ttAAACCATTTTACTAAACCAGTCATTTTATTAGACAt                                                            <  1:3734753/38‑1 (MQ=38)
                                       ttAAACCATTTTACTAAACCAGTCATTTTATTAGACAt                                                            <  1:26667/38‑1 (MQ=38)
                                       ttAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat                                                         <  1:2705755/41‑1 (MQ=255)
                                        tAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAAt                                      <  1:5585337/59‑1 (MQ=255)
                                         aaaCCATTTTACTAAACCAGTCATTTTATTAGACATat                                                          >  1:208707/1‑37 (MQ=37)
                                         aaaCCATTTTACTAAACCAGTCATTTTATTAGACATAtttta                                                      >  1:2586244/1‑40 (MQ=37)
                                         aaaCCATTTTACTAAACCAGTCATTTTATTAGACATAtttta                                                      >  1:4194764/1‑40 (MQ=37)
                                          aaCCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAg                        <  1:822315/71‑1 (MQ=255)
                                          aaCCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAg                        <  1:1135570/71‑1 (MQ=255)
                                          aaCCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAg                        <  1:5210032/71‑1 (MQ=255)
                                           aCCATTTCACTAAACCAGTCATTTTGTTAGACATAatta                                                       >  1:1228825/1‑39 (MQ=255)
                                           aCCATTTCACTAAACCAGTCATTTTGTTAGACATAatta                                                       >  1:3855462/1‑39 (MQ=255)
                                           aCCATTTCACTAAACCAGTCATTTTGTTAGACATAatta                                                       >  1:3903845/1‑39 (MQ=255)
                                           aCCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTg                             <  1:5656010/65‑1 (MQ=255)
                                               tttCACTAAACCAGTCATTTTGTTAGACATAatta                                                       >  1:1697838/1‑35 (MQ=255)
                                               tttCACTAAACCAGTCATTTTGTTAGACATAatta                                                       >  1:1317535/1‑35 (MQ=255)
                                               tttCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaagaa                                        <  1:805918/50‑1 (MQ=255)
                                               tttCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATgg                   <  1:324009/71‑1 (MQ=255)
                                                 tCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGaaa                                       >  1:1662990/1‑49 (MQ=255)
                                                  cACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTcc                >  1:2260131/1‑71 (MQ=255)
                                                    cTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTg                             <  1:1574227/56‑1 (MQ=255)
                                                    cTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGa              <  1:1615037/71‑1 (MQ=255)
                                                     tAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGa              >  1:4667678/1‑70 (MQ=255)
                                                      aaaCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTa                                    <  1:4349363/47‑1 (MQ=255)
                                                      aaaCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGaa             >  1:1268701/1‑70 (MQ=255)
                                                       aaCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGaaa                                       >  1:4647215/1‑43 (MQ=255)
                                                       aaCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGaaa                                       >  1:5162719/1‑43 (MQ=255)
                                                       aaCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGaaa                                       >  1:3400259/1‑43 (MQ=255)
                                                       aaCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGaaa            <  1:2944775/70‑1 (MQ=255)
                                                            gTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCtt                              <  1:2945207/47‑1 (MQ=255)
                                                            gTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAat      <  1:5500449/71‑1 (MQ=255)
                                                             tCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGc                                  <  1:1190865/42‑1 (MQ=255)
                                                             tCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTgg                            <  1:3166608/48‑1 (MQ=255)
                                                             tCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTgg                            <  1:908996/48‑1 (MQ=255)
                                                             tCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATg                    <  1:2905249/56‑1 (MQ=255)
                                                             tCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAata     <  1:1167003/71‑1 (MQ=255)
                                                               aTTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCt                               >  1:2719502/1‑43 (MQ=255)
                                                               aTTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAatat    >  1:1868999/1‑70 (MQ=255)
                                                               aTTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAatat    >  1:5624157/1‑70 (MQ=255)
                                                               aTTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAatat    >  1:2534053/1‑70 (MQ=255)
                                                               aTTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAatat    <  1:4326628/70‑1 (MQ=255)
                                                                 tttGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCt                               >  1:2000319/1‑41 (MQ=255)
                                                                 tttGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCt                               >  1:1543214/1‑41 (MQ=255)
                                                                  ttGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCa  >  1:3122699/1‑69 (MQ=255)
                                                                   tGTTAGACATAATTATTACCTTTTGAAGAAATTAGccc                                >  1:13992/1‑38 (MQ=255)
                                                                   tGTTAGACATAATTATTACCTTTTGAAGAAATTAGccc                                >  1:4148454/1‑38 (MQ=255)
                                                                    |                                                                  
CCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCA  >  Exported/980228‑980362

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: