Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,094,567 T→C 33.3% I604V (ATT→GTT)  ydiJ ← predicted FAD‑linked oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,094,5670TC33.3% 62.4 / 52.4 57I604V (ATT→GTT) ydiJpredicted FAD‑linked oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (20/18);  new base C (10/9);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GCCGCGATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGT  >  minE/1094502‑1094635
                                                                 |                                                                    
gccgcGATCCGCCAACAAACGCAGCCATTCACGCCCCAGCGTTGCGCGCCCTTTCGGTGAATGAACCCGGt                                                                 >  1:3110066/1‑71 (MQ=255)
gccgcGATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAACCCGGt                                                                 >  1:3010894/1‑71 (MQ=255)
gccgcGATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAACCCGGt                                                                 >  1:3442040/1‑71 (MQ=255)
  cgcgATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTc                                                               <  1:3589212/71‑1 (MQ=255)
  cgcgATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTc                                                               <  1:2244312/71‑1 (MQ=255)
  cgcgATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTc                                                               <  1:3794971/71‑1 (MQ=255)
  cgcgATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTc                                                               <  1:2121667/71‑1 (MQ=255)
        ccGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAACCCGGTTCTGGGt                                                          >  1:1474944/1‑70 (MQ=255)
         cGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAACCCGGTTCTGGGTg                                                         <  1:1300538/70‑1 (MQ=255)
         cGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAACCCGGTTCTGGGTg                                                         <  1:3778496/70‑1 (MQ=255)
                aaaCGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAACCCGGTTCTGGGTGATCTTCAt                                                 >  1:2952103/1‑71 (MQ=255)
                aaaCGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAACCCGGTTCTGGGTGATCTTCAt                                                 >  1:2285726/1‑71 (MQ=255)
                aaaCGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAACCCGGTTCTGGGTGATCTTCAt                                                 >  1:2143816/1‑71 (MQ=255)
                aaaCGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAACCCGGTTCTGGGTGATCTTCAt                                                 >  1:1487982/1‑71 (MQ=255)
                 aaCGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATc                                                <  1:769142/71‑1 (MQ=255)
                 aaCGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATc                                                <  1:1838316/71‑1 (MQ=255)
                    gCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAACCCGGTTCTGGGTg                                                         <  1:3124679/59‑1 (MQ=255)
                       gCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGa                                          <  1:3821289/71‑1 (MQ=255)
                       gCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGa                                          <  1:4093009/71‑1 (MQ=255)
                         cATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGaca                                        <  1:3592776/71‑1 (MQ=255)
                            tCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACAt                                     >  1:4013628/1‑71 (MQ=255)
                            tCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACAt                                     >  1:2711837/1‑71 (MQ=255)
                             cACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATc                                    >  1:3791646/1‑71 (MQ=255)
                             cACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATc                                    >  1:2252052/1‑71 (MQ=255)
                             cACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATc                                    >  1:2748763/1‑71 (MQ=255)
                              aCGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCg                                   <  1:1640670/71‑1 (MQ=255)
                              aCGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCg                                   <  1:1044575/71‑1 (MQ=255)
                              aCGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAACCCGGTTCTGGGTGATCTTCATCGACGGACACATCg                                   <  1:2361188/71‑1 (MQ=255)
                                gCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCgg                                  <  1:1892313/70‑1 (MQ=255)
                                    cAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTAc                             <  1:995342/71‑1 (MQ=255)
                                    cAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTAc                             <  1:881642/71‑1 (MQ=255)
                                    cAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTAc                             <  1:386748/71‑1 (MQ=255)
                                     aGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACg                            >  1:547268/1‑71 (MQ=255)
                                     aGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACg                            >  1:3070344/1‑71 (MQ=255)
                                     aGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACg                            >  1:4098458/1‑71 (MQ=255)
                                      gCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTc                                                   >  1:47764/1‑47 (MQ=255)
                                          tGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCAt                       >  1:1403616/1‑71 (MQ=255)
                                           gcgcgcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCgg                                  >  1:2279447/1‑59 (MQ=255)
                                           gcgcgcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGa                                 >  1:3682355/1‑60 (MQ=255)
                                           gcgcgcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGa                                 >  1:2473337/1‑60 (MQ=255)
                                           gcgcgcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGa                                 >  1:2200209/1‑60 (MQ=255)
                                            cgcgcCCTTTCTGTGAATGAACCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCa                     <  1:3936564/71‑1 (MQ=255)
                                            cgcgcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATc                      <  1:2543157/70‑1 (MQ=255)
                                            cgcgcCCTTTCGGTGAATGAACCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCa                     <  1:1350096/71‑1 (MQ=255)
                                            cgcgcCCTTTCGGTGAATGAACCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCa                     <  1:893668/71‑1 (MQ=255)
                                              cgcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCaaa                   <  1:693901/71‑1 (MQ=255)
                                              cgcCCTTTCGGTGAATGAACCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCAtaa                     >  1:1624598/1‑67 (MQ=255)
                                                  cTTTCGGTGAATGAACCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTg               <  1:3064006/71‑1 (MQ=255)
                                                  cTTTCGGTGAATGAACCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTg               <  1:586221/71‑1 (MQ=255)
                                                     tCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTa                              >  1:3630044/1‑53 (MQ=255)
                                                     tCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGaa             >  1:1157582/1‑70 (MQ=255)
                                                     tCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAg            >  1:390214/1‑71 (MQ=255)
                                                     tCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAg            >  1:3425914/1‑71 (MQ=255)
                                                     tCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAg            >  1:293875/1‑71 (MQ=255)
                                                     tCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAg            >  1:2122327/1‑71 (MQ=255)
                                                     tCGGTGAATGAACCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAg            >  1:108809/1‑71 (MQ=255)
                                                       ggTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGc           <  1:3390226/70‑1 (MQ=255)
                                                           aatgaatCCGGTTCTGGGTGATCTTCATCGACGGACACATCg                                   <  1:4125754/42‑1 (MQ=255)
                                                           aatgaatCCGGTTCTGGGTGAGCTTCATCGACGGACACATCg                                   <  1:385002/42‑1 (MQ=255)
                                                               aatCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGt  <  1:2370766/71‑1 (MQ=255)
                                                               aatCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGt  <  1:2126987/71‑1 (MQ=255)
                                                               aatCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGt  <  1:1991017/71‑1 (MQ=255)
                                                                 |                                                                    
GCCGCGATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGT  >  minE/1094502‑1094635

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