Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 385,897 A→G 18.3% L369L (TTA→TTG) ‡ cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported385,8970AG18.3% 220.5 / 42.7 109L369L (TTA→TTG) ‡cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base A (40/49);  new base G (9/11);  total (49/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

Only 100 of 127 total aligned reads displayed.

TTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCA  >  Exported/385828‑385967
                                                                     |                                                                      
ttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTaa                                                                       >  1:25075/1‑71 (MQ=255)
  ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTa                                                                        <  1:4011975/68‑1 (MQ=255)
  ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTa                                                                        <  1:4899299/68‑1 (MQ=255)
  ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTa                                                                        <  1:2760664/68‑1 (MQ=255)
  ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTa                                                                        <  1:3034470/68‑1 (MQ=255)
  ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTa                                                                        <  1:2013463/68‑1 (MQ=255)
  ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAAt                                                                      >  1:1746415/1‑70 (MQ=255)
     ggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATgg                                                                    <  1:5649198/69‑1 (MQ=255)
     ggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATgg                                                                    <  1:2247715/69‑1 (MQ=255)
      gACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAAt                                                                 <  1:433986/71‑1 (MQ=255)
        ccACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATcc                                                               <  1:4484983/71‑1 (MQ=255)
        ccACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATcc                                                               <  1:2874289/71‑1 (MQ=255)
          aCGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCtt                                                             <  1:375380/71‑1 (MQ=255)
           caccGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTGATGGAATCCtt                                                             <  1:1252535/68‑1 (MQ=255)
           caccGAAGATGCTGGCGAACGAAGGGCAGACGTGCTTTATCGGCTACGGTGGGATGTTGATGGAATCCtt                                                             <  1:5587177/68‑1 (MQ=255)
           cGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTc                                                            <  1:1682076/71‑1 (MQ=255)
            gCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCg                                                           <  1:898131/71‑1 (MQ=255)
            gCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCg                                                           >  1:5626112/1‑71 (MQ=255)
             ccGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATGCTTCGt                                                          <  1:6528/71‑1 (MQ=255)
             ccGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGt                                                          <  1:4641613/71‑1 (MQ=255)
               gAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATcc                                                               >  1:4508440/1‑64 (MQ=255)
                    tGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTaa                                                                       >  1:4227598/1‑51 (MQ=255)
                    tGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTaa                                                                       >  1:5374554/1‑51 (MQ=255)
                    tGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTaa                                                                       >  1:3160814/1‑51 (MQ=255)
                    tGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTaa                                                                       >  1:3732506/1‑51 (MQ=255)
                     gCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATcc                                                               >  1:3680521/1‑58 (MQ=255)
                     gCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTa                                                  >  1:3148184/1‑71 (MQ=255)
                     gCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTa                                                  >  1:2766293/1‑71 (MQ=255)
                     gCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTa                                                  >  1:1189088/1‑71 (MQ=255)
                     gCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTa                                                  >  1:2620599/1‑71 (MQ=255)
                     gCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTa                                                  >  1:2123505/1‑71 (MQ=255)
                        ggCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTaa                                                                       <  1:3026065/47‑1 (MQ=255)
                         gCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGgc                                              >  1:789094/1‑71 (MQ=255)
                            aacgaaGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCt                                           <  1:2447684/71‑1 (MQ=255)
                              cgaaGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAAt                                                                 >  1:4405294/1‑47 (MQ=255)
                              cgaaGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAAt                                                                 >  1:4450174/1‑47 (MQ=255)
                              cgaaGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATc                                                                >  1:3750771/1‑48 (MQ=255)
                               gaaGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGa                                                     >  1:5062498/1‑58 (MQ=255)
                               gaaGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGa                                                     >  1:620664/1‑58 (MQ=255)
                               gaaGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGt                                        <  1:1472331/71‑1 (MQ=255)
                               gaaGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGt                                        <  1:4094316/71‑1 (MQ=255)
                                aaGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGtt                                       <  1:4389761/71‑1 (MQ=255)
                                   ggCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTcc                                    <  1:5588546/71‑1 (MQ=255)
                                   ggCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTcc                                    <  1:1768261/71‑1 (MQ=255)
                                     cAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGc                                  >  1:4717433/1‑71 (MQ=255)
                                        gCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTg                               <  1:5210306/71‑1 (MQ=255)
                                         agTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGt                              <  1:171864/70‑1 (MQ=255)
                                           tGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCtt                                                             <  1:5316669/38‑1 (MQ=255)
                                           tGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCtt                                                             <  1:3634666/38‑1 (MQ=255)
                                           tGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCtt                                                             <  1:4061539/38‑1 (MQ=255)
                                           tGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTg                                                         <  1:63155/42‑1 (MQ=255)
                                             cTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGa                                                     >  1:1709021/1‑44 (MQ=255)
                                             cTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatca                          <  1:2988911/71‑1 (MQ=255)
                                             cTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatca                          >  1:796552/1‑71 (MQ=255)
                                             cTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatc                           <  1:382975/70‑1 (MQ=255)
                                               tcaTCGGCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTa                                                  >  1:4385773/3‑45 (MQ=255)
                                               tcaTCGGCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCatcatc                        >  1:4048213/3‑71 (MQ=255)
                                               tcaTCGGCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCatcatc                        >  1:5459176/3‑71 (MQ=255)
                                               tcaTCGGCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCatcatc                        >  1:2478888/3‑71 (MQ=255)
                                               tcaTCGGCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCatcatc                        >  1:1024988/3‑71 (MQ=255)
                                               tcaTCGGCTACGGAGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCatcatc                        >  1:1412700/3‑71 (MQ=255)
                                                tATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGAt                                                    <  1:714986/42‑1 (MQ=255)
                                                tATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGAt                                                    <  1:4782108/42‑1 (MQ=255)
                                                tATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTcc                                    >  1:5249769/1‑58 (MQ=255)
                                                caTCGGCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCatcatc                        >  1:2165164/2‑70 (MQ=255)
                                                 aTCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTa                                                  >  1:2587767/1‑43 (MQ=255)
                                                 aTCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTa                                                  >  1:400964/1‑43 (MQ=255)
                                                 aTCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCg                       <  1:2290757/70‑1 (MQ=255)
                                                 aTCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCg                       <  1:2067750/70‑1 (MQ=255)
                                                  tCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTc                                     >  1:2316601/1‑55 (MQ=255)
                                                  tCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGa                      <  1:3811882/70‑1 (MQ=255)
                                                  tCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGa                      <  1:564175/70‑1 (MQ=255)
                                                  tCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGa                      <  1:342877/70‑1 (MQ=255)
                                                   cGGCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGAt                                                    <  1:5372381/39‑1 (MQ=12)
                                                    ggCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatcatc                        >  1:1404056/1‑66 (MQ=255)
                                                    ggCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatcatc                        >  1:4375191/1‑66 (MQ=255)
                                                    ggCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatcatc                        >  1:3924332/1‑66 (MQ=255)
                                                    ggCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGaac                    <  1:401342/70‑3 (MQ=255)
                                                    ggCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAAcc                   >  1:3009516/1‑71 (MQ=255)
                                                         cGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCggg                <  1:1775997/69‑1 (MQ=255)
                                                            tGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTg                               <  1:987350/51‑1 (MQ=255)
                                                             gggATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTa          >  1:4061144/1‑71 (MQ=255)
                                                               gATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatcatc                        <  1:2527693/55‑1 (MQ=255)
                                                               gATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAtt        <  1:3964663/71‑1 (MQ=255)
                                                                aTGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAttt       >  1:1417581/1‑71 (MQ=255)
                                                                 tGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGt                                        <  1:1261406/37‑1 (MQ=255)
                                                                  gTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTc                                     <  1:987297/39‑1 (MQ=255)
                                                                  gTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTc                                     <  1:560446/39‑1 (MQ=255)
                                                                  gTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTg     >  1:1918695/1‑71 (MQ=255)
                                                                   ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGa                      <  1:229017/53‑1 (MQ=255)
                                                                   ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGa                      <  1:2767376/53‑1 (MQ=255)
                                                                   ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGa                      <  1:782660/53‑1 (MQ=255)
                                                                   ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGAt                     >  1:3080136/1‑54 (MQ=255)
                                                                   ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGc    >  1:3382033/1‑71 (MQ=255)
                                                                    tAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatc                           >  1:536425/1‑47 (MQ=255)
                                                                     aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAttt       >  1:957108/1‑66 (MQ=255)
                                                                     aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGc    >  1:3199603/1‑69 (MQ=255)
                                                                     aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCa  >  1:4313481/1‑71 (MQ=255)
                                                                     aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCa  >  1:2477081/1‑71 (MQ=255)
                                                                     aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCa  >  1:442668/1‑71 (MQ=255)
                                                                     |                                                                      
TTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCA  >  Exported/385828‑385967

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: