Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,643,860 A→G 26.8% R9R (CGT→CGC purC ← phosphoribosylaminoimidazole‑succinocarboxamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,643,8600AG26.8% 82.7 / 32.5 51R9R (CGT→CGCpurCphosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base A (13/24);  new base G (7/6);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTATCATTGCGGAATTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTT  >  minE/1643791‑1643926
                                                                     |                                                                  
cGTATCATTGCGGAATTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCAc                                                                   <  1:1588269/71‑1 (MQ=255)
 gTATCATTGCGGAATTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACcgcg                                                                  <  1:1861304/71‑1 (MQ=255)
 gTATCATTGCGGAATTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCAc                                                                   <  1:1117565/70‑1 (MQ=255)
  tATCATTGCGGAATTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGa                                                                 >  1:1332969/1‑71 (MQ=255)
    tCATTGCGGAATTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACcgcg                                                                  <  1:363138/68‑1 (MQ=255)
      aTTGCGGAATTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATAc                                                              >  1:798737/1‑70 (MQ=255)
      aTTGCGGAATTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATAc                                                              <  1:1863382/70‑1 (MQ=255)
      aTTGCGGAATTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATAc                                                              >  1:313957/1‑70 (MQ=255)
              aTTCGAGCACCAACAGGTCCGGGTTTTCCGTTCTGTATACGGTTTTCGCTTTACCACg                                                                  <  1:952215/58‑1 (MQ=255)
              aTTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACg                                                                  <  1:471719/58‑1 (MQ=255)
                   aGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAgcttgc                                                  <  1:646647/69‑1 (MQ=255)
                    gCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCttt                                               >  1:1270773/1‑71 (MQ=255)
                    gCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCttt                                               >  1:1549920/1‑71 (MQ=255)
                        cAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCAc                                                                   <  1:863743/47‑1 (MQ=255)
                        cAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAgcttgct                                                 <  1:1534442/65‑1 (MQ=255)
                        cAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCa                                           <  1:2142210/71‑1 (MQ=255)
                        cAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCa                                           <  1:1141701/71‑1 (MQ=255)
                        cAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCa                                           <  1:1716930/71‑1 (MQ=255)
                         aaCAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCGCGATACa                                                             >  1:1023249/1‑52 (MQ=255)
                         aaCAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGa                                                                 >  1:613048/1‑48 (MQ=255)
                         aaCAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTAC‑‑CGATACAACTCAGCTTGCTTTTGCATCt                                        >  1:1515664/1‑71 (MQ=255)
                          aCAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCGCGATACAACTCAGCTTGCTTTTGCAt                                          >  1:483284/1‑70 (MQ=255)
                            aGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCtt                                       <  1:1754235/71‑1 (MQ=255)
                               tCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCGCGATACAACTCAGCTTGCTTTTGCATCTTTAt                                    <  1:2091631/71‑1 (MQ=255)
                               tCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTAt                                    <  1:1124471/71‑1 (MQ=255)
                                    gTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAgc                                                      <  1:222922/48‑1 (MQ=255)
                                     ttttCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCt                                        <  1:839773/61‑1 (MQ=255)
                                      tttCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCtttt                                              >  1:1849371/1‑54 (MQ=255)
                                       ttCCGTACTGTATACGGTTTTCGCTTTACCGCGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTg                            >  1:282281/1‑71 (MQ=255)
                                          cGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTAt                                    <  1:325915/60‑1 (MQ=255)
                                            tGCTGTATACGGTTTTCGCTTTACCGCGTTACAACTCAgcttg                                                   <  1:1661025/43‑1 (MQ=38)
                                              cTGTATACGGTTTTCGCTTTACCGCGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGtgtg                      <  1:2062144/70‑1 (MQ=255)
                                                 tataCGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAAtt                  <  1:225346/71‑1 (MQ=255)
                                                   taCGGTTTTCGCTTTACCGCGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTaa                >  1:382649/1‑71 (MQ=255)
                                                   taCGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATc                                         <  1:1410156/46‑1 (MQ=255)
                                                   taCGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAAtt                  <  1:1386275/69‑1 (MQ=255)
                                                    aCGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAAc               >  1:2093117/1‑71 (MQ=255)
                                                     cGGTTTTCGCTTTACCGCGATACAACTCAGCTTGCTTTTGCATCTTTATCa                                  >  1:844622/1‑51 (MQ=255)
                                                     cGGTTTTCGCTTTACCGCGATACAACTCAGCTTGCTTTTGCATCTTTATCa                                  >  1:742207/1‑51 (MQ=255)
                                                     cGGTTTTCGCTTTACCACGATACAACTCAgcttgct                                                 <  1:2174388/36‑1 (MQ=255)
                                                     cGGTTTTCGCTTTACCACGATACAACTCAgcttgct                                                 <  1:38920/36‑1 (MQ=255)
                                                     cGGTTTTCGCTTTACCACGATACAACTCAgcttgct                                                 <  1:1242199/36‑1 (MQ=255)
                                                     cGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGaa                    >  1:1368447/1‑65 (MQ=255)
                                                     cGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAAc               >  1:1648400/1‑70 (MQ=255)
                                                     cGGTTTTAGCTTTACCGCGATACAACTCAGCTTGCTTTTGCATCTTTAt                                    >  1:1408593/1‑49 (MQ=39)
                                                      ggTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGc                                            <  1:2183515/40‑1 (MQ=255)
                                                        ttttCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTTAATTAACGttt           >  1:989222/1‑71 (MQ=255)
                                                          ttCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATc                                   >  1:2061398/1‑45 (MQ=255)
                                                           tCGCTTTACCGCGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTa         <  1:717983/70‑1 (MQ=255)
                                                               tttACCGCGATACAACTCAGCTTGCTTTTGCATCTTTATCACTcc                              <  1:2201408/45‑1 (MQ=255)
                                                                 tACCGCGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCt   <  1:1057296/70‑1 (MQ=255)
                                                                 tACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATc    <  1:560663/69‑1 (MQ=255)
                                                                 tACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCtt  >  1:1957350/1‑71 (MQ=255)
                                                                     |                                                                  
CGTATCATTGCGGAATTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTTTTAAAATCTT  >  minE/1643791‑1643926

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: