Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 241,121 G→A 37.6% N1034N (AAC→AAT sbcC ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE241,1210GA37.6% 74.8 / 77.5 122N1034N (AAC→AATsbcCexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (41/35);  new base A (30/16);  total (71/51)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 122 total aligned reads displayed.

TTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCTTCTACGTGGCTAAT  >  minE/241056‑241183
                                                                 |                                                              
ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                                             <  1:2720103/69‑1 (MQ=255)
ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                                             <  1:280635/69‑1 (MQ=255)
ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                                             <  1:2899339/69‑1 (MQ=255)
ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                                             <  1:2281758/69‑1 (MQ=255)
ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                                             <  1:717226/69‑1 (MQ=255)
ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                                             <  1:1260393/69‑1 (MQ=255)
ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                                             <  1:1817645/69‑1 (MQ=255)
ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                                             <  1:3231/69‑1 (MQ=255)
ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                                             <  1:479375/69‑1 (MQ=255)
ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                                             <  1:1919259/69‑1 (MQ=255)
ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                                             <  1:190353/69‑1 (MQ=255)
ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                                             <  1:3228755/69‑1 (MQ=255)
ttCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTg                                                             >  1:1781231/1‑69 (MQ=255)
 tCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTATGCTGTAGCCCAGTCCATTg                                                             <  1:514840/68‑1 (MQ=255)
  cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGGCGTTGAt                                                           >  1:1408251/1‑69 (MQ=255)
  cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGa                                                            >  1:1629080/1‑68 (MQ=255)
  cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                                           >  1:1020227/1‑69 (MQ=255)
  cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                                           >  1:2182684/1‑69 (MQ=255)
  cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                                           >  1:2865706/1‑69 (MQ=255)
  cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                                           >  1:1778231/1‑69 (MQ=255)
  cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                                           >  1:3248476/1‑69 (MQ=255)
  cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                                           >  1:856964/1‑69 (MQ=255)
  cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                                           >  1:455152/1‑69 (MQ=255)
  cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCACAGGCCGTTGAt                                                           >  1:3376484/1‑69 (MQ=255)
   aTTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATc                                                          >  1:3125601/1‑69 (MQ=255)
                 tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                                              <  1:3007322/67‑1 (MQ=255)
                 tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                                              <  1:2567791/67‑1 (MQ=255)
                 tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                                              <  1:2797596/67‑1 (MQ=255)
                 tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                                              <  1:730892/67‑1 (MQ=255)
                 tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                                              <  1:512247/67‑1 (MQ=255)
                 tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                                              <  1:2101935/67‑1 (MQ=255)
                 tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                                              <  1:2987739/67‑1 (MQ=255)
                 tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                                              <  1:396518/67‑1 (MQ=255)
                 tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                                              <  1:3239289/67‑1 (MQ=255)
                 tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGa                                             <  1:2005819/68‑1 (MQ=255)
                 tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGa                                             <  1:3429108/68‑1 (MQ=255)
                 tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGa                                             <  1:779620/68‑1 (MQ=255)
                 tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGa                                             <  1:2895176/68‑1 (MQ=255)
                 tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGa                                             <  1:2458974/68‑1 (MQ=255)
                      aTTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCa                                       >  1:1765779/1‑69 (MQ=255)
                      aTTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCa                                       >  1:1934868/1‑69 (MQ=255)
                      aTTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCa                                       >  1:1858116/1‑69 (MQ=255)
                      aTTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCa                                       >  1:1033669/1‑69 (MQ=255)
                      aTTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCa                                       >  1:1759204/1‑69 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACtt                                                >  1:3314314/1‑57 (MQ=37)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGAt                                            >  1:233559/1‑61 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCt                                          >  1:92309/1‑63 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc                                     >  1:3324311/1‑68 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc                                     >  1:1550864/1‑68 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc                                     >  1:395627/1‑68 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc                                     >  1:1362266/1‑68 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc                                     >  1:3071400/1‑68 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc                                     >  1:52681/1‑68 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc                                     >  1:360481/1‑68 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc                                     >  1:2929200/1‑68 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc                                     >  1:2281905/1‑68 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc                                     >  1:2876781/1‑68 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc                                     >  1:1609019/1‑68 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc                                     >  1:2851452/1‑68 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc                                     >  1:2648132/1‑68 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc                                     >  1:89025/1‑68 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAc                                      >  1:1095361/1‑67 (MQ=255)
                         tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCAAGTCCATTGATCTTCTTCACTTTGa                                             >  1:2030315/1‑60 (MQ=255)
                         tCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCAcc                                     >  1:1640536/1‑68 (MQ=255)
                           aCTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCgg                                   <  1:1158434/68‑1 (MQ=255)
                           aCTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCgg                                   <  1:2771199/68‑1 (MQ=255)
                             tGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                                              <  1:1272222/55‑1 (MQ=255)
                               cAAACGTACTTTCCAGTTTGCTGTAGCCCAGGTCGTTGATCtt                                                        >  1:3042714/1‑43 (MQ=255)
                               cAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAAt                                >  1:655197/1‑67 (MQ=255)
                               cAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATAc                              >  1:2584604/1‑69 (MQ=255)
                               cAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATAc                              >  1:498491/1‑69 (MQ=255)
                               cAAACGTACTTTCCAGTTTGCTGTAGACCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATAc                              >  1:574309/1‑69 (MQ=255)
                                    gTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGctct                         <  1:1849505/69‑1 (MQ=255)
                                    gTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGctct                         <  1:2007360/69‑1 (MQ=255)
                                       cTTTCCAGTTTGCTGTACCCCAGGCCGTTGATCTTTTTCACttt                                               >  1:1715077/1‑44 (MQ=255)
                                         ttCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGc                                        >  1:1533721/1‑49 (MQ=255)
                                         ttCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGc                                        >  1:63468/1‑49 (MQ=255)
                                         ttCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGc                                        >  1:2602322/1‑49 (MQ=255)
                                         ttCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGc                                        >  1:3327722/1‑49 (MQ=255)
                                         ttCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGc                                        >  1:1844986/1‑49 (MQ=255)
                                         ttCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGc                                        >  1:2390377/1‑49 (MQ=255)
                                         ttCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGc                                        >  1:2346315/1‑49 (MQ=255)
                                         ttCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGc                                        >  1:3114638/1‑49 (MQ=255)
                                         ttCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGc                                        >  1:306933/1‑49 (MQ=255)
                                                 tGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCTTCTa            <  1:2291991/69‑1 (MQ=255)
                                                    tGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCTTCTACg          <  1:1973234/68‑1 (MQ=255)
                                                      tAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGcct               >  1:2312918/1‑59 (MQ=255)
                                                      tAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCCTc              >  1:1009058/1‑62 (MQ=255)
                                                      tAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCCTc              >  1:2608184/1‑62 (MQ=255)
                                                      tAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCCTc              >  1:2475765/1‑62 (MQ=255)
                                                      tAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCCTc              >  1:3213541/1‑62 (MQ=255)
                                                      tAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCCTc              >  1:1656509/1‑62 (MQ=255)
                                                      tAGCCCAGTCCATTGATCTTCCTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCCTc              >  1:3086990/1‑62 (MQ=255)
                                                        gCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCTTCTACGTGGCt     >  1:45251/1‑69 (MQ=255)
                                                           cAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCTTCTACGTGGCTAAt  <  1:149525/69‑1 (MQ=255)
                                                           cAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCTTCTACGTGGCTAAt  <  1:654519/69‑1 (MQ=255)
                                                           cAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCTTCTACGTGGCTAAt  <  1:477613/69‑1 (MQ=255)
                                                           cAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCTTCTACGTGGCTAAt  <  1:1242142/69‑1 (MQ=255)
                                                           cAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCTTCTACGTGGCTAAt  <  1:105831/69‑1 (MQ=255)
                                                           cAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCTTCTACGTGGCTAAt  <  1:1818973/69‑1 (MQ=255)
                                                                 |                                                              
TTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGCTCTTTCATCGCTTCTACGTGGCTAAT  >  minE/241056‑241183

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: