Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 385,957 G→T 17.4% V389V (GTG→GTT) ‡ cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported385,9570GT17.4% 135.0 / 13.4 53V389V (GTG→GTT) ‡cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (18/25);  new base T (6/3);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

GGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTG  >  Exported/385889‑386024
                                                                    |                                                                   
gggATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTa                                                                   >  1:4061144/1‑71 (MQ=255)
  gATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAtt                                                                 <  1:3964663/71‑1 (MQ=255)
   aTGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATGGATCCGGGCGTGTAttt                                                                >  1:319888/1‑71 (MQ=255)
   aTGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAttt                                                                >  1:1417581/1‑71 (MQ=255)
     gTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTg                                                              >  1:1918695/1‑71 (MQ=255)
      ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGc                                                             >  1:3382033/1‑71 (MQ=255)
        aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAttt                                                                >  1:957108/1‑66 (MQ=255)
        aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGc                                                             >  1:3199603/1‑69 (MQ=255)
        aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCa                                                           >  1:442668/1‑71 (MQ=255)
        aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCa                                                           >  1:4313481/1‑71 (MQ=255)
        aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCa                                                           >  1:2477081/1‑71 (MQ=255)
         aTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCa                                                           >  1:2778469/1‑70 (MQ=255)
         aTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCAt                                                          <  1:4172047/71‑1 (MQ=255)
         aTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCAt                                                          <  1:3215696/71‑1 (MQ=255)
          tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAttt                                                                >  1:4449895/1‑64 (MQ=255)
          tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                         <  1:4695503/71‑1 (MQ=255)
          tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                         <  1:3588898/71‑1 (MQ=255)
          tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                         <  1:5190215/71‑1 (MQ=255)
          tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                         <  1:5561103/71‑1 (MQ=255)
           ggAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTCTATTTTGCCATg                                                         >  1:5022132/1‑70 (MQ=255)
            gAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGa                                                        <  1:1944747/70‑1 (MQ=255)
               tCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACac                                                    <  1:1128206/71‑2 (MQ=255)
                 cTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGcc                                                  >  1:144419/1‑71 (MQ=255)
                 cTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACAccc                                                  >  1:5577617/1‑71 (MQ=255)
                   tCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCg                                                >  1:1968933/1‑71 (MQ=255)
                    cGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc                                                 <  1:5523451/69‑1 (MQ=255)
                    cGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc                                                 >  1:4800594/1‑69 (MQ=255)
                    cGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc                                               >  1:2396515/1‑70 (MQ=255)
                    cGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc                                               >  1:1954198/1‑70 (MQ=255)
                    cGTGGCGATTACGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc                                                 <  1:1744653/69‑1 (MQ=255)
                        gCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc                                               <  1:1177082/67‑2 (MQ=255)
                        gCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc                                               <  1:1170488/67‑2 (MQ=255)
                           aTTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGt                                        <  1:5356536/71‑1 (MQ=255)
                             tATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACAccc                                                  >  1:2232338/1‑59 (MQ=255)
                             tATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACAccc                                                  >  1:122627/1‑59 (MQ=255)
                               tGGCGCTGGTTTCCGCCTGTATCATCGAGCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCt                                     <  1:2370925/70‑1 (MQ=255)
                                    cTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATg                                             >  1:4033427/1‑57 (MQ=255)
                                      ggTTTCCGCCTGTATCATCGATCCTGGCGTGTATTTTGCCATGAACAGccc                                                 <  1:46340/51‑1 (MQ=255)
                                      ggTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc                                                 <  1:1721911/51‑1 (MQ=255)
                                          tCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                         <  1:5483642/39‑1 (MQ=255)
                                          tCCGCCTGTATCATCGATCCGGGCGTGTATTTAGCCATg                                                         <  1:1252400/39‑1 (MQ=255)
                                              ccTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGAcg                     >  1:3031944/1‑71 (MQ=255)
                                                  tATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATg                                             <  1:764686/43‑1 (MQ=255)
                                                   atcatcGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATg                                             <  1:3170218/42‑1 (MQ=255)
                                                    tcatcGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGAt               >  1:3087739/1‑71 (MQ=255)
                                                           tCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCg        <  1:3520844/71‑1 (MQ=255)
                                                           tCCGGGCGTGTATTTTACCATTAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCg        <  1:1715959/71‑1 (MQ=255)
                                                            ccGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGc       >  1:5260507/1‑71 (MQ=255)
                                                             cGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCt      <  1:4650254/71‑1 (MQ=255)
                                                               ggCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTc    <  1:336755/71‑1 (MQ=255)
                                                                gCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCt   <  1:3744493/71‑1 (MQ=255)
                                                                 cGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTg  <  1:3751628/71‑1 (MQ=255)
                                                                 cGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTg  <  1:48219/71‑1 (MQ=255)
                                                                    |                                                                   
GGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTG  >  Exported/385889‑386024

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: