Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,076,795 C→T 7.4% H121Y (CAT→TAT)  ebgC → cryptic beta‑D‑galactosidase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,076,7950CT 7.4% 256.1 / 7.4 81H121Y (CAT→TAT) ebgCcryptic beta‑D‑galactosidase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (28/46);  new base T (2/4);  total (30/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATG  >  minE/2076725‑2076864
                                                                      |                                                                     
ccGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATcc                                                                       <  1:1407652/71‑1 (MQ=255)
     aaTATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATg                                                                    >  1:658304/1‑69 (MQ=255)
      atatTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGaa                                                                  <  1:3494709/70‑1 (MQ=255)
       tatTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCAt                                                                     <  1:2617193/66‑1 (MQ=255)
             aaaGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATc                                                           <  1:1296954/70‑1 (MQ=255)
              aaGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCgg                                                         >  1:221060/1‑71 (MQ=255)
               aGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGt                                                        >  1:1490019/1‑71 (MQ=255)
                 gCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGttt                                                      <  1:4367377/71‑1 (MQ=255)
                     cGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATc                                                   <  1:155103/70‑1 (MQ=255)
                        agaAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGc                                                <  1:2915783/70‑1 (MQ=255)
                        agaAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCa                                               <  1:202113/71‑1 (MQ=255)
                         gaAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCaa                                              >  1:4079187/1‑71 (MQ=255)
                         gaAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCa                                               <  1:66286/70‑1 (MQ=255)
                          aaaCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAAt                                             <  1:4123844/71‑1 (MQ=255)
                           aaCCGTTGATGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAg                                                                 <  1:1320291/50‑1 (MQ=255)
                           aaCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAg                                                                 <  1:10642/50‑1 (MQ=255)
                           aaCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAg                                                                 <  1:242351/50‑1 (MQ=255)
                           aaCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAg                                                                 <  1:4663567/50‑1 (MQ=255)
                           aaCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAg                                                                 <  1:640330/50‑1 (MQ=255)
                           aaCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAg                                                                 <  1:1682937/50‑1 (MQ=255)
                           aaCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAg                                                                 <  1:3788417/50‑1 (MQ=255)
                           aaCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAg                                                                 <  1:444224/50‑1 (MQ=255)
                           aaCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAg                                                                 <  1:826452/50‑1 (MQ=255)
                           aaCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATc                                                           >  1:4529866/1‑56 (MQ=255)
                           aaCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATc                                                           >  1:3791314/1‑56 (MQ=255)
                            aCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATaa                                           >  1:4079891/1‑71 (MQ=255)
                            aCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATaa                                           >  1:640475/1‑71 (MQ=255)
                            aCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATaa                                           >  1:500865/1‑71 (MQ=255)
                            aCCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATaa                                           >  1:766466/1‑71 (MQ=255)
                              cggtGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAAcg                                         <  1:3939220/68‑1 (MQ=255)
                              cgatGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCTATGAAGCGTATCGGTTTATCTGCAATAAcg                                         >  1:1863360/4‑71 (MQ=255)
                              cGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCTATGAAGCGTATCGGTTTATCTGCAATAAcg                                         >  1:1432724/1‑71 (MQ=255)
                              cGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAAcg                                         <  1:3075428/71‑1 (MQ=255)
                              cGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAAcg                                         <  1:4594457/71‑1 (MQ=255)
                              cGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAAcg                                         <  1:2193699/71‑1 (MQ=255)
                                  gatgTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGt                                     >  1:110624/4‑71 (MQ=255)
                                  gAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCgg                                      <  1:1682702/70‑1 (MQ=255)
                                      tCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGt                                                              <  1:3120025/42‑1 (MQ=255)
                                      tCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGt                                                              <  1:3164316/42‑1 (MQ=255)
                                      tCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGt                                                              <  1:1129294/42‑1 (MQ=255)
                                      tCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGt                                                              <  1:1253082/42‑1 (MQ=255)
                                      tCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGt                                                              <  1:3634891/42‑1 (MQ=255)
                                      tCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaa                                  >  1:4035999/1‑70 (MQ=255)
                                       ccACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaaa                                 >  1:1389369/1‑70 (MQ=255)
                                        cACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCt                                                  >  1:75883/1‑52 (MQ=255)
                                        cACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCt                                                  >  1:4629882/1‑52 (MQ=255)
                                        cACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCt                                                  >  1:4424742/1‑52 (MQ=255)
                                        cACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaaaa                                >  1:2791001/1‑70 (MQ=255)
                                        cACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaaaa                                >  1:1291235/1‑70 (MQ=255)
                                        cACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaaaa                                >  1:1696713/1‑70 (MQ=255)
                                          cGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAg                             <  1:1948543/71‑1 (MQ=255)
                                          cGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAg                             <  1:3783815/71‑1 (MQ=255)
                                          cGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAg                             <  1:3070280/71‑1 (MQ=255)
                                                   aaTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATaa                                           >  1:3569571/1‑48 (MQ=255)
                                                   aaTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATaa                                           >  1:65749/1‑48 (MQ=255)
                                                   aaTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaaaa                                <  1:1752563/59‑1 (MQ=255)
                                                   aaTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTtctc                     >  1:684668/1‑70 (MQ=255)
                                                   aaTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTtctc                     >  1:2615191/1‑70 (MQ=255)
                                                   aaTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCa                    <  1:4598779/71‑1 (MQ=255)
                                                   aaTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCa                    <  1:4553688/71‑1 (MQ=255)
                                                   aaTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCa                    <  1:740368/71‑1 (MQ=255)
                                                   aaTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCa                    <  1:3718357/71‑1 (MQ=255)
                                                   aaTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCa                    <  1:3683430/71‑1 (MQ=255)
                                                   aaTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCa                    <  1:3082215/71‑1 (MQ=255)
                                                   aaTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCa                    <  1:3043458/71‑1 (MQ=255)
                                                   aaTCGTTATTTGCGATATCAATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCa                    <  1:2307741/71‑1 (MQ=255)
                                                    aTCGTTATTTGCGATATCTATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCaa                   <  1:1142256/71‑1 (MQ=255)
                                                    aTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaaaa                                <  1:1441456/58‑1 (MQ=255)
                                                    aTCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaaaa                                <  1:1691035/58‑1 (MQ=255)
                                                     tCGTTATTTGCGATATCCATGAAGCGTATCGGTTTa                                                     >  1:3530660/1‑36 (MQ=255)
                                                     tCGTTATTTGCGATATCCATGAAGCGTATCGGTTTa                                                     >  1:710093/1‑36 (MQ=255)
                                                      cGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaaaaaa                              >  1:1520558/1‑58 (MQ=255)
                                                       gTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGt                <  1:858489/71‑1 (MQ=255)
                                                         tATTTGCGATATCTATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaaaa                                <  1:745255/53‑1 (MQ=255)
                                                         tATTTGCGATATCTATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaaaa                                <  1:3469653/53‑1 (MQ=255)
                                                         tATTTGCGATATCTATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaaaa                                <  1:3642658/53‑1 (MQ=255)
                                                          aTTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaa                                  <  1:3380108/50‑1 (MQ=255)
                                                          aTTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaa                                  <  1:1041805/50‑1 (MQ=255)
                                                          aTTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCaa                                  <  1:3247653/50‑1 (MQ=255)
                                                            ttGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAAcg                                         <  1:4351614/41‑1 (MQ=255)
                                                                    tCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTtctc                     >  1:652034/1‑53 (MQ=255)
                                                                    tCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTtctc                     >  1:4249890/1‑53 (MQ=255)
                                                                     ccATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGaa      <  1:3388676/67‑1 (MQ=255)
                                                                     ccATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGaa      <  1:4507186/67‑1 (MQ=255)
                                                                     ccATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGaa      <  1:1140609/67‑1 (MQ=255)
                                                                     ccATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAg     >  1:1317513/1‑68 (MQ=255)
                                                                     ccATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATg  >  1:4042474/1‑71 (MQ=255)
                                                                      cATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCa              >  1:4017381/1‑58 (MQ=255)
                                                                      cATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCa              >  1:673519/1‑58 (MQ=255)
                                                                      cATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCa              >  1:3373034/1‑58 (MQ=255)
                                                                      cATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCa              >  1:171612/1‑58 (MQ=255)
                                                                      cATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATg  <  1:2733321/70‑1 (MQ=255)
                                                                      |                                                                     
CCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATATCCATGAAGCGTATCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATG  >  minE/2076725‑2076864

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: