Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 386,427 C→G 27.5% T546S (ACT→AGT) ‡ cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported386,4270CG27.5% 132.3 / 52.6 80T546S (ACT→AGT) ‡cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base C (32/26);  new base G (18/4);  total (50/30)
Fisher's exact test for biased strand distribution p-value = 3.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

CGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTG  >  Exported/386358‑386497
                                                                     |                                                                      
cgcTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACAct                                                                       >  1:4569652/1‑71 (MQ=255)
cgcTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACAct                                                                       <  1:4234938/71‑1 (MQ=255)
cgcTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACAct                                                                       <  1:3681459/71‑1 (MQ=255)
cgcTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACAct                                                                       <  1:5171616/71‑1 (MQ=255)
  cTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACActct                                                                     >  1:1258156/1‑71 (MQ=255)
  cTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACActct                                                                     >  1:1434803/1‑71 (MQ=255)
   tgtgCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACAct                                                                       >  1:2042437/1‑68 (MQ=255)
   tgtgCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACAct                                                                       >  1:4495123/1‑68 (MQ=255)
     tgCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCtg                                                                    <  1:4539334/69‑1 (MQ=255)
     tgCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGACCCGCTGGGCGGCATTAACACTCtg                                                                    <  1:5365051/69‑1 (MQ=255)
         tACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCg                                                              <  1:4755316/71‑1 (MQ=255)
           cTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCtgtg                                                                  >  1:2113715/1‑65 (MQ=255)
             ggCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGt                                                          >  1:5405944/1‑71 (MQ=255)
              ctgTGGGGCTACCTGCTGTATCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGtt                                                         >  1:4246239/3‑71 (MQ=255)
               cGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGttt                                                        <  1:1226240/71‑1 (MQ=255)
                 tGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTgg                                                      <  1:461061/71‑1 (MQ=255)
                 tGGGGCTACCTGCTGTATCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTgg                                                                 >  1:3887993/1‑60 (MQ=255)
                 tGGGGCTACCTGCTGTATCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTgg                                                                 >  1:1617627/1‑60 (MQ=255)
                 tGGGGCTACCTGCTGTATCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGttcg                                                       <  1:4717093/70‑3 (MQ=255)
                 tGGGGCTACCTGCTGTATCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCgg                                                      >  1:631643/1‑71 (MQ=255)
                    ggCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAt                                                   <  1:143491/71‑1 (MQ=255)
                    ggCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAt                                                   >  1:3791246/1‑71 (MQ=255)
                        aCTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACActc                                                                      <  1:2358133/48‑1 (MQ=255)
                         cTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCtg                                                                    >  1:2041384/1‑49 (MQ=255)
                         cTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCt                                                            <  1:3705264/57‑1 (MQ=255)
                          ttCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGt                                                          <  1:4154420/58‑1 (MQ=255)
                          ttCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCaa                                             >  1:3665008/1‑71 (MQ=255)
                           tcctccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAAc                                            >  1:568502/1‑71 (MQ=255)
                            cctccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGttt                                                        >  1:1398909/1‑58 (MQ=255)
                            cctccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTa                                                    >  1:405259/1‑62 (MQ=255)
                              tgtaTCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTAt                                                   >  1:469441/4‑61 (MQ=255)
                               ccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTa                                                    >  1:5199583/1‑59 (MQ=255)
                               ccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAt                                                   >  1:4503111/1‑60 (MQ=255)
                               ccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAt                                                   >  1:2478397/1‑60 (MQ=255)
                                taTCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGa                                        <  1:4894737/69‑1 (MQ=255)
                                cATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAtt                                                  <  1:675295/60‑1 (MQ=255)
                                 aTCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGAt                                       >  1:3166269/1‑70 (MQ=255)
                                 aTCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATg                                      >  1:5089241/1‑71 (MQ=255)
                                  tCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCt                                                            >  1:5307376/1‑48 (MQ=21)
                                  tCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCt                                                            >  1:4005202/1‑48 (MQ=21)
                                  tCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCt                                                            >  1:2601680/1‑48 (MQ=21)
                                  tCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCgg                                                      >  1:21154/1‑54 (MQ=25)
                                  tCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGc                                     <  1:3424525/71‑1 (MQ=255)
                                  tCAGGGCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGc                                     >  1:512031/1‑71 (MQ=255)
                                  tCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTgg                                                      <  1:4409241/54‑1 (MQ=255)
                                  tCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTgg                                                      <  1:5247843/54‑1 (MQ=255)
                                  tCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAAc                                            <  1:2615699/64‑1 (MQ=255)
                                  tCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAg                                         <  1:3803587/67‑1 (MQ=255)
                                  tCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGc                                     >  1:3092143/1‑71 (MQ=255)
                                    aGGTCGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGCTg                                   >  1:4025979/1‑71 (MQ=19)
                                    aGGGCGTGGTCGATCCGCTGGGCGGGATCAACACCCTGTGGCCGCTGTTCGGCATCGCCAACCAGATGCTg                                   >  1:598713/1‑71 (MQ=255)
                                    aGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTa                                                    >  1:4872730/1‑54 (MQ=255)
                                    aGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAtt                                                  >  1:2362930/1‑56 (MQ=255)
                                    aGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTg                                   >  1:764011/1‑71 (MQ=255)
                                    aGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTg                                   >  1:892867/1‑71 (MQ=255)
                                        cGTGGTTGATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGCTGGCAg                               >  1:5376531/1‑71 (MQ=25)
                                        cGTGGTCGATCCGCTGGGCGGCATTAACACTCtgt                                                                   >  1:315017/1‑35 (MQ=255)
                                        cGTGGTCGATCCGCTGGGCGGCATTAACACTCtgt                                                                   >  1:2767976/1‑35 (MQ=255)
                                        cGTGGTCGATCCGCTGGGCGGCATTAACACTCtgt                                                                   >  1:2761479/1‑35 (MQ=255)
                                        cGTGGTCGATCCGCTGGGCGGCATTAACACTCtgt                                                                   >  1:3416126/1‑35 (MQ=255)
                                         gTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAgg                              <  1:705262/71‑1 (MQ=255)
                                         gTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAgg                              >  1:5403478/1‑71 (MQ=255)
                                            gTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGttt                                                        >  1:793803/1‑42 (MQ=255)
                                            gTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGttt                                                        >  1:5442634/1‑42 (MQ=255)
                                            gTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGttt                                                        >  1:1523083/1‑42 (MQ=255)
                                            gTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGt                                                     >  1:4778075/1‑45 (MQ=255)
                                             ttgATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAAc                                            >  1:943000/3‑53 (MQ=16)
                                             ttgATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGCTgg                                  <  1:3797837/61‑1 (MQ=19)
                                              cGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTg                                                       >  1:4795395/1‑41 (MQ=255)
                                               gATCCGCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGt                                                     >  1:3596631/1‑42 (MQ=19)
                                               gATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTAt                                                   <  1:1805577/44‑1 (MQ=255)
                                               gATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGgc                        <  1:2984461/71‑1 (MQ=255)
                                                   cGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCt                     <  1:5492450/70‑1 (MQ=255)
                                                   cGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCt                     <  1:5444251/70‑1 (MQ=255)
                                                    gCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGa                                        >  1:2204434/1‑50 (MQ=255)
                                                    gCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGa                                        >  1:146458/1‑50 (MQ=255)
                                                    gCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGCTgg                                  >  1:4957401/1‑56 (MQ=16)
                                                    gCTGGGCGGCGTTAAGAGCCTGTGGCCGCTGTTCGGTATCTCCAACCAGATGCTgg                                  >  1:1378796/1‑56 (MQ=16)
                                                    gCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGAt                           <  1:2398335/63‑1 (MQ=255)
                                                    gCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCt                     <  1:2758244/69‑1 (MQ=255)
                                                    gCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCt                     <  1:146810/69‑1 (MQ=255)
                                                     cTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGAt                           <  1:490666/62‑1 (MQ=255)
                                                        ggcggcATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATg                 <  1:4941133/69‑1 (MQ=255)
                                                        ggcggcATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGc                <  1:1693235/70‑1 (MQ=255)
                                                               ttAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTgg                                  <  1:3298866/45‑1 (MQ=255)
                                                               ttAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGcc        <  1:4829369/71‑1 (MQ=255)
                                                                aaaCACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATg                                      >  1:4734756/2‑40 (MQ=255)
                                                                 aaCACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgt      >  1:3141143/1‑71 (MQ=255)
                                                                 aaCACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgt      >  1:1956186/1‑71 (MQ=255)
                                                                  acacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGc                                 >  1:1600008/1‑43 (MQ=255)
                                                                  acacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCa                                >  1:891339/1‑44 (MQ=255)
                                                                  acacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCa                                >  1:976486/1‑44 (MQ=255)
                                                                  acacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCa                                >  1:2459854/1‑44 (MQ=255)
                                                                  acacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCtgt           >  1:5453234/1‑65 (MQ=255)
                                                                  acacTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtg     >  1:22249/1‑71 (MQ=255)
                                                                    acTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtgg  <  1:2690552/71‑2 (MQ=255)
                                                                     ctctGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCgtggtg  <  1:4502642/71‑1 (MQ=255)
                                                                     |                                                                      
CGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTG  >  Exported/386358‑386497

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: