Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,977,661 T→C 37.3% A107A (GCT→GCC endA → DNA‑specific endonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,977,6610TC37.3% 102.4 / 114.5 174A107A (GCT→GCCendADNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base T (48/61);  new base C (12/53);  total (60/114)
Fisher's exact test for biased strand distribution p-value = 5.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 174 total aligned reads displayed.

TCCCGCCTGGCAGTTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCGGTGAG  >  minE/1977598‑1977724
                                                               |                                                               
tCCTGCCTGGCAATTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAg                                                             <  1:2190982/68‑1 (MQ=255)
tCCTGCCTGGCAATTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAg                                                             <  1:1016460/68‑1 (MQ=255)
tCCTGCCTGGCAATTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAg                                                             <  1:2253654/68‑1 (MQ=255)
tCCTGCCTGGCAATTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAg                                                             <  1:1299929/68‑1 (MQ=255)
tCCTGCCTGGCAATTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAg                                                             <  1:212477/68‑1 (MQ=255)
tCCCGCCTGGCAGTTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAg                                                             <  1:3360546/68‑1 (MQ=255)
 cctgCCTGGCAATTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGAt                                                           >  1:3279709/4‑69 (MQ=255)
 cctgCCTGGCAATTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGAt                                                           >  1:393844/4‑69 (MQ=255)
 cctgCCTGGCAATTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGAt                                                           >  1:2275601/4‑69 (MQ=255)
 cctgCCTGGCAATTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGAt                                                           >  1:1347218/4‑69 (MQ=255)
 cctgCCTGGCAATTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGAt                                                           >  1:234160/4‑69 (MQ=255)
 cctgCCTGGCAATTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGAt                                                           >  1:2287174/4‑69 (MQ=255)
 cccGCCTGGCAGTTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGAt                                                           >  1:2562089/1‑69 (MQ=255)
 cccGCCTGGCAGTTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGAt                                                           >  1:3075075/1‑69 (MQ=255)
 cccGCCTGGCAGTTCGGTCACCAGCGCCAGTGCTGGCAGGAGGTGGACGTAAAAACTGCGCTAAAGATc                                                          >  1:56857/1‑69 (MQ=255)
              tCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGATc                                                          <  1:170920/57‑1 (MQ=255)
              tCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGATc                                                          <  1:3110976/57‑1 (MQ=255)
              tCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGATc                                                          <  1:2502014/57‑1 (MQ=255)
              tCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGATc                                                          <  1:2621409/57‑1 (MQ=255)
              tCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGATc                                                          <  1:2752223/57‑1 (MQ=255)
              tCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGATc                                                          <  1:777226/57‑1 (MQ=255)
              tCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGATc                                                          <  1:1385774/57‑1 (MQ=255)
              tCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGATc                                                          <  1:1247536/57‑1 (MQ=255)
               cGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTAt                                                  >  1:740173/1‑64 (MQ=255)
               cGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATc                                                 >  1:2986777/1‑65 (MQ=255)
               cGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATc                                                 >  1:1845442/1‑65 (MQ=255)
               cGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATc                                                 >  1:1948746/1‑65 (MQ=255)
               cGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATc                                                 >  1:2641147/1‑65 (MQ=255)
               cGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCAATc                                                 >  1:2799235/1‑65 (MQ=255)
                      cAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGaa                                      <  1:3176150/69‑1 (MQ=255)
                      cAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGaa                                      >  1:2894710/1‑69 (MQ=255)
                      cAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGaa                                      >  1:833108/1‑69 (MQ=255)
                      cAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGa                                       <  1:1767217/68‑1 (MQ=255)
                      cAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGa                                       <  1:2462552/68‑1 (MQ=255)
                       aGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAATTGCGCCAAAGATc                                                          <  1:2202912/48‑1 (MQ=39)
                           ccAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGa                                 <  1:925781/69‑1 (MQ=255)
                           ccAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGa                                 <  1:2504773/69‑1 (MQ=255)
                           ccAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGa                                 <  1:1649485/69‑1 (MQ=255)
                           ccAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGa                                 <  1:966673/69‑1 (MQ=255)
                           ccAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGa                                 <  1:1612419/69‑1 (MQ=255)
                            cAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGat                                <  1:1362580/69‑1 (MQ=255)
                                gCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAg                                                             >  1:1726278/1‑36 (MQ=255)
                                 cTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCt                                                    >  1:1485080/1‑44 (MQ=255)
                                 cTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGc                            >  1:64577/1‑68 (MQ=255)
                                 cTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGc                            >  1:2955833/1‑68 (MQ=255)
                                 cTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGc                            >  1:199950/1‑68 (MQ=255)
                                 cTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGc                            >  1:1142324/1‑68 (MQ=255)
                                   ggCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATa                         <  1:1386682/69‑1 (MQ=255)
                                   ggCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATa                         <  1:3363280/69‑1 (MQ=255)
                                   ggCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATa                         <  1:151785/69‑1 (MQ=255)
                                    gCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATa                         <  1:2856732/68‑1 (MQ=255)
                                          gggtggtCGTAAAATCTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCgt                                 <  1:710054/53‑2 (MQ=255)
                                           ggTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAATCGa                                 >  1:1451849/1‑53 (MQ=255)
                                           ggTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGa                                 <  1:3425322/53‑1 (MQ=255)
                                           ggTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGa                                 <  1:690662/53‑1 (MQ=255)
                                           ggTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGa                                 <  1:2282314/53‑1 (MQ=255)
                                           ggTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGa                                 <  1:3262266/53‑1 (MQ=255)
                                           ggTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGa                                 <  1:2853439/53‑1 (MQ=255)
                                           ggTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGa                                 >  1:352012/1‑53 (MQ=255)
                                           ggTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGa                                 <  1:54876/53‑1 (MQ=255)
                                             tGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGa                                 <  1:374649/51‑1 (MQ=255)
                                                 cGTAAAAATTGCGCCAAAGATCCGGTCTATCGCAAGATGGAAAGCGatat                              <  1:577468/50‑1 (MQ=38)
                                                 cGTAAAAATTGCGCCAAAGATCCGGTCTATCGCAAGATGGAAAGCGatat                              <  1:943526/50‑1 (MQ=38)
                                                 cGTAAAAATTGCGCCAAAGATCCGGTCTATCGCAAGATGGAAAGCGatat                              <  1:3240102/50‑1 (MQ=38)
                                                 cGTAAAAATTGCGCCAAAGATCCGGTCTATCGCAAGATGGAAAGCGatat                              <  1:797179/50‑1 (MQ=38)
                                                 cGTAAAAATTGCGCCAAAGATCCGGTCTATCGCAAGATGGAAAGCGatat                              <  1:53134/50‑1 (MQ=38)
                                                 cGTAAAAATTGCGCCAAAGATCCGGTCTATCGCAAGATGGAAAGCGatat                              <  1:2669695/50‑1 (MQ=38)
                                                 cGTAAAAATTGCGCCAAAGATCCGGTCTATCGCAAGATGGAAAGCGatat                              <  1:2971750/50‑1 (MQ=38)
                                                 cGTAAAAATTGCGCCAAAGATCCGGTCTATCGCAAGATGGAAAGCGatat                              <  1:1403409/50‑1 (MQ=38)
                                                 cGTAAAAATTGCGCCAAAGATCCGGTCTATCGCAAGATGGAAAGCGatat                              <  1:1677092/50‑1 (MQ=38)
                                                 cGTAAAAATTGCGCCAAAGATCCGGTCTATCGCAAGATGGAAAGCGatat                              <  1:1889839/50‑1 (MQ=38)
                                                 cGTAAAAATTGCGCCAAAGATCCGGTCTATCGCAAGATGGAAAGCGatat                              <  1:2061628/50‑1 (MQ=38)
                                                  gTAAAAATTGCGCCAAAGATCCGGTCTATCGCAAGATGGAAAGCGatat                              <  1:2775780/49‑1 (MQ=38)
                                                     aaaaCTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATa                         <  1:3157555/51‑1 (MQ=255)
                                                     aaaaCTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATa                         <  1:324234/51‑1 (MQ=255)
                                                     aaaaCTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATa                         <  1:1056137/51‑1 (MQ=255)
                                                     aaaaCTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATa                         <  1:2773752/51‑1 (MQ=255)
                                                     aaaaCTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATa                         <  1:2761177/51‑1 (MQ=255)
                                                     aaaaCTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATa                         <  1:3371602/51‑1 (MQ=255)
                                                     aaaaCTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCg       >  1:1027730/1‑69 (MQ=255)
                                                     aaaaCTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCg       >  1:149088/1‑69 (MQ=255)
                                                     aaaaCTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCg       >  1:2008063/1‑69 (MQ=255)
                                                     aaaaCTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCg       >  1:2219334/1‑69 (MQ=255)
                                                        aCTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAAcc                      >  1:2557516/1‑51 (MQ=255)
                                                         cTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTc        >  1:338694/1‑64 (MQ=255)
                                                         cTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCggtga   >  1:1506857/1‑69 (MQ=255)
                                                         cTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCggtga   <  1:873439/69‑1 (MQ=255)
                                                         cTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCggtga   >  1:1082361/1‑69 (MQ=255)
                                                         cTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCggtga   <  1:1303868/69‑1 (MQ=255)
                                                         cTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCggtga   <  1:3405892/69‑1 (MQ=255)
                                                         cTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCggtga   >  1:46979/1‑69 (MQ=255)
                                                         cTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCggtga   <  1:1622531/69‑1 (MQ=255)
                                                         cTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCggtga   <  1:2685442/69‑1 (MQ=255)
                                                         cTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCggtga   <  1:2004816/69‑1 (MQ=255)
                                                         cTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCggtga   <  1:551398/69‑1 (MQ=255)
                                                         cTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCggtga   <  1:2446091/69‑1 (MQ=255)
                                                          tGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATa                         >  1:640856/1‑46 (MQ=255)
                                                          tGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCggtgag  >  1:1068386/1‑69 (MQ=255)
                                                          tGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCggtgag  >  1:2907667/1‑69 (MQ=255)
                                                          tGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCggtgag  >  1:3033827/1‑69 (MQ=255)
                                                               |                                                               
TCCCGCCTGGCAGTTCGGTCACCAGCGCCAGTGCTGGCAGGACGGTGGACGTAAAAACTGCGCTAAAGATCCGGTCTATCGCAAGATGGAAAGCGATATGCATAACCTGCAGCCGTCAGTCGGTGAG  >  minE/1977598‑1977724

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: