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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported385,2380CA32.2% 95.6 / 41.9 62L150M (CTG→ATG) cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base C (22/20);  new base A (10/10);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCA  >  Exported/385169‑385307
                                                                      |                                                                     
gCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGc                                                                       <  1:2879691/70‑1 (MQ=255)
gCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCt                                                                      >  1:2012251/1‑71 (MQ=255)
 ccGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTg                                                                     <  1:5132589/71‑1 (MQ=255)
   ggtgcggtgcAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGt                                                                   >  1:3763399/1‑71 (MQ=255)
    gtgcggtgcAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTc                                                                  <  1:1909713/71‑1 (MQ=255)
        ggtgcAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAaag                                                              <  1:1055205/71‑1 (MQ=255)
         gtgcAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGCTCGCTCGCTGGGTGAGCTGGTCAaaga                                                             <  1:1892611/71‑1 (MQ=255)
              ggACTTTATGGTGCTGTTTATCTCCTCTCGCCGTAATGGCGCATCTCTTGGTGAGATGATCAAAGAagag                                                         <  1:1528555/70‑1 (MQ=255)
                aTTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATg                                                      >  1:3716572/1‑71 (MQ=255)
                 tttCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATgg                                                     >  1:3351490/1‑71 (MQ=255)
                 tttCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATgg                                                     <  1:3741813/71‑1 (MQ=255)
                   tCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAaagaa                                                            >  1:1069574/1‑62 (MQ=255)
                   tCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAaagaa                                                            >  1:741553/1‑62 (MQ=255)
                   tCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAaagaa                                                            >  1:4164555/1‑62 (MQ=255)
                   tCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATgggg                                                   >  1:1642879/1‑71 (MQ=255)
                   tCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATgggg                                                   >  1:2193180/1‑71 (MQ=255)
                    taTGGTGCTGTTTATCTCCTCTCGCCGTAATGGCGCATCTCTTGGTGAGATGATCAAAGAagag                                                         >  1:4429809/2‑64 (MQ=255)
                     aTGGTGCTGTTTATCTCCTCTCGCCGTAATGGCGCATCTCTTGGTGAGATGATCAAAGAAGAGATgggac                                                  <  1:4817021/70‑3 (MQ=255)
                     aTGGTGCTGTTTATCTCCTCTCGCCGTAATGGCGCATCTCTTGGTGAGATGATCAAAGAAGAGATgggac                                                  <  1:5493565/70‑3 (MQ=255)
                        gTGCTGTTTATCTCCTCTCGCCGTAATGGCGCATCTCTTGGTGAGATGATCAAAGAAGAGATGGGACCag                                               <  1:5512658/70‑2 (MQ=255)
                           cTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCaaa                                                               <  1:2023765/51‑1 (MQ=255)
                             gtttgtttCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAaaga                                                             >  1:3286166/1‑51 (MQ=255)
                              tttgtttCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCaaa                                                               <  1:2671661/48‑1 (MQ=255)
                              tttgtttCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCaaa                                                               <  1:4307058/48‑1 (MQ=255)
                              tttgtttCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAAccgccg                                         <  1:4726214/70‑1 (MQ=255)
                               ttgtttCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCaaa                                                               >  1:3585794/1‑47 (MQ=255)
                                tgtttCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCaa                                                                >  1:3947781/1‑45 (MQ=255)
                                tgtttCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCaa                                                                >  1:1913823/1‑45 (MQ=255)
                                tgtttCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAgaga                                                        >  1:5420126/1‑53 (MQ=255)
                                tgtttCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCgggg                                      >  1:4774703/1‑71 (MQ=255)
                                  tttCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAagaaga                                                          <  1:2877875/49‑1 (MQ=255)
                                    tCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAaagaa                                                            >  1:390412/1‑45 (MQ=255)
                                     cTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAAccg                                            <  1:37237/60‑1 (MQ=255)
                                         gcgcCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGcc                                                 <  1:122803/51‑1 (MQ=255)
                                         ctcgccgtaatggcgcatctcttGGTGAGATGATCAaag                                                              >  1:1745845‑M1/24‑39 (MQ=255)
                                         ctcgccgtaatggcgcatctcttGGTGAGATGATCAaag                                                              >  1:3188056‑M1/24‑39 (MQ=255)
                                         ctcgccgtaatggcgcatctcttGGTGAGATGATCAaag                                                              >  1:2846759‑M1/24‑39 (MQ=255)
                                         ctcgccgtaatggcgcatctcttGGTGAGATGATCAaag                                                              >  1:5621722‑M1/24‑39 (MQ=255)
                                         ctcgccgtaatggcgcatctcttGGTGAGATGATCAAagaaga                                                          >  1:3695404‑M1/24‑43 (MQ=255)
                                         ctcgccgtaatggcgcatctcttGGTGAGATGATCAAagaaga                                                          >  1:3106716‑M1/24‑43 (MQ=255)
                                          tcgccgtaatggcgcatctcttGGTGAGATGATCAAAGAagag                                                         <  1:4256615‑M1/21‑1 (MQ=255)
                                            gccgtaatggcgcatctcttGGTGAGATGATCAaagaa                                                            <  1:475511‑M1/18‑1 (MQ=255)
                                            gccgtaatggcgcatctcttGGTGAGATGATCAAAGAAGAGATGGGACcag                                               <  1:857612‑M1/31‑2 (MQ=255)
                                            ccGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCtgg                          <  1:5085332/71‑1 (MQ=255)
                                                taatggcgcatctcttGGTGAGATGATCAAagaaga                                                          >  1:3596421‑M1/17‑36 (MQ=255)
                                                gACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGc                      >  1:968373/1‑71 (MQ=255)
                                                gACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGc                      >  1:2406027/1‑71 (MQ=255)
                                                 aCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATgggg                                                   >  1:351411/1‑41 (MQ=255)
                                                  atggcgcatctcttGGTGAGATGATCAAAGAAGAGATGGGACcag                                               <  1:4448854‑M1/31‑2 (MQ=255)
                                                  atggcgcatctcttGGTGAGATGATCAAAGAAGAGATGGGACcag                                               <  1:4267046‑M1/31‑2 (MQ=255)
                                                  atggcgcatctcttGGTGAGATGATCAAAGAAGAGATGGGACcag                                               <  1:511134‑M1/31‑2 (MQ=255)
                                                  atggcgcatctcttGGTGAGATGATCAAAGAAGAGATGGGACCa                                                >  1:2087755‑M1/15‑44 (MQ=255)
                                                   tggcgcatctcttGGTGAGATGATCAAAGAAgaga                                                        >  1:4054219‑M1/14‑35 (MQ=255)
                                                            gctGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGa           >  1:4633605/1‑71 (MQ=255)
                                                             ctGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGAt          <  1:1695858/71‑1 (MQ=255)
                                                              tGGGTGAGCTGGTCAAAGAAGAGATTGGGCCAACCGCCGGGGTGATTg                                <  1:3204449/48‑1 (MQ=255)
                                                              tGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTg                                <  1:101812/48‑1 (MQ=255)
                                                               gggTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCtt                >  1:1233697/1‑63 (MQ=255)
                                                                   gAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGt                                     <  1:5362216/38‑1 (MQ=255)
                                                                   gAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGt                                     <  1:443114/38‑1 (MQ=255)
                                                                   gAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGt    >  1:2884073/1‑71 (MQ=255)
                                                                    aGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGt    <  1:5165371/70‑1 (MQ=255)
                                                                     gCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCa  <  1:634859/71‑1 (MQ=255)
                                                                      |                                                                     
GCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCA  >  Exported/385169‑385307

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: