Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,102 G→T 31.0% G104G (GGG→GGT cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,1020GT31.0% 90.2 / 25.7 51G104G (GGG→GGTcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (16/19);  new base T (6/10);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCC  >  minE/385034‑385171
                                                                    |                                                                     
ccGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGcc                                                                     >  1:489017/1‑71 (MQ=255)
 cGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCg                                                                    >  1:4572003/1‑71 (MQ=255)
  gACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCgg                                                                   <  1:554019/71‑1 (MQ=255)
  gACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCgg                                                                   <  1:1661392/71‑1 (MQ=255)
    ctgACAAGAAAGTGCTGTTCGGGCACCATTTTGCGGCCATTGCCGGAGCAGGGCCGCTGGTGGGGcc                                                                     <  1:3712220/65‑1 (MQ=255)
     ggACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTAc                                                                >  1:1046367/1‑71 (MQ=255)
        cAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTgg                                                             <  1:2712877/71‑1 (MQ=255)
         aaGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTggc                                                            >  1:3787673/1‑71 (MQ=255)
               gTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCa                                                      <  1:1289279/71‑1 (MQ=255)
               gTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCa                                                      <  1:1296041/71‑1 (MQ=255)
                 gttGTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCgccgcg                                                        <  1:1649463/65‑1 (MQ=255)
                  cTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAAt                                                   >  1:868934/1‑71 (MQ=255)
                   tGTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCa                                                      >  1:3884430/1‑67 (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)
                           cACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGt                                                                  <  1:1850941/47‑1 (MQ=255)
                           cACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGt                                                                  <  1:625748/47‑1 (MQ=255)
                            aCCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTg                                         >  1:4949266/1‑71 (MQ=255)
                            aCCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCt                                          <  1:3629173/70‑1 (MQ=255)
                             ccACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCt                                              <  1:273080/65‑1 (MQ=255)
                             ccACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCt                                              <  1:4808523/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)
                                         cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACt                                                               <  1:4810106/36‑1 (MQ=255)
                                          aTTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACt                                                               <  1:1446467/35‑1 (MQ=255)
                                          aTTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACt                                                               <  1:3462893/35‑1 (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:2519884/1‑59 (MQ=255)
                                                   gCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGgatgat                              >  1:3264642/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:449019/1‑71 (MQ=255)
                                                               gtggGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTgct      >  1:3439234/1‑71 (MQ=255)
                                                               gtggGTCCGGTTCGCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGgct                        <  1:3781756/53‑1 (MQ=255)
                                                                   ggCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGcc  <  1:5507184/71‑1 (MQ=255)
                                                                    |                                                                     
CCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCC  >  minE/385034‑385171

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: