Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,520,434 G→A W44* (TGG→TAG)  ligB → DNA ligase, NAD(+)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,520,4340GA98.3% 193.8 / ‑3.9 60W44* (TGG→TAG) ligBDNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (34/25);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TCCGGCCTGGTCGCCAGCCAGAGCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTGGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGGTGTTTACGATCAGTTAAGCGC  >  minE/2520369‑2520500
                                                                 |                                                                  
tCCGGCCTGGTCGCCAGCCAGAGCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTagg                                                                  <  1:1264625/68‑1 (MQ=255)
tCCGGCCTGGTCGCCAGCCAGAGCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGAc                                                                <  1:322581/70‑1 (MQ=255)
 ccGGCCTGGTCGCCAGCCAGAGCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGa                                                              >  1:622422/1‑71 (MQ=255)
    gCCTGGTCGCCAGCCAGAGCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGa                                                           >  1:1847919/1‑71 (MQ=255)
      cTGGTCGCCAGCCAGAGCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGac                                                          >  1:795795/1‑70 (MQ=255)
         gTCGCCAGCCAGAGCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGactact                                                      >  1:885458/1‑71 (MQ=255)
         gTCGCCAGCCAGAGCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGactact                                                      >  1:1235886/1‑71 (MQ=255)
             ccagccagAGCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTgg                                                    >  1:410139/1‑69 (MQ=255)
              cagccagAGCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTggaag                                                 <  1:101218/71‑1 (MQ=255)
               agccagAGCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTggaagg                                                <  1:469906/71‑1 (MQ=255)
               agccagAGCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTggaag                                                 >  1:252948/1‑70 (MQ=255)
                gccagAGCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTggaagga                                               <  1:133031/71‑1 (MQ=255)
                 ccagAGCAGAGGAAGAAATTTCCCGCCTGCAACAGCAAATTAAACAGTAGGACGAt                                                             <  1:754097/56‑1 (MQ=255)
                   agagCACAGGAAGAAATTTCCCGCCTGCAACCGCAAATAAAACAGTAGGACGATGACTACTggaaggaag                                             >  1:1879795/1‑70 (MQ=255)
                    gagCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTggaaggaagg                                            >  1:1102076/1‑70 (MQ=255)
                     agCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTggaaggaagga                                           >  1:1234877/1‑70 (MQ=255)
                     agCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTggaaggaagga                                           >  1:2313448/1‑70 (MQ=255)
                       cacaGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGaaa                                         >  1:794305/1‑70 (MQ=255)
                       cacaGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGaaa                                         >  1:1547632/1‑70 (MQ=255)
                         caGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGa                                                           >  1:996969/1‑50 (MQ=255)
                         caGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAg                                      >  1:1656567/1‑71 (MQ=255)
                          aGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGaaa                                         <  1:306404/67‑1 (MQ=255)
                           ggAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGaaa                                         >  1:1854763/1‑66 (MQ=255)
                           ggAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTg                                    <  1:376972/71‑1 (MQ=255)
                             aagaaATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAg                                  <  1:1012380/71‑1 (MQ=255)
                             aagaaATTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAg                                  <  1:1490611/71‑1 (MQ=255)
                              agaaATTTCCCGCCTGCAACAGCAAATAAAACAGTGGGACGATGACTACTGGAAGTAAGGAAAAAGTGAgg                                 <  1:1539475/71‑1 (MQ=255)
                                 aaTTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTgg                              <  1:1607/71‑1 (MQ=255)
                                 aaTTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTg                               >  1:1387778/1‑70 (MQ=255)
                                  aTTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGactac                                                       <  1:2208185/45‑1 (MQ=255)
                                  aTTTCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTgg                              >  1:873439/1‑70 (MQ=255)
                                   tttCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTgg                              <  1:1530434/69‑1 (MQ=255)
                                   tttCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTgg                              <  1:484629/69‑1 (MQ=255)
                                    ttCCCGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTgg                              >  1:1530094/1‑68 (MQ=255)
                                      cccGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGa                          <  1:1608637/70‑1 (MQ=255)
                                      cccGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGa                          <  1:584296/70‑1 (MQ=255)
                                      cccGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGAc                         >  1:1193822/1‑71 (MQ=255)
                                       ccGCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACg                        <  1:1258603/71‑1 (MQ=255)
                                         gCCTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGgt                      >  1:2011744/1‑71 (MQ=255)
                                          ccTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGgtg                     >  1:1133452/1‑71 (MQ=255)
                                          ccTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGgtg                     >  1:1404394/1‑71 (MQ=255)
                                           cTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGaaaaa                                       >  1:1368736/1‑52 (MQ=255)
                                           cTGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGgtgt                    >  1:1665582/1‑71 (MQ=255)
                                           cTGCAACAGCAAATAAAACAGTAGGACGATGACAACTGGAAGGAAGGaaaa                                        >  1:26616/1‑51 (MQ=255)
                                            tGCAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGgtg                     >  1:402087/1‑69 (MQ=255)
                                              cAACAGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGGTGtt                   <  1:662954/69‑1 (MQ=255)
                                                  aGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGa                          <  1:1362362/58‑1 (MQ=255)
                                                  aGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGa                          <  1:1255418/58‑1 (MQ=255)
                                                  aGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGGTGTTTAc                >  1:2337066/1‑68 (MQ=255)
                                                  aGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGGTGTTTACGAt             >  1:1274576/1‑71 (MQ=255)
                                                  aGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGGTGTTTACGAt             >  1:2005861/1‑71 (MQ=255)
                                                  aGCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGGTGTTTACGAt             >  1:143953/1‑71 (MQ=255)
                                                  aGCAAAGAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGGTGTTTACGAt             >  1:24/1‑71 (MQ=255)
                                                   gCAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGGTGTTTACGAt             >  1:1555646/1‑70 (MQ=255)
                                                    cAAATAAAACAGTAGGACGATGACTACTGGAAGGAAGGaaa                                         <  1:1127320/41‑1 (MQ=255)
                                                       aTAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTgg                              <  1:428758/49‑1 (MQ=255)
                                                       aTAAAACAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGGTGTTTACGATCAGtt        <  1:1661627/71‑1 (MQ=255)
                                                         aaaaCAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGa                                   <  1:1864265/42‑1 (MQ=255)
                                                            aCAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGGTGTTTACGATCAGTTAAgcg   >  1:1719036/1‑71 (MQ=255)
                                                             cAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGGTGTTTACGATCAGTTAAgcgc  <  1:2277404/71‑1 (MQ=255)
                                                             cAGTAGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGGTGTTTACGATCAGTTAAgcgc  <  1:1864858/71‑1 (MQ=255)
                                                                 |                                                                  
TCCGGCCTGGTCGCCAGCCAGAGCACAGGAAGAAATTTCCCGCCTGCAACAGCAAATAAAACAGTGGGACGATGACTACTGGAAGGAAGGAAAAAGTGAGGTGGAAGACGGTGTTTACGATCAGTTAAGCGC  >  minE/2520369‑2520500

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: