Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,094,789 C→T 10.4% G21G (GGG→GGA)  tdcG ← L‑serine dehydratase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,094,7890CT10.4% 208.3 / 11.5 77G21G (GGG→GGA) tdcGL‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base C (26/43);  new base T (6/2);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 6.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TGGCTCGTCGCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCA  >  minE/2094721‑2094856
                                                                    |                                                                   
tGGCTCGTCGCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGccc                                                                   <  1:4898725/71‑1 (MQ=255)
tGGCTCGTCGCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGccc                                                                   >  1:4839957/1‑71 (MQ=255)
     cgtcgCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACgt                                                              >  1:4200169/1‑70 (MQ=255)
      gtcgCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACgg                                                              >  1:3269073/1‑70 (MQ=255)
       tcgCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGt                                                             <  1:853704/70‑1 (MQ=255)
        cgcgGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTAt                                                           >  1:1860267/1‑71 (MQ=255)
           ggTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGt                                                             >  1:696469/1‑66 (MQ=255)
           ggTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCg                                                        <  1:1441762/71‑1 (MQ=255)
            gTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGa                                                       <  1:699096/71‑1 (MQ=255)
                ataaGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTg                                                   <  1:1327468/71‑1 (MQ=255)
                  aaGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTg                                                   >  1:4823722/1‑69 (MQ=255)
                   aGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGa                                                 <  1:3861389/70‑1 (MQ=255)
                   aGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGa                                                 <  1:1862327/70‑1 (MQ=255)
                   aGCCGCTACTTT‑CAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCa                                                                 >  1:70190/1‑53 (MQ=255)
                           cTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCaa                                         <  1:4652208/70‑1 (MQ=255)
                           cTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAAt                                        <  1:4763655/71‑1 (MQ=255)
                           cTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAAt                                        >  1:5313544/1‑71 (MQ=255)
                             ttCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATcc                                      >  1:5114875/1‑71 (MQ=255)
                              tCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATcc                                      <  1:1358814/70‑1 (MQ=255)
                              tCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATcc                                      >  1:5202564/1‑70 (MQ=255)
                              tCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGAACCAATcc                                      >  1:4056668/1‑70 (MQ=255)
                              tCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCACACGGTATGCGAACTGGAGGGACCAATcc                                      >  1:491368/1‑70 (MQ=255)
                               ccAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACt                                                    <  1:2402142/55‑1 (MQ=255)
                                cAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGc                                                         >  1:5368591/1‑49 (MQ=255)
                                      aTCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCg                                                        <  1:2889229/44‑1 (MQ=255)
                                       tCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaa                            <  1:4121931/71‑1 (MQ=255)
                                       tCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGa                             >  1:3134508/1‑70 (MQ=255)
                                         aaTAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaaaa                          <  1:4592910/71‑1 (MQ=255)
                                             aaaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCg                                                        <  1:5395617/37‑1 (MQ=255)
                                             aaaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCg                                                        <  1:3587636/37‑1 (MQ=255)
                                             aaaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGa                                                 >  1:1471810/1‑44 (MQ=255)
                                             aaaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATc                      <  1:5076281/71‑1 (MQ=255)
                                              aaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                     <  1:1765634/71‑1 (MQ=255)
                                              aaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                     <  1:3942990/71‑1 (MQ=255)
                                              aaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                     <  1:1802592/71‑1 (MQ=255)
                                               aCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                     <  1:4506995/70‑1 (MQ=255)
                                               aCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGa                    <  1:218863/71‑1 (MQ=255)
                                                 tttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAAt                                        >  1:4217547/1‑49 (MQ=255)
                                                  ttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                 >  1:5202255/1‑71 (MQ=255)
                                                  ttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                 >  1:3337653/1‑71 (MQ=255)
                                                  ttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                 >  1:2724242/1‑71 (MQ=255)
                                                   tttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAAt                                  >  1:3459377/1‑53 (MQ=255)
                                                   tttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAAt                                  >  1:5331385/1‑53 (MQ=255)
                                                       cGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACt             <  1:3990211/70‑1 (MQ=255)
                                                       cGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACt             <  1:4448373/70‑1 (MQ=255)
                                                        ggCATTCATTGGCCCCACGGTATGCGAACTGGAGGGAc                                            >  1:249946/1‑38 (MQ=255)
                                                        ggCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTaa           <  1:3967628/71‑1 (MQ=255)
                                                         gctttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTAAt          >  1:3441019/4‑71 (MQ=255)
                                                          ctttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCTAAta                 >  1:569464/3‑62 (MQ=255)
                                                          ctttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAAta                 >  1:4994076/3‑62 (MQ=255)
                                                          ctttcattGGTCCAACGGTATGAGAACTGGAAGGGCAAATGCCGATTTTGAAAata                        >  1:2413606/3‑56 (MQ=255)
                                                           tttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAAta                 >  1:2722853/2‑61 (MQ=255)
                                                           attcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCa        >  1:1953523/1‑71 (MQ=255)
                                                            ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaa                            <  1:975789/50‑1 (MQ=255)
                                                            ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaa                            <  1:2567559/50‑1 (MQ=255)
                                                            ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaa                            <  1:3304151/50‑1 (MQ=255)
                                                            ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                     <  1:4444206/57‑1 (MQ=255)
                                                            ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                     <  1:1480598/57‑1 (MQ=255)
                                                            ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                 <  1:5178284/61‑1 (MQ=255)
                                                            ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                 <  1:2780132/61‑1 (MQ=255)
                                                            ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                 <  1:4988320/61‑1 (MQ=255)
                                                            ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                 <  1:4937353/61‑1 (MQ=255)
                                                            ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                 <  1:4680636/61‑1 (MQ=255)
                                                            ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                 <  1:4967718/61‑1 (MQ=255)
                                                            ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCat       <  1:1549673/71‑1 (MQ=255)
                                                             tcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATc                      <  1:201504/55‑1 (MQ=255)
                                                             tcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCata      >  1:4600281/1‑71 (MQ=255)
                                                             tcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCata      >  1:185753/1‑71 (MQ=255)
                                                              cattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGaaaa                          >  1:2021963/1‑50 (MQ=255)
                                                               attGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCatat     <  1:3790774/70‑1 (MQ=255)
                                                               attGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCatat     <  1:3673656/70‑1 (MQ=255)
                                                               attGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATc    <  1:3152344/71‑1 (MQ=255)
                                                               attGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATc    <  1:2866449/71‑1 (MQ=255)
                                                                ttGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTAATCatagga   <  1:1323247/71‑4 (MQ=255)
                                                                ttGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATc    <  1:4339447/70‑1 (MQ=255)
                                                                 tGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTAATCatagga   <  1:3516884/70‑4 (MQ=255)
                                                                  ggCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAatat                       >  1:335444/1‑49 (MQ=255)
                                                                  ggCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAat                         >  1:3137799/1‑47 (MQ=255)
                                                                  ggCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCca  <  1:1537781/70‑1 (MQ=255)
                                                                    |                                                                   
TGGCTCGTCGCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCA  >  minE/2094721‑2094856

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑