Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,218 G→A 45.6% L280L (CTG→CTA kdsA → 3‑deoxy‑D‑manno‑octulosonate 8‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE821,2180GA45.6% 11.6 / 59.6 46L280L (CTG→CTAkdsA3‑deoxy‑D‑manno‑octulosonate 8‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (6/19);  new base A (9/12);  total (15/31)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

GGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCA  >  minE/821149‑821286
                                                                     |                                                                    
ggCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTg                                                                      >  1:3535663/1‑70 (MQ=255)
 gCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGa                                                                    <  1:2042952/71‑1 (MQ=255)
  cTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGAt                                                                   <  1:247021/71‑1 (MQ=255)
     aaCTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTAGATa                                                                  >  1:2837884/1‑69 (MQ=255)
          gAACGGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTAGATACCAGc                                                             <  1:2462716/69‑1 (MQ=255)
          gAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTAGATACCAGc                                                             <  1:3133148/69‑1 (MQ=255)
          gAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTAGATACCAGc                                                             <  1:1440703/69‑1 (MQ=255)
          gAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTAGATACCAGc                                                             <  1:1126486/69‑1 (MQ=255)
                ttCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATc                                                     <  1:108779/71‑1 (MQ=255)
                ttCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTAGATACCAGCAAGTaa                                                       <  1:3849632/69‑1 (MQ=255)
                ttCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTAGATACCAGCAAGTAATc                                                     >  1:2396890/1‑71 (MQ=255)
                      aaGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTAGATACCAGCAAGTAATCtttt                                                 <  1:3886230/69‑1 (MQ=255)
                                aaGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTg                                                                      <  1:2032651/38‑1 (MQ=255)
                                     aTTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGc                                                             <  1:1133044/42‑1 (MQ=255)
                                     aTTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAg                                <  1:2575250/71‑1 (MQ=255)
                                     aTTGATGATCTGGTGAAAGGTTTCGAAGAACTAGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAg                                <  1:3253185/71‑1 (MQ=255)
                                      ttGATGATCTGGTGAAAGGTTTCGAAGAACTGGAt                                                                   <  1:3056723/35‑1 (MQ=255)
                                      ttGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAg                                <  1:2201916/70‑1 (MQ=255)
                                       tgatgaTCTGGTGAAAGGTTTCGAAGAACTAGATAc                                                                 >  1:1258018/1‑36 (MQ=255)
                                        gatgatCTGGTTAAAGGTTTCGAAGAACTAGATACCAg                                                              <  1:1694600/38‑1 (MQ=37)
                                         atgatCTGGTGAAAGGTTTCGAAGAACTAGATACCAGCAAGTAATCtttt                                                 >  1:1154843/1‑50 (MQ=255)
                                         atgatCTGGTGAAAGGTTTCGAAGAACTAGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTAtt                            <  1:3381630/71‑1 (MQ=255)
                                         atgatCTGGTGAAAGGTTTCGAAGAACTAGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTAt                             >  1:3106363/1‑70 (MQ=255)
                                          tgatCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTa                           <  1:1448882/71‑1 (MQ=255)
                                                ggTGAAAGGTTTCGAAGAACTAGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCg                        <  1:3266827/68‑1 (MQ=255)
                                                    aaaGGTTTCGAAGAACTGGATACCAGCAAGTAATc                                                     <  1:2562324/35‑1 (MQ=255)
                                                    aaaGGTTTCGAAGAACTAGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCg                 >  1:3213534/1‑71 (MQ=255)
                                                      aGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTa                              >  1:206447/1‑56 (MQ=255)
                                                      aGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGt                >  1:514790/1‑70 (MQ=255)
                                                        gTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTaa             <  1:1384130/71‑1 (MQ=255)
                                                         tttCGAATAACTGGATACCAGCAAGTAATCTTTTTTGCTTg                                          >  1:1720918/1‑41 (MQ=255)
                                                         tttCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGtt               >  1:2710250/1‑68 (MQ=255)
                                                           tCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGAc          >  1:3190807/1‑71 (MQ=255)
                                                             gaagaaCTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATaa                                  <  1:1915414/45‑1 (MQ=255)
                                                             gaagaaCTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACtt        <  1:145530/71‑1 (MQ=255)
                                                              aagaaCTAGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATa                                   <  1:13090/43‑1 (MQ=255)
                                                               agaaCTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGc                         >  1:3970874/1‑52 (MQ=255)
                                                               agaaCTAGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACtttt      <  1:2691959/71‑1 (MQ=255)
                                                                gaaCTAGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACttttt     >  1:99741/1‑71 (MQ=255)
                                                                gaaCTAGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACttttt     >  1:1849774/1‑71 (MQ=255)
                                                                 aaCTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTa                           <  1:3661286/48‑1 (MQ=255)
                                                                 aaCTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACtttttt    <  1:3438667/71‑1 (MQ=255)
                                                                 aaCTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACtttttt    <  1:1963310/71‑1 (MQ=255)
                                                                  aCTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTa                           <  1:3187481/47‑1 (MQ=255)
                                                                  aCTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACtttttt    <  1:1297129/70‑1 (MQ=255)
                                                                  aCTAGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTa                              >  1:2646525/1‑44 (MQ=255)
                                                                   cTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCa  <  1:2602030/71‑1 (MQ=255)
                                                                    tGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATa                                   <  1:2667518/37‑1 (MQ=255)
                                                                     ggATACCAGCAAGTAATCTTTTTTGCTTGAAAAATaaa                                 >  1:3894447/1‑38 (MQ=255)
                                                                     ggATACCAGCAAGTAATCTTTTTTGCTTGAAAAATaaa                                 >  1:85489/1‑38 (MQ=255)
                                                                     |                                                                    
GGCTAAACTGGAACCGTTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAGCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCA  >  minE/821149‑821286

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: