Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 241,084 C→T 42.5% V1047M (GTG→ATG)  sbcC ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE241,0840CT42.5% 29.7 / 64.2 87V1047M (GTG→ATG) sbcCexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (22/28);  new base T (21/16);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

TAACGACCAGCAAGGCCAAGAGATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGG  >  minE/241016‑241150
                                                                    |                                                                  
tAACGACCAGCAAGGCCAAGAGATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCAc                                                                    >  1:1070228/1‑69 (MQ=255)
tAACGACCAGCAAGGCCAAGAGATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCAc                                                                    >  1:980244/1‑69 (MQ=255)
tAACGACCAGCAAGGCCAAGAGATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCAc                                                                    >  1:1161485/1‑69 (MQ=255)
tAACGACCAGCAAGGCCAAGAGATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCAc                                                                    >  1:845369/1‑69 (MQ=255)
tAACGACCAGCAAGGCCAAGAGATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCAc                                                                    >  1:3009110/1‑69 (MQ=255)
tAACGACCAGCAAGGCCAAGAGATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCAc                                                                    >  1:2963162/1‑69 (MQ=255)
tAACGACCAGCAAGGCCAAGAGATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCAc                                                                    >  1:2061174/1‑69 (MQ=255)
                  agagATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAg                                                  <  1:3252389/69‑1 (MQ=255)
                  agagATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAg                                                  <  1:2529938/69‑1 (MQ=255)
                  agagATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAg                                                  <  1:2538386/69‑1 (MQ=255)
                  agagATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAg                                                  <  1:3422236/69‑1 (MQ=255)
                  agagATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAg                                                  <  1:485944/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:1817645/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:2720103/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:479375/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:190353/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:1919259/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:3228755/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:2281758/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:3231/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:1816241/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:280635/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:2824884/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:612912/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:717226/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:1260393/69‑1 (MQ=255)
                                        ttCATTATCCTGCTGAATAGTTATTTCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTg                            <  1:2899339/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:3050543/1‑69 (MQ=255)
                                          cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                          >  1:2865706/1‑69 (MQ=255)
                                          cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                          >  1:3248476/1‑69 (MQ=255)
                                          cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                          >  1:856964/1‑69 (MQ=255)
                                          cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                          >  1:1746819/1‑69 (MQ=255)
                                          cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                          >  1:1778231/1‑69 (MQ=255)
                                          cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                          >  1:455152/1‑69 (MQ=255)
                                          cATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGAt                          >  1:2182684/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:396518/67‑1 (MQ=255)
                                                         tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg             <  1:3385839/67‑1 (MQ=255)
                                                         tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg             <  1:3007322/67‑1 (MQ=255)
                                                         tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg             <  1:2987739/67‑1 (MQ=255)
                                                         tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg             <  1:512247/67‑1 (MQ=255)
                                                         tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg             <  1:730892/67‑1 (MQ=255)
                                                         tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg             <  1:3239289/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:2895176/68‑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:2005819/68‑1 (MQ=255)
                                                         tAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGa            <  1:2200648/68‑1 (MQ=255)
                                                         tAGTTATTTAACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTg             <  1:1859314/67‑1 (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:1858116/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:2851452/1‑68 (MQ=255)
                                                                 tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc    >  1:2876781/1‑68 (MQ=255)
                                                                 tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc    >  1:89025/1‑68 (MQ=255)
                                                                 tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc    >  1:1362266/1‑68 (MQ=255)
                                                                 tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc    >  1:2929200/1‑68 (MQ=255)
                                                                 tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc    >  1:1550864/1‑68 (MQ=255)
                                                                 tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc    >  1:1609019/1‑68 (MQ=255)
                                                                 tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc    >  1:52681/1‑68 (MQ=255)
                                                                 tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc    >  1:2648132/1‑68 (MQ=255)
                                                                 tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc    >  1:3071400/1‑68 (MQ=255)
                                                                 tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc    >  1:395627/1‑68 (MQ=255)
                                                                 tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc    >  1:360481/1‑68 (MQ=255)
                                                                 tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc    >  1:3206014/1‑68 (MQ=255)
                                                                 tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc    >  1:2281905/1‑68 (MQ=255)
                                                                 tCATTGCAAACGCACTTTCCAGTTTGCTGTAGCCCAGTCCATTGATCTTCTTCACTTTGATCTGCAcc    >  1:3324311/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:2771199/68‑1 (MQ=255)
                                                                   aCTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCgg  <  1:157309/68‑1 (MQ=255)
                                                                   aCTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCgg  <  1:1158434/68‑1 (MQ=255)
                                                                    |                                                                  
TAACGACCAGCAAGGCCAAGAGATAATTCGCCCCTCTGTATTCATTATCCTGCTGAATAGTTATTTCACTGCAAACGTACTTTCCAGTTTGCTGTAGCCCAGGCCGTTGATCTTTTTCACTTTGATCTGCACCGG  >  minE/241016‑241150

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: