Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,989 T→C 7.5% intergenic (+1757/‑22) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,9890TC 7.5% 149.7 / 3.4 53intergenic (+1757/‑22)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (16/33);  new base C (2/2);  total (18/35)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAATGGCGCTTATGGTGCTAG  >  minE/2572920‑2573056
                                                                     |                                                                   
cAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACtt                                                                    >  1:497022/1‑71 (MQ=255)
       taGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAAc                                                              <  1:3280511/70‑1 (MQ=255)
        agagGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAAc                                                              >  1:1297163/1‑69 (MQ=255)
           ggCTTAAGGAAATTCGTAAGATGTCAGCCACTATACCGATATAAATAATAAGACTCACCTGCAAACCAGAc                                                         >  1:3006214/1‑71 (MQ=255)
             cTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACgg                                                       <  1:2647293/71‑1 (MQ=255)
                   gagaTTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAAtt                                                  >  1:493408/1‑70 (MQ=255)
                   gagaTTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAAtt                                                  >  1:2010998/1‑70 (MQ=255)
                    agaTTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAAc                                                              <  1:3300571/57‑1 (MQ=255)
                    agaTTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTa                                                <  1:3157943/71‑1 (MQ=255)
                     gaTTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGa                                                          <  1:3619598/60‑1 (MQ=255)
                     gaTTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGa                                                          <  1:2893990/60‑1 (MQ=255)
                        tCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAatga                                            >  1:2972118/1‑71 (MQ=255)
                          gTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAatgatg                                          <  1:619818/71‑1 (MQ=255)
                                 atCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAAcc                                                             >  1:3292384/1‑45 (MQ=255)
                                   cAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAAcc                                                             >  1:3438555/1‑43 (MQ=255)
                                    aGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGct                                 <  1:3494122/70‑1 (MQ=255)
                                    aGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGct                                 <  1:2876749/70‑1 (MQ=255)
                                    aGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGct                                 <  1:1330900/70‑1 (MQ=255)
                                     gCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGcttt                               <  1:3489300/71‑1 (MQ=255)
                                      ccACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGctttc                              >  1:2100725/1‑71 (MQ=255)
                                       cACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGctttct                             >  1:796904/1‑71 (MQ=255)
                                        aCTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGctttctt                            <  1:286569/71‑1 (MQ=255)
                                        aCTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGctttctt                            <  1:2493715/71‑1 (MQ=255)
                                          tataCCGATATAAATAATAAGACTCACTTGCAAACCa                                                            <  1:2606018/37‑1 (MQ=255)
                                          tataCCGATATAAACAATAAGACTCACCTGCAAACCAGACGGTAATTTAATGATGATGAACGctttcttt                           <  1:1235599/70‑1 (MQ=255)
                                            taCCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACg                                   <  1:2604146/60‑1 (MQ=255)
                                             aCCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCgg                       >  1:3105408/1‑71 (MQ=255)
                                             aCCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCgg                       >  1:855299/1‑71 (MQ=255)
                                              ccGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCgg                       <  1:2597802/70‑1 (MQ=255)
                                              ccGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGc                      <  1:1770654/71‑1 (MQ=255)
                                               cGATATAAATAATAAGACTCACTTGCAAACCAGACGGTa                                                     >  1:23755/1‑39 (MQ=255)
                                                gATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCaa                    <  1:2366424/71‑1 (MQ=255)
                                                 atatAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAAt                   >  1:585878/1‑71 (MQ=255)
                                                 atatAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAAt                   >  1:1117127/1‑71 (MQ=255)
                                                    taaataaTAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGc                      <  1:1297051/65‑1 (MQ=255)
                                                    taaataaTAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGc                      <  1:2498539/65‑1 (MQ=255)
                                                    taaataaTAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAATgg                 <  1:3401984/70‑1 (MQ=255)
                                                    taaataaTAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAATgg                 <  1:2258209/70‑1 (MQ=255)
                                                      aataatAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGctttctt                            <  1:2322626/57‑1 (MQ=255)
                                                         aataaGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGc                      <  1:2806604/60‑1 (MQ=255)
                                                         aataaGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGc                      <  1:667705/60‑1 (MQ=255)
                                                          ataaGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAATGGCGCTTa           <  1:2540605/70‑1 (MQ=255)
                                                           taaGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAAc                                    <  1:1331837/44‑1 (MQ=255)
                                                           taaGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAATGGCGCTTa           <  1:2497575/69‑1 (MQ=255)
                                                           taaGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAATGGCGCTTATg         <  1:3523933/71‑1 (MQ=255)
                                                           taaGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAATGGCGCTTATg         <  1:3603025/71‑1 (MQ=255)
                                                           taaGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAATGGCGCTTATg         >  1:3577324/1‑71 (MQ=255)
                                                           taaGACTCACCTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAATGGCGCTTATg         >  1:2314387/1‑71 (MQ=255)
                                                              gACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAATGGCGCTTATGGTg      >  1:109426/1‑71 (MQ=255)
                                                                 tCACCTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAATGGCGCt             <  1:2604881/61‑1 (MQ=255)
                                                                   aCTTGCAAACCAGACGGTAATTTAATGATGATGAACGctt                                >  1:3498164/1‑40 (MQ=255)
                                                                   aCTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAATGGCGCTTATGGTGCTAg  <  1:2251272/70‑1 (MQ=255)
                                                                   aCTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAATGGCGCTTATGGTGCTAg  <  1:1642018/70‑1 (MQ=255)
                                                                    cTTGCAAACCAGACGGTAATTTAATGATGATGAACGCtttctttc                          <  1:873274/45‑1 (MQ=255)
                                                                     ttGCAAACCAGACGGTAATTTAATGATGATGAACGctttcttt                           >  1:309388/1‑43 (MQ=255)
                                                                     ttGCAAACCAGACGGTAATTTAATGATGATGAACGctttcttt                           >  1:3310740/1‑43 (MQ=255)
                                                                     ttGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGc                      <  1:2611738/48‑1 (MQ=255)
                                                                     |                                                                   
CAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTTAATGATGATGAACGCTTTCTTTCCGGCAATGGCGCTTATGGTGCTAG  >  minE/2572920‑2573056

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: