Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,211,333 G→T 12.8% Q7K (CAG→AAG) ‡ cspC ← stress protein, member of the CspA‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,211,3330GT12.8% 279.1 / 16.3 93Q7K (CAG→AAG) ‡cspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base G (35/46);  new base T (4/8);  total (39/54)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGA  >  Exported/1211264‑1211399
                                                                     |                                                                  
acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                                                                   <  1:3971008/71‑1 (MQ=255)
acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                                                                   <  1:3624785/71‑1 (MQ=255)
acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                                                                   >  1:2078914/1‑71 (MQ=255)
   cgtCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCt                                                                >  1:3944241/3‑71 (MQ=12)
   caTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCt                                                                <  1:1152790/71‑1 (MQ=255)
   caTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCt                                                                <  1:4781149/71‑1 (MQ=255)
    aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCtt                                                               <  1:4970200/71‑1 (MQ=255)
       tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTg                                                                    >  1:5654221/1‑63 (MQ=255)
       tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTaa                                                            >  1:4453267/1‑71 (MQ=255)
         tgctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTaa                                                            <  1:3695205/69‑1 (MQ=255)
         tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                             <  1:3454643/68‑1 (MQ=255)
         tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                             <  1:3678572/68‑1 (MQ=255)
         tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                             <  1:2164328/68‑1 (MQ=255)
          gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCt                                                         >  1:3595011/1‑71 (MQ=255)
          gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCt                                                         >  1:5137044/1‑71 (MQ=255)
           ctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCtt                                                        <  1:5041135/71‑1 (MQ=255)
           ctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCt                                                         <  1:3743808/70‑1 (MQ=255)
            tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                             <  1:402581/65‑1 (MQ=255)
            tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                             <  1:1711364/65‑1 (MQ=255)
              ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGc                                                     >  1:1797797/1‑71 (MQ=255)
                aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa                                                   <  1:5384946/71‑1 (MQ=255)
                 tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCt                                                         <  1:3575/64‑1 (MQ=255)
                 tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCt                                                         <  1:2535722/64‑1 (MQ=255)
                 tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                                                  <  1:1346808/71‑1 (MQ=255)
                 tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                                                  >  1:5457350/1‑71 (MQ=255)
                   aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCttt                                                       >  1:4610099/1‑64 (MQ=255)
                   aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt                                                 >  1:3981886/1‑70 (MQ=255)
                     ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCTTAGACAtcaa                                               <  1:4366361/70‑4 (MQ=255)
                     ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCTTAGACAtcaa                                               <  1:3496766/70‑4 (MQ=255)
                     ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCTTAGACAtcaa                                               <  1:3476175/70‑4 (MQ=255)
                     ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCTTAGACAtca                                                <  1:3627065/69‑3 (MQ=255)
                     ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCTTAGACAtca                                                <  1:3136418/69‑3 (MQ=255)
                       ggAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa                                                   <  1:2943115/64‑1 (MQ=255)
                       ggAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa                                            <  1:3067404/71‑1 (MQ=255)
                       ggAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa                                            <  1:3113148/71‑1 (MQ=255)
                         aGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt                                                 >  1:3768443/1‑64 (MQ=255)
                         aGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCttag                                                      >  1:986/1‑57 (MQ=16)
                         aGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCtta                                                       >  1:810974/1‑57 (MQ=16)
                          gTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTAcc                                                                 >  1:3366984/1‑47 (MQ=9)
                           tAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCtt                                                               <  1:4339395/48‑1 (MQ=255)
                           tAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCtt                                                               <  1:1822953/48‑1 (MQ=255)
                           tAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa                                            <  1:4087591/67‑1 (MQ=255)
                            aaTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCttt                                                              <  1:3264398/48‑1 (MQ=255)
                             aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTaa                                                            >  1:172183/1‑49 (MQ=255)
                             aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTaa                                                            >  1:4980506/1‑49 (MQ=255)
                             aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                      >  1:3116939/1‑71 (MQ=255)
                             aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                      >  1:3877377/1‑71 (MQ=255)
                             aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                      <  1:2504374/71‑1 (MQ=255)
                             aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                      >  1:3092134/1‑71 (MQ=255)
                             aTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                      <  1:3240707/71‑1 (MQ=255)
                              tGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                                                  <  1:4771719/58‑1 (MQ=255)
                              tGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                                                  <  1:5568382/58‑1 (MQ=255)
                              tGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCt                                     <  1:541882/71‑1 (MQ=255)
                                aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                                          >  1:4413746/1‑64 (MQ=255)
                                aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                                          >  1:1909686/1‑64 (MQ=255)
                                aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                                          >  1:2473115/1‑64 (MQ=255)
                                aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                                          >  1:2466670/1‑64 (MQ=255)
                                aaGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa                                          >  1:560782/1‑64 (MQ=255)
                                  gCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCt                                     >  1:2943908/1‑67 (MQ=255)
                                   ccAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTTcc                                                    <  1:2447139/51‑1 (MQ=255)
                                    cAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGAt                                >  1:2456088/1‑70 (MQ=255)
                                     aaaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATc                                                          <  1:2119062/43‑1 (MQ=255)
                                     aaaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGcc                                                    <  1:3424188/49‑1 (MQ=255)
                                     aaaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                                                  >  1:2056603/1‑51 (MQ=255)
                                     aaaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTg                              <  1:3300760/71‑1 (MQ=255)
                                      aaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                             <  1:56089/39‑1 (MQ=255)
                                      aaaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                             <  1:4757809/39‑1 (MQ=255)
                                       aaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt                            >  1:292420/1‑71 (MQ=255)
                                       aaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt                            <  1:2991673/71‑1 (MQ=255)
                                       aaCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt                            <  1:4897250/71‑1 (MQ=255)
                                        aCCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAAt                                                           >  1:3057925/1‑39 (MQ=255)
                                          cTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTa                                   <  1:4466759/61‑1 (MQ=255)
                                          cTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTa                                   <  1:455692/61‑1 (MQ=255)
                                          cTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttc                         <  1:1593635/71‑1 (MQ=255)
                                          cTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttc                         <  1:1600336/71‑1 (MQ=255)
                                           tttAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttct                        >  1:492893/1‑71 (MQ=255)
                                            ttAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttc                         >  1:4366700/1‑69 (MQ=255)
                                                    gTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCg                 >  1:5641283/1‑69 (MQ=255)
                                                      tGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg             >  1:746258/1‑71 (MQ=255)
                                                       gAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg             >  1:726041/1‑70 (MQ=255)
                                                       gAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg             >  1:3792355/1‑70 (MQ=255)
                                                       gAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg             <  1:4010527/70‑1 (MQ=255)
                                                           cACTTAACCTGACCTTTAATCTTTGCCATTTGaaa                                            <  1:2944527/35‑1 (MQ=255)
                                                           cACTTAACCTGACCTTTAATCTTTGCCATTTGaaa                                            <  1:2895290/35‑1 (MQ=255)
                                                             cTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt                            <  1:5642193/49‑1 (MQ=255)
                                                             cTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt                            <  1:2023398/49‑1 (MQ=255)
                                                             cTTAACATGACCGTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACa      >  1:6902/1‑71 (MQ=255)
                                                               tAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCg                 >  1:2381631/1‑58 (MQ=255)
                                                                aaCCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                      >  1:4456357/1‑36 (MQ=255)
                                                                aaCCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtttt                          >  1:2220450/1‑48 (MQ=255)
                                                                 aCCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAg   <  1:2264444/70‑1 (MQ=255)
                                                                 aCCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGa  <  1:3641857/71‑1 (MQ=255)
                                                                    tGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtttt                          >  1:5311271/1‑44 (MQ=255)
                                                                    tGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGcc                   <  1:83661/51‑1 (MQ=255)
                                                                     |                                                                  
ACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGA  >  Exported/1211264‑1211399

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: