Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,581,504 C→T 9.1% L577L (TTG→TTA ligA ← DNA ligase, NAD(+)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,581,5040CT 9.1% 159.6 / 6.6 55L577L (TTG→TTAligADNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base C (21/29);  new base T (2/3);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGAC  >  minE/1581435‑1581572
                                                                     |                                                                    
aCCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGcc                                                                      <  1:1510174/70‑1 (MQ=255)
      aaaCGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGcc                                                                      <  1:1403922/64‑1 (MQ=255)
         cGGGCTGTCAATCTCTTCCGCGTTGATAACGATCTGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGc                                                            >  1:1491087/1‑71 (MQ=255)
         cGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGc                                                            >  1:1285524/1‑71 (MQ=255)
         cGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGc                                                            <  1:1313848/71‑1 (MQ=255)
          gggCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCt                                                           <  1:2069877/71‑1 (MQ=255)
               gTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCtgatga                                                      >  1:1387533/1‑71 (MQ=255)
                tCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCtgatga                                                      <  1:73935/70‑1 (MQ=255)
                tCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGAc                                                     <  1:1552527/71‑1 (MQ=255)
                tCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGAc                                                     <  1:156135/71‑1 (MQ=255)
                tCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGAc                                                     <  1:1635350/71‑1 (MQ=255)
                tCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGAc                                                     <  1:647479/71‑1 (MQ=255)
                tCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGAc                                                     <  1:1637980/71‑1 (MQ=255)
                 cAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGAc                                                     <  1:181281/70‑1 (MQ=255)
                    tctcTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACATTCCGCCAACAGCTCGCTGATGACATTg                                                 >  1:736801/1‑71 (MQ=255)
                     ctctTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCTAACAGCTCGCTGATGACAtt                                                  >  1:675570/1‑69 (MQ=255)
                     ctctTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCtg                                                          <  1:356189/61‑1 (MQ=255)
                      tctTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCTAACAGCTCGCTGATGACATTGCg                                               >  1:1099102/1‑71 (MQ=255)
                          ccGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCtga                                                         >  1:1724355/1‑57 (MQ=255)
                          ccGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCtga                                                         >  1:1264316/1‑57 (MQ=255)
                           cgcgTTATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGAGGTTGCt                                         <  1:101318/71‑1 (MQ=255)
                            gcgTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCt                                         <  1:2204314/71‑1 (MQ=255)
                             cgTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCt                                         <  1:150329/70‑1 (MQ=255)
                                 gATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCTAACAGCTCGCTGATTACATTGCGGTTGCTttc                                      <  1:1438878/69‑1 (MQ=255)
                                 gATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCTAACAGCTCGCTGATGACATTGCGGTTGCTttc                                      <  1:1506802/69‑1 (MQ=255)
                                 gATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTttctt                                    >  1:37700/1‑71 (MQ=255)
                                 gATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTttctt                                    >  1:1102731/1‑71 (MQ=255)
                                 gATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTttctt                                    >  1:1454267/1‑71 (MQ=255)
                                  aTAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCt                                         <  1:1217930/65‑1 (MQ=255)
                                  aTAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCt                                         <  1:903582/65‑1 (MQ=255)
                                   tAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTttc                                      <  1:489572/67‑1 (MQ=255)
                                   tAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTttc                                      <  1:2178314/67‑1 (MQ=255)
                                     aCGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCgg                                 >  1:272000/1‑70 (MQ=255)
                                           ggCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAaaga                           <  1:793588/70‑1 (MQ=255)
                                           ggCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAaaga                           <  1:221210/70‑1 (MQ=255)
                                            gcgcAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAaagaag                         <  1:1295381/71‑1 (MQ=255)
                                            gcgcAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAaagaag                         <  1:1900630/71‑1 (MQ=255)
                                            gcgcAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAaagaag                         <  1:2180791/71‑1 (MQ=255)
                                            gcgcAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAaagaag                         <  1:615550/71‑1 (MQ=255)
                                                aGGCGAGTGAACACCTTCCGCCTACAGCTCGCtgatga                                                      <  1:1551784/38‑1 (MQ=38)
                                                    cAGTGAACACCTTCCGCTAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAAc                 <  1:1634699/71‑1 (MQ=255)
                                                    cAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGaa                  >  1:324424/1‑70 (MQ=255)
                                                    cAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAAc                 >  1:1066710/1‑71 (MQ=255)
                                                    cAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAAc                 >  1:109269/1‑71 (MQ=255)
                                                       tGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGaa                  >  1:1508361/1‑67 (MQ=255)
                                                       tGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGaa                  >  1:957851/1‑67 (MQ=255)
                                                         aaCACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTgg             <  1:1235782/70‑1 (MQ=255)
                                                           cacCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTttc                                      >  1:1845349/1‑43 (MQ=255)
                                                           cacCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTttc                                      >  1:1445928/1‑43 (MQ=255)
                                                           cacCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTttc                                      >  1:1311038/1‑43 (MQ=255)
                                                            acCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATg         >  1:515944/1‑71 (MQ=255)
                                                                 ccGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAaagaatt                        >  1:452135/1‑49 (MQ=255)
                                                                 ccGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAAc     <  1:461868/70‑1 (MQ=255)
                                                                 ccGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAAc     >  1:1495995/1‑70 (MQ=255)
                                                                  cGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACg    <  1:943648/70‑1 (MQ=255)
                                                                   gCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGAc  <  1:1330319/71‑1 (MQ=255)
                                                                   gCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGAc  <  1:946586/71‑1 (MQ=255)
                                                                     |                                                                    
ACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGAC  >  minE/1581435‑1581572

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: