Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,120 A→G 35.2% Q110Q (CAA→CAG cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,1200AG35.2% 46.6 / 46.7 54Q110Q (CAA→CAGcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base A (18/17);  new base G (6/13);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTT  >  minE/385052‑385188
                                                                    |                                                                    
cTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAAt                                                                    >  1:868934/1‑71 (MQ=255)
 tGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAAt                                                                    <  1:5489501/70‑1 (MQ=255)
  gTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATg                                                                   <  1:5612341/70‑1 (MQ=255)
  gTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATg                                                                   <  1:5279289/70‑1 (MQ=255)
    ttggTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGc                                                                >  1:230749/3‑71 (MQ=255)
     cGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGc                                                                >  1:2273881/1‑70 (MQ=255)
      ggTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCt                                                               <  1:1631204/70‑1 (MQ=255)
       gTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGc                                                                >  1:204562/1‑68 (MQ=255)
          aCCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTg                                                          >  1:4949266/1‑71 (MQ=255)
          aCCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCt                                                           <  1:3629173/70‑1 (MQ=255)
           ccACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCt                                                               <  1:4808523/65‑1 (MQ=255)
           ccACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCt                                                               <  1:273080/65‑1 (MQ=255)
           ccACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGc                                                         >  1:3509648/1‑71 (MQ=255)
              ttttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCt                                                               <  1:647690/62‑1 (MQ=255)
              ttttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCgg                                                      >  1:373704/1‑71 (MQ=255)
              ttttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCgg                                                      >  1:5385301/1‑71 (MQ=255)
                ttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCggg                                                     >  1:2082745/1‑70 (MQ=255)
                ttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCggg                                                     >  1:283420/1‑70 (MQ=255)
                      ccATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATc                                              >  1:724685/1‑71 (MQ=255)
                           gCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTa                                                              >  1:2431413/1‑50 (MQ=255)
                            ccGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGgctg                                        <  1:3752196/71‑1 (MQ=255)
                             cGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGTGCAGATGGGCTACCTGCCTGGCACGCTGTGgctgc                                       <  1:2802760/71‑1 (MQ=255)
                                 gCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGgatgat                                               >  1:3264642/1‑59 (MQ=255)
                                 gCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGgatgat                                               >  1:2519884/1‑59 (MQ=255)
                                  cAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCt                                  <  1:1548247/71‑1 (MQ=255)
                                    ggTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCgggg                                                    <  1:455446/51‑1 (MQ=255)
                                      tCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGcg                                  >  1:4242236/1‑66 (MQ=255)
                                        cGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTAGCTGGggt                             <  1:4945212/70‑1 (MQ=255)
                                         gCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGggtgg                           <  1:2738382/71‑1 (MQ=255)
                                          cTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGggtcgt                          <  1:3437747/71‑1 (MQ=255)
                                             gtggGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGgct                                         <  1:499600/53‑1 (MQ=255)
                                             gtggGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTgct                       >  1:3439234/1‑71 (MQ=255)
                                             gtggGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTgct                       >  1:449019/1‑71 (MQ=255)
                                             gtggGTCCGGTTCGCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGgct                                         <  1:3781756/53‑1 (MQ=255)
                                                ggtccGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGc                    >  1:3274463/4‑71 (MQ=255)
                                                 ggCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGcc                   <  1:5507184/71‑1 (MQ=255)
                                                      gTACTGGCGGCGCAAATGGGCTACCTGCCGGGgatg                                                 <  1:1668014/36‑1 (MQ=255)
                                                      gTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCg                                    <  1:1648172/49‑1 (MQ=255)
                                                       tACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGggtgg                           <  1:5656408/57‑1 (MQ=255)
                                                       tACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGCCTGCTCGCTGGggtgg                           <  1:3690331/57‑1 (MQ=255)
                                                        aCTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTgg                                >  1:4122070/1‑51 (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)
                                                                  cAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTggg                               <  1:1156919/42‑1 (MQ=255)
                                                                  cAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTggg                               <  1:4314084/42‑1 (MQ=255)
                                                                  cAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAttt  <  1:1112480/71‑1 (MQ=255)
                                                                    aaTGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGtggtgg                           <  1:4952236/44‑1 (MQ=255)
                                                                    aaTGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGggtgg                           <  1:890215/44‑1 (MQ=255)
                                                                    |                                                                    
CTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTT  >  minE/385052‑385188

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: