Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,111 G→C 27.6% L107L (CTG→CTC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,1110GC27.6% 74.1 / 28.9 47L107L (CTG→CTCcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (11/23);  new base C (5/8);  total (16/31)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTG  >  minE/385042‑385180
                                                                     |                                                                     
cAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTgg                                                                      <  1:29184/71‑1 (MQ=255)
 aaGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTggc                                                                     >  1:1271111/1‑71 (MQ=255)
 aaGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTgg                                                                      >  1:3370545/1‑70 (MQ=255)
     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:1832261/1‑71 (MQ=255)
                                 cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCt                                     <  1:2075286/71‑1 (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:696918/1‑71 (MQ=255)
                                               gTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCgggg                       >  1:3133188/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:2969541/71‑1 (MQ=255)
                                                                   cTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggt   >  1:51912/1‑71 (MQ=255)
                                                                     ggcggcGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggtg  >  1:2952552/1‑70 (MQ=255)
                                                                     |                                                                     
CAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTG  >  minE/385042‑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: