Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 241,093 T→C 38.1% T1044A (ACG→GCG)  sbcC ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE241,0930TC38.1% 48.0 / 65.0 97T1044A (ACG→GCG) sbcCexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base C (21/16);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGAGATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGCTCT  >  minE/241034‑241160
                                                           |                                                                   
agagATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAg                                                            <  1:3252389/69‑1 (MQ=255)
agagATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAg                                                            <  1:2529938/69‑1 (MQ=255)
agagATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAg                                                            <  1:485944/69‑1 (MQ=255)
agagATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAg                                                            <  1:3422236/69‑1 (MQ=255)
agagATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAg                                                            <  1:2538386/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:2720103/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:2899339/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:2281758/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:3228755/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:3231/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:2824884/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:1919259/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:190353/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:1817645/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:1816241/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:479375/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:280635/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:612912/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:1260393/69‑1 (MQ=255)
                      ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                                      <  1:717226/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:455152/1‑69 (MQ=255)
                        cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                    >  1:2865706/1‑69 (MQ=255)
                        cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                    >  1:856964/1‑69 (MQ=255)
                        cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                    >  1:3050543/1‑69 (MQ=255)
                        cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                    >  1:2182684/1‑69 (MQ=255)
                        cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                    >  1:1020227/1‑69 (MQ=255)
                        cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                    >  1:3248476/1‑69 (MQ=255)
                        cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                    >  1:1778231/1‑69 (MQ=255)
                        cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                                    >  1:1746819/1‑69 (MQ=255)
                        cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCACAGGCCGTTGAt                                    >  1:3376484/1‑69 (MQ=255)
                         aTTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATc                                   >  1:3125601/1‑69 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACtttt                       <  1:2615244/67‑2 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                       <  1:2987739/67‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                       <  1:3007322/67‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                       <  1:3239289/67‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                       <  1:3385839/67‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                       <  1:396518/67‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                       <  1:512247/67‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                       <  1:730892/67‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                       <  1:999853/67‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                       <  1:2797596/67‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                       <  1:2567791/67‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                       <  1:2318164/67‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                       <  1:2101935/67‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGa                      <  1:779620/68‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGa                      <  1:3429108/68‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGa                      <  1:2458974/68‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGa                      <  1:2895176/68‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGa                      <  1:2005819/68‑1 (MQ=255)
                                       tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGa                      <  1:2200648/68‑1 (MQ=255)
                                       tAGTTATTTAACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg                       <  1:1859314/67‑1 (MQ=255)
                                            aTTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCa                >  1:1858116/1‑69 (MQ=255)
                                            aTTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCa                >  1:1033669/1‑69 (MQ=255)
                                            aTTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCa                >  1:1759204/1‑69 (MQ=255)
                                            aTTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCa                >  1:1765779/1‑69 (MQ=255)
                                            aTTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCa                >  1:1934868/1‑69 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTccat                                        >  1:934760/1‑40 (MQ=25)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACtt                         >  1:3314314/1‑57 (MQ=37)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGAt                     >  1:233559/1‑61 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCt                   >  1:92309/1‑63 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:2648132/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:1550864/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:1609019/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:1362266/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:52681/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:395627/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:360481/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:89025/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:2851452/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:3324311/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:2876781/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:2929200/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:3206014/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:2281905/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc              >  1:3071400/1‑68 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAc               >  1:1095361/1‑67 (MQ=255)
                                               tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCAAGTCCATTGATCTTCTTCACTTTGa                      >  1:2030315/1‑60 (MQ=255)
                                               tCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAg                                             >  1:259199/1‑37 (MQ=255)
                                               tCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCAcc              >  1:1640536/1‑68 (MQ=255)
                                               tCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCAcc              >  1:432273/1‑68 (MQ=255)
                                                 aCTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCgg            <  1:157309/68‑1 (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:498491/1‑69 (MQ=255)
                                                     cAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATAc       >  1:2584604/1‑69 (MQ=255)
                                                     cAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATAc       >  1:2494677/1‑69 (MQ=255)
                                                     cAAACGTACTTTCCAGTTTGCTGTAGACCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATAc       >  1:574309/1‑69 (MQ=255)
                                                          gTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGctct  <  1:1154965/69‑1 (MQ=255)
                                                          gTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGctct  <  1:1849505/69‑1 (MQ=255)
                                                          gTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGctct  <  1:2007360/69‑1 (MQ=255)
                                                           |                                                                   
AGAGATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGGAATACGCTCT  >  minE/241034‑241160

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: