Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,114 G→C 26.1% A108A (GCG→GCC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,1140GC26.1% 75.2 / 26.9 46A108A (GCG→GCCcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (11/23);  new base C (4/8);  total (15/31)
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.85e-01

AAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTG  >  minE/385047‑385180
                                                                   |                                                                  
acgtgtTGTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCgccgcg                                                                 >  1:2778902/3‑71 (MQ=255)
aaGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTggcggc                                                                  <  1:952885/70‑1 (MQ=255)
     ttGTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGAt                                                            <  1:2429710/70‑1 (MQ=255)
     cTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCaaa                                                             >  1:1152441/1‑70 (MQ=255)
        ttCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATggg                                                         >  1:491594/1‑71 (MQ=255)
        ttCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATg                                                           <  1:2390870/69‑1 (MQ=255)
            gTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTAc                                                     <  1:788121/71‑1 (MQ=255)
               aCCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTg                                                  >  1:633270/1‑71 (MQ=255)
               aCCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTg                                                  <  1:3351358/71‑1 (MQ=255)
                ccACTTCGCCGCTATCGCCGGTGCTGGTTCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGc                                                 <  1:3578533/71‑1 (MQ=255)
                ccACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGc                                                 <  1:2743310/71‑1 (MQ=255)
                ccACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATAGGCTACCTGc                                                 <  1:1676088/71‑1 (MQ=255)
                    tttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCggg                                             >  1:2037630/1‑71 (MQ=255)
                     ttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCgg                                              <  1:2008806/69‑1 (MQ=255)
                     ttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGAGGCGCAAATGGGCTACCTGCCggg                                             <  1:938978/70‑1 (MQ=255)
                        cGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGgatg                                         <  1:3329337/71‑1 (MQ=255)
                          gCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGgatgat                                       <  1:2355728/71‑1 (MQ=255)
                           ccATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGcc                                                <  1:1011411/61‑1 (MQ=255)
                            cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATTGGCTACCTGCCGGGGATGATCt                                     <  1:3008884/71‑1 (MQ=255)
                            cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCt                                     <  1:2075286/71‑1 (MQ=255)
                            cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCt                                     >  1:1832261/1‑71 (MQ=255)
                            cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCt                                     <  1:3106573/71‑1 (MQ=255)
                                 ccGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTAc                                                     >  1:3760422/1‑50 (MQ=255)
                                  cGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCag                                                              <  1:1982117/40‑2 (MQ=12)
                                    gAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGcc                                                <  1:961999/52‑1 (MQ=255)
                                    gAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGgctgct                              <  1:3736903/70‑1 (MQ=255)
                                         ggTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTgg                        <  1:2535290/71‑1 (MQ=255)
                                          gTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGcgg                         <  1:2764584/69‑3 (MQ=255)
                                          gTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCgggg                       >  1:3133188/1‑71 (MQ=255)
                                          gTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCgggg                       >  1:696918/1‑71 (MQ=255)
                                          gTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGc                           <  1:450810/67‑1 (MQ=255)
                                              gCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGggtgg                   >  1:2151925/1‑71 (MQ=255)
                                               cTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTg                         <  1:2541821/64‑1 (MQ=255)
                                                tggtggGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGtt                 >  1:1105524/1‑71 (MQ=255)
                                                     gggCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGgctgc                               <  1:532101/52‑1 (MQ=255)
                                                     gggCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGc            <  1:3899209/71‑1 (MQ=255)
                                                        ccGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCcgg         >  1:1007658/1‑71 (MQ=255)
                                                        ccGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCcgg         >  1:642165/1‑71 (MQ=255)
                                                        ccGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcgg         <  1:1456049/71‑1 (MQ=255)
                                                        ccGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcgg         <  1:887719/71‑1 (MQ=255)
                                                         cGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCcgg         <  1:1106663/70‑1 (MQ=255)
                                                           gTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGgct                                 <  1:3381243/44‑1 (MQ=255)
                                                             aCTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTggg                       <  1:3414683/52‑1 (MQ=255)
                                                             aCTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTggg                       <  1:3049909/52‑1 (MQ=255)
                                                             aCTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcgg    >  1:3593323/1‑71 (MQ=255)
                                                              cTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGgct                                 >  1:2079447/1‑41 (MQ=255)
                                                              cTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggt   >  1:51912/1‑71 (MQ=255)
                                                              cTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggt   <  1:2969541/71‑1 (MQ=255)
                                                                ggcggcGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggtg  >  1:2952552/1‑70 (MQ=255)
                                                                 gcggcgCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggtg  >  1:3180320/1‑69 (MQ=255)
                                                                 gcggcgCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggtg  <  1:920041/69‑1 (MQ=255)
                                                                   |                                                                  
AAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTG  >  minE/385047‑385180

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: