Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 385,960 T→C 18.3% Y390Y (TAT→TAC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported385,9600TC18.3% 100.6 / 15.7 49Y390Y (TAT→TACcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (16/24);  new base C (6/3);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGC  >  Exported/385891‑386030
                                                                     |                                                                      
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:4313481/1‑71 (MQ=255)
      aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCa                                                                 >  1:442668/1‑71 (MQ=255)
      aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCa                                                                 >  1:2477081/1‑71 (MQ=255)
       aTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCa                                                                 >  1:2778469/1‑70 (MQ=255)
       aTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCAt                                                                <  1:3215696/71‑1 (MQ=255)
       aTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCAt                                                                <  1:4172047/71‑1 (MQ=255)
        tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAttt                                                                      >  1:4449895/1‑64 (MQ=255)
        tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                               <  1:5561103/71‑1 (MQ=255)
        tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                               <  1:5190215/71‑1 (MQ=255)
        tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                               <  1:3588898/71‑1 (MQ=255)
        tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                               <  1:4695503/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:4800594/1‑69 (MQ=255)
                  cGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc                                                       <  1:5523451/69‑1 (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:1170488/67‑2 (MQ=255)
                      gCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc                                                     <  1:1177082/67‑2 (MQ=255)
                         aTTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGt                                              <  1:5356536/71‑1 (MQ=255)
                           tATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACAccc                                                        >  1:122627/1‑59 (MQ=255)
                           tATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACAccc                                                        >  1:2232338/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)
                                                                   tATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcgc    >  1:3079273/1‑71 (MQ=255)
                                                                   tATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcg     >  1:963691/1‑70 (MQ=255)
                                                                     ttttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGTCGGATGTGGTCGCTTCTGCcgcg   <  1:1938472/70‑1 (MQ=255)
                                                                     ttttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGc             <  1:2033383/60‑1 (MQ=255)
                                                                     ttttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcgcg   >  1:19775/1‑70 (MQ=255)
                                                                     ttttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcgcg   <  1:4090815/70‑1 (MQ=255)
                                                                     ttttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcgcg   <  1:365381/70‑1 (MQ=255)
                                                                     ttttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCgcgc  <  1:5174180/71‑1 (MQ=255)
                                                                     |                                                                      
GATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGC  >  Exported/385891‑386030

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: