Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 385,847 A→C 23.8% M353L (ATG→CTG)  cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported385,8470AC23.8% 138.7 / 67.8 101M353L (ATG→CTG) cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base A (32/45);  new base C (14/10);  total (46/55)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

Only 100 of 109 total aligned reads displayed.

TTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGT  >  Exported/385778‑385911
                                                                     |                                                                
tttATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                                                                  <  1:5116961/70‑1 (MQ=255)
tttATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                                                                  <  1:1576820/70‑1 (MQ=255)
tttATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                                                                  <  1:1661703/70‑1 (MQ=255)
      aCCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTgg                                                            >  1:3700686/1‑70 (MQ=255)
      aCCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGc                                                           >  1:4549970/1‑71 (MQ=255)
      aCCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGc                                                           >  1:2658140/1‑71 (MQ=255)
       ccATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGct                                                          <  1:3375521/71‑2 (MQ=255)
          tCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTgg                                                            <  1:4625930/66‑1 (MQ=255)
          tCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaac                                                       >  1:5309445/1‑71 (MQ=255)
          tCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaac                                                       >  1:1952477/1‑71 (MQ=255)
           cGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaat                                                       >  1:1707696/1‑69 (MQ=255)
            gCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatg                                                      >  1:5027682/1‑68 (MQ=255)
            gCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatg                                                      <  1:5170441/70‑3 (MQ=255)
            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                                                      <  1:2238044/70‑1 (MQ=255)
            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                                                      <  1:3113432/70‑1 (MQ=255)
            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                                                      <  1:2474919/70‑1 (MQ=255)
            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                                                      <  1:4227084/70‑1 (MQ=255)
            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                                                      <  1:133727/70‑1 (MQ=255)
            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcg                                                          >  1:3418768/1‑66 (MQ=255)
             ccTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatga                                                     >  1:1528704/1‑70 (MQ=255)
             ccTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatga                                                     >  1:4632422/1‑70 (MQ=255)
              cTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAg                                                   >  1:2125439/1‑71 (MQ=255)
                 tgGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaaacc                                                 <  1:89464/70‑4 (MQ=255)
                 tgGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGct                                                          <  1:1943278/61‑2 (MQ=255)
                 tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                                                                  >  1:3948685/1‑53 (MQ=255)
                 tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                                                                  >  1:4199975/1‑53 (MQ=255)
                 tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                                                                  >  1:2359027/1‑53 (MQ=255)
                 tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                                                                  >  1:5612503/1‑53 (MQ=255)
                 tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                                                                  >  1:2217885/1‑53 (MQ=255)
                  ggCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATg                                                                >  1:5336243/1‑54 (MQ=255)
                  ggCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATg                                                                >  1:4360710/1‑54 (MQ=255)
                  ggCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATg                                                                >  1:2834309/1‑54 (MQ=255)
                  ggCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTAGGGACCACGCCGAAGATg                                                                >  1:2771463/1‑53 (MQ=255)
                   gtgcgGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaa                                                        >  1:36651/3‑61 (MQ=255)
                   gtgcgGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaa                                                        >  1:5412112/3‑61 (MQ=255)
                   gtgcgGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGc                                                           >  1:3204708/3‑58 (MQ=255)
                         tgtCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                                                      <  1:4439594/57‑1 (MQ=255)
                         tgtCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGc                                        <  1:510735/71‑1 (MQ=255)
                           tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                                                      <  1:1964989/55‑1 (MQ=255)
                           tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                                      >  1:4033926/1‑71 (MQ=255)
                           tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                                      <  1:3530569/71‑1 (MQ=255)
                           tCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaat                                                       <  1:256970/54‑2 (MQ=255)
                             tGGCTTCCCTGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTa                                    >  1:4829790/1‑71 (MQ=255)
                             tGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTa                                    >  1:383174/1‑71 (MQ=255)
                             tGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaaac                                                  >  1:3451681/1‑55 (MQ=255)
                               tcTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGAGAACGAAGGGCAGGCGTGATTTATc                                  >  1:1509315/2‑71 (MQ=255)
                               gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCttt                                     <  1:4386772/68‑1 (MQ=255)
                               gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCttt                                     <  1:2768365/68‑1 (MQ=255)
                               gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATc                                  >  1:889385/1‑71 (MQ=255)
                               gCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaa                                                        >  1:3400887/1‑49 (MQ=255)
                                    cATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTa                             <  1:385124/71‑1 (MQ=255)
                                    cATGCGCTGATCTCTTCTGGGACCACGCAGAAGATGCTGGCGAACGAAGGGCAGGc                                            >  1:4279267/1‑56 (MQ=255)
                                     acgcgcTGATCTCTTCTGGTACGACGCCAAAACTGCTGGCTAATGAAACCGACGCGCGTTTCATCGGCTAc                            <  1:5073468/69‑1 (MQ=255)
                                     acgcgcTGATCTCTTCCGGTACGACGCCAAAACTGCTGGCTAATGAAACCGACGCGCGTTTCATCGGCTAc                            <  1:3511640/69‑1 (MQ=255)
                                     acgcgcTGATCTCTTCCGGTACGACGCCAAAACTGCTGGCTAATGAAACCGACGCGCGTTTCATCGGCTAc                            <  1:3148854/69‑1 (MQ=255)
                                     aTGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTAc                            <  1:5662606/71‑1 (MQ=255)
                                       gcgcTGATCTCTTCCGGTACGACGCCAAAACTGCTGGCTAATGAAACCGACGCGCGTTTCATCGGCTACgg                          >  1:1947638/1‑71 (MQ=255)
                                        cgcTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTAt                                   >  1:1328776/1‑61 (MQ=255)
                                        cgcTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACgg                          >  1:516288/1‑70 (MQ=255)
                                        cgcTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGt                         >  1:1336817/1‑71 (MQ=255)
                                         gcTGATCTCTTCTGGGACCACGCCGAAGATGCAGGcgaacga                                                     <  1:4413063/42‑1 (MQ=255)
                                          cTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg                       <  1:4844150/71‑1 (MQ=255)
                                          cTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg                       <  1:1495428/71‑1 (MQ=255)
                                          cTGATCTCTTCCGGTACGACGCCAAAACTGCTGGCTAATGAAACCGACGCGCGTTTCATCGGCTACggcgc                       <  1:4059700/71‑4 (MQ=255)
                                           tGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg                       <  1:3426164/70‑1 (MQ=255)
                                            gATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacgaa                                                    >  1:50421/1‑40 (MQ=255)
                                            gATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacgaa                                                    >  1:4492211/1‑40 (MQ=255)
                                             aTCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGAt                    >  1:3010674/1‑71 (MQ=255)
                                               ctctTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACgg                          >  1:3112958/1‑63 (MQ=255)
                                               ctctTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGt                  >  1:1357813/1‑71 (MQ=255)
                                                tcttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                                      <  1:1965076/50‑1 (MQ=255)
                                                tcttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                                      <  1:4180707/50‑1 (MQ=255)
                                                 cttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGt                  >  1:4686747/1‑69 (MQ=255)
                                                 cttctGCGACCACGCCGAAGATGCTGGAGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGt                  >  1:2410232/1‑69 (MQ=255)
                                                  ttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTaa               >  1:25075/1‑71 (MQ=255)
                                                   tctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGtt                 <  1:2364841/68‑1 (MQ=255)
                                                    ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTa                <  1:2760664/68‑1 (MQ=255)
                                                    ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTa                <  1:4899299/68‑1 (MQ=255)
                                                    ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTa                <  1:2013463/68‑1 (MQ=255)
                                                    ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTa                <  1:3034470/68‑1 (MQ=255)
                                                    ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTa                <  1:4011975/68‑1 (MQ=255)
                                                       ggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg                       <  1:761870/58‑1 (MQ=255)
                                                       ggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATgg            <  1:2247715/69‑1 (MQ=255)
                                                       ggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATgg            <  1:5649198/69‑1 (MQ=255)
                                                        gACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAAt         <  1:440833/71‑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:409881/68‑1 (MQ=255)
                                                             caccGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTGATGGAATCCtt     <  1:1252535/68‑1 (MQ=255)
                                                             caccGAAGATGCTGGCGAACGAAGGGCAGACGTGCTTTATCGGCTACGGTGGGATGTTGATGGAATCCtt     <  1:5587177/68‑1 (MQ=255)
                                                             cGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTc    >  1:3637276/1‑71 (MQ=255)
                                                             cGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTc    <  1:1682076/71‑1 (MQ=255)
                                                              gCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCg   >  1:5626112/1‑71 (MQ=255)
                                                              gCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCg   <  1:898131/71‑1 (MQ=255)
                                                               ccGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATGCTTCGt  <  1:6528/71‑1 (MQ=255)
                                                               ccGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGt  <  1:4641613/71‑1 (MQ=255)
                                                                cGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg                       <  1:906365/49‑1 (MQ=255)
                                                                 gAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATcc       >  1:4508440/1‑64 (MQ=255)
                                                                     |                                                                
TTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGT  >  Exported/385778‑385911

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: