Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 613,789 C→G 16.9% pseudogene (79/987 nt) cstA → carbon starvation protein; b0598_2

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606613,7890C→G16.9% 2.2 ‑3.5 ‑0.040pseudogene (79/987 nt)cstAcarbon starvation protein; b0598_2
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (17/16):  total (24/16)
Fisher's exact test strand distribution p-value = 2.95e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.89e-01

GCGTTTACTTTGCCATGAACAGCCCAATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCT  >  REL606/613754‑613824
                                   |                                   
gCGTTTACTTTGCCATGAACAGCCCAATGGCGGTGc                                     <  1:35741/36‑1 (MQ=255)
  gTTTACTTTGCCATGAACAGCCCAATGGCGGTGCTg                                   >  1:4657568/1‑36 (MQ=255)
  gTTTACTTTGCCATGAACAGCCCAATGGCGGTGCTg                                   <  1:3178941/36‑1 (MQ=255)
   tttACTTTGCCATGAACAGCCCAATGGCGGTGGTgg                                  >  1:4984944/1‑36 (MQ=255)
     tACTTTGCCATGAACAGCCCAATGGCGGTGGTgggt                                >  1:5997016/1‑34 (MQ=38)
     tACTTTGCCATGAACAGCCCAATGGCGGTGCTGGCt                                <  1:5611352/36‑1 (MQ=255)
     tACTTTGCCATGAACAGCCCAATGGCGGTGCTGGCt                                <  1:2008721/36‑1 (MQ=255)
      aCTTTGCCATGAACAGCCCAATGGCGGTGGTGGGTc                               >  1:6036233/1‑36 (MQ=37)
       cTTTGCCATGAACAGCCCAATGGCGGTGCTGGCTcc                              <  1:5037947/36‑1 (MQ=255)
       cTTTGCCATGAACAGCCCAATGGCGGTGCTGGCTcc                              <  1:5892448/36‑1 (MQ=255)
        tttGCCATGAACAGCCCAATGGCGGGGGTGGGTCCg                             >  1:5926647/1‑36 (MQ=25)
          tGCCATGAACAGCCCAATGGCGGTGCTGGCTCCGGc                           <  1:2058743/36‑1 (MQ=255)
           gCCATGAACAGCCCAATGGCGGTGGTGGGTCCGGcc                          >  1:3408785/1‑35 (MQ=37)
           gCCATGAACAGCCCAATGGCGGTGCTGGCTCCGGCa                          <  1:8448310/36‑1 (MQ=255)
           gCCATGAACAGCCCAATGGCGGGGCTGGGGCCGGcg                          >  1:2511623/1‑35 (MQ=21)
            ccATGAACAGCCCAATGGCGGGGGTTGCTCCGGccg                         >  1:5425359/1‑34 (MQ=21)
            ccATGAACAGCCCAATGGCGGGGCTGGCGTCGGccg                         >  1:318393/1‑34 (MQ=18)
              aTGAACAGCCCAATGGCGGGGCTGGCTCCGGCCggg                       >  1:5777319/1‑36 (MQ=37)
               tGAACAGCCCAATGGCGGTGGTGGGTCCGGCAgggg                      >  1:1209918/1‑35 (MQ=37)
               tGAACAGCCCAATGGCGGTGCTGGCTCCGGCAgggg                      >  1:7880674/1‑35 (MQ=255)
                 aaCAGCCCAATGGCGGTGCTTGGTCCGGCCgggggg                    >  1:991956/1‑33 (MQ=19)
                 aaCAGCCCAATGGCGGTGCTGGCTCCGGCGgggggg                    >  1:1843798/1‑33 (MQ=37)
                 aaCAGCCCAATGGCGGTGCTGGCTCCGGCAGGGAcg                    <  1:2338691/36‑1 (MQ=255)
                  aCAGCCCAATGGCGGTGCTGGCTCCGGCAGGGAcgg                   <  1:5935742/36‑1 (MQ=255)
                  aCAGCCCAATGGCGGGGCTGGCTCCGGCTGGGAcgg                   >  1:8337018/1‑36 (MQ=37)
                  aCAGCCCAATGGCGGGGCCGGCTCCCGCGGGGGcgg                   >  1:6986357/1‑36 (MQ=16)
                    aGCCCAAGGGCGGTGCTGGCTCCGGCAGGGAcggcg                 <  1:4821368/36‑1 (MQ=255)
                     gCCCAATGGCGGTGCTGGCTCCGGCAGGGAcggcgg                <  1:2923034/36‑1 (MQ=255)
                     gCCCAATGGCGGTGCTGGCTCCGGCAGGGAcggcgg                <  1:6538324/36‑1 (MQ=255)
                        cAATGGCGGTGCTGGCTCCGGCAGGGACGGCGGAtg             <  1:8147655/36‑1 (MQ=255)
                             gCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGtg        >  1:2303855/1‑35 (MQ=255)
                             gCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTc        <  1:5346566/36‑1 (MQ=255)
                             gCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTc        >  1:4404858/1‑36 (MQ=255)
                             gCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTc        >  1:1795734/1‑36 (MQ=255)
                               ggTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGc      >  1:3842211/1‑36 (MQ=255)
                               ggTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGc      >  1:5720133/1‑36 (MQ=255)
                               ggTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGc      >  1:3881159/1‑36 (MQ=255)
                               ggTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGc      >  1:1277200/1‑36 (MQ=255)
                                gTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCggt     >  1:2132288/1‑34 (MQ=255)
                                 tGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCtt    <  1:4330048/36‑1 (MQ=255)
                                  gCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTc   <  1:2943772/36‑1 (MQ=255)
                                   cTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCt  >  1:7033093/1‑36 (MQ=255)
                                   |                                   
GCGTTTACTTTGCCATGAACAGCCCAATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCT  >  REL606/613754‑613824

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 13 ≤  ATCG  < 28 ≤  ATCG  < 37 ≤  ATCG