Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 385,165 T→G 30.8% V125V (GTT→GTG cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported385,1650TG30.8% 37.1 / 45.8 52V125V (GTT→GTGcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (14/22);  new base G (4/12);  total (18/34)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.50e-01

GTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGT  >  Exported/385097‑385234
                                                                    |                                                                     
gtggGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTgct                                                                     >  1:3439234/1‑71 (MQ=255)
gtggGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTgct                                                                     >  1:449019/1‑71 (MQ=255)
   ggtccGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGc                                                                  >  1:3274463/4‑71 (MQ=255)
    ggCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGcc                                                                 <  1:5507184/71‑1 (MQ=255)
           aCTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcg                                                           >  1:4901219/1‑70 (MQ=255)
            cTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTtctc                                                                    >  1:2271265/1‑60 (MQ=255)
              ggcggcGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTtctc                                                                    >  1:3012489/1‑58 (MQ=255)
               gcggcgCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggtg                                                        >  1:1625793/1‑69 (MQ=255)
               gccgcgcAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCcggtgcggt                                                         <  1:1569900/65‑1 (MQ=255)
               gccgcgcAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCcggtgcggt                                                         <  1:12417/65‑1 (MQ=255)
                  gcgcAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAg                                                     <  1:2642605/69‑1 (MQ=255)
                   cgcAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGac                                                  <  1:4336456/71‑2 (MQ=255)
                   cgcAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGac                                                  <  1:3892571/71‑2 (MQ=255)
                    gcAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAtt                                                 >  1:5187398/1‑71 (MQ=255)
                     cagaTGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGActt                                                <  1:4055455/68‑1 (MQ=255)
                     cAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAttt                                                <  1:1112480/71‑1 (MQ=255)
                       gaTGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGActtta                                              <  1:5612495/70‑3 (MQ=255)
                       gaTGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGActtta                                              >  1:3683789/2‑69 (MQ=255)
                       gaTGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGActtta                                              <  1:1917284/70‑3 (MQ=255)
                         tGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGACTTTATg                                            <  1:4973445/71‑1 (MQ=255)
                         tGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTAGCCGGTGCGGTGCAGGATTTCATg                                            >  1:5038970/1‑71 (MQ=255)
                              tACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGACTTTATGGTGCt                                       >  1:268300/1‑71 (MQ=255)
                              tACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGACTTTATGGTGCt                                       >  1:2429487/1‑71 (MQ=255)
                                ccTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCtgt                                     <  1:3102020/71‑1 (MQ=255)
                                 cTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCtgtt                                    >  1:4242997/1‑71 (MQ=255)
                                  tGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcg                                                                >  1:5497278/1‑42 (MQ=255)
                                  tGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcg                                                                >  1:3477614/1‑42 (MQ=255)
                                    ccGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAttt                                                <  1:808785/56‑1 (MQ=255)
                                    ccGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCtgtttgt                                 >  1:4970713/1‑71 (MQ=255)
                                      ggggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCa                                              >  1:890358/1‑56 (MQ=255)
                                      ggggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCa                                              >  1:3114377/1‑56 (MQ=255)
                                       gggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtg                                                             <  1:4303072/40‑1 (MQ=255)
                                       gggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtg                                                             <  1:4647609/40‑1 (MQ=255)
                                       gggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtg                                                             <  1:247791/40‑1 (MQ=255)
                                        ggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCt                             <  1:416204/71‑1 (MQ=255)
                                        ggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCt                             <  1:3221073/71‑1 (MQ=255)
                                         gatgatCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCt                             <  1:356398/70‑1 (MQ=255)
                                            gatCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTAcgc                         >  1:4146158/1‑71 (MQ=255)
                                             atCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAttt                                                <  1:3860205/47‑1 (MQ=255)
                                             atCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAttt                                                <  1:3721332/47‑1 (MQ=255)
                                                tGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGTTGCAGGAtt                                                 <  1:3479728/43‑1 (MQ=255)
                                                  gctgctGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGACTTTATGGTGCTGTTTATCTCCTCTCGCCGt                    <  1:3090677/70‑1 (MQ=255)
                                                  gctgctGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGACTTTATGGTGCTGTTTATCTCCTCTCGCCGt                    <  1:1542498/70‑1 (MQ=255)
                                                  gctgctGGCGGGGGGCGTGCTGGCCGGTGCGGTTCAGGACTTTATGGTGCTGTTTATCTCCTCTCGCCGt                    <  1:500199/70‑1 (MQ=255)
                                                  gctgctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGcc                      <  1:4305117/68‑1 (MQ=255)
                                                  gctgctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGt                    <  1:4418362/70‑1 (MQ=255)
                                                  gctgctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGt                    <  1:3200192/70‑1 (MQ=255)
                                                  gctgctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGt                    <  1:2458152/70‑1 (MQ=255)
                                                   ctgctgGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGACTTTATGGTGCtgttt                                   <  1:1922410/54‑1 (MQ=255)
                                                        cGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGtc             >  1:5490728/1‑71 (MQ=255)
                                                         gCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGtcg            >  1:3046235/1‑71 (MQ=255)
                                                            ggggTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGt              <  1:5310234/66‑1 (MQ=255)
                                                             gggTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGtcgctc         <  1:1595512/70‑1 (MQ=255)
                                                              ggtggtTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGtcgctcgc       >  1:3698662/1‑71 (MQ=255)
                                                               gtggtTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTcgctcgct      <  1:1056214/71‑1 (MQ=255)
                                                                  gtTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTa                            <  1:1259021/46‑1 (MQ=255)
                                                                    tctcGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGt  >  1:4562062/1‑70 (MQ=255)
                                                                    tctcGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGt  >  1:1874767/1‑70 (MQ=255)
                                                                    tctcGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGt  >  1:4414246/1‑70 (MQ=255)
                                                                    |                                                                     
GTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGT  >  Exported/385097‑385234

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: