Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,605 T→C 15.2% intergenic (‑142/+312) ldrA ← / ← ldrC toxic polypeptide, small/toxic polypeptide, small

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE821,6050TC15.2% 90.1 / 11.0 46intergenic (‑142/+312)ldrA/ldrCtoxic polypeptide, small/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base T (22/17);  new base C (3/4);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

GTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGT  >  minE/821537‑821675
                                                                    |                                                                      
ttCTGGTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTc                                                                      <  1:2163854/69‑1 (MQ=11)
     gTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGc                                                                  >  1:1714470/1‑70 (MQ=35)
      gttCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGc                                                                  <  1:4008734/68‑1 (MQ=255)
       ttCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAc                                                               <  1:4248633/71‑1 (MQ=21)
         cAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                                                              >  1:2561353/1‑70 (MQ=35)
         cAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                                                              <  1:848013/70‑1 (MQ=35)
         cAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                                                              <  1:1348151/70‑1 (MQ=35)
          aaGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCa                                                             <  1:3345401/70‑1 (MQ=25)
            gATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATg                                                          >  1:2076834/1‑71 (MQ=34)
            gATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATg                                                          <  1:2913665/71‑1 (MQ=21)
             aTTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGc                                                         >  1:2321658/1‑71 (MQ=21)
             aTTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGc                                                         <  1:2163303/71‑1 (MQ=21)
               tAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGcca                                                       >  1:1831334/1‑71 (MQ=21)
               tAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGcca                                                       >  1:2765249/1‑71 (MQ=21)
                 gCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGccac                                                      <  1:1062496/70‑1 (MQ=35)
                       gTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGTGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGa                                                >  1:4432181/1‑70 (MQ=255)
                        ttCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGAt                                               <  1:4036196/70‑1 (MQ=25)
                        ttCTGTTGTAAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCtt                                                                        >  1:863469/1‑45 (MQ=255)
                         tCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAc                                                               >  1:5281153/1‑53 (MQ=16)
                          cTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATa                                              >  1:559051/1‑69 (MQ=14)
                          cTGTTGTCAGGTTGTACCTATCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGc                                                         >  1:205662/1‑58 (MQ=11)
                                   ggTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCa                                   <  1:4729267/71‑1 (MQ=25)
                                   aggtttACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGc                                    <  1:1827786/67‑1 (MQ=25)
                                     ttgtACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAAc                                 <  1:2443606/68‑1 (MQ=255)
                                       gtACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACa                                <  1:5293143/69‑1 (MQ=255)
                                           ctctCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTgc                            >  1:4618189/1‑70 (MQ=40)
                                             ctcAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGcc                         >  1:5364258/1‑71 (MQ=39)
                                               cAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCt                        >  1:829715/1‑70 (MQ=40)
                                               cAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCt                        >  1:3491571/1‑70 (MQ=40)
                                               cAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCt                        >  1:4857575/1‑70 (MQ=40)
                                                aaCGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTc                       <  1:2141498/70‑1 (MQ=40)
                                                aaCGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCa                      >  1:613276/1‑71 (MQ=39)
                                                aaCGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCa                      <  1:1142705/71‑1 (MQ=39)
                                                aaCGTGCGGGGGTTTTCTCTCTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCa                      >  1:4978412/1‑71 (MQ=255)
                                                 aCGTGCGGGGGTTTTCTCTCTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCAc                     >  1:216469/1‑71 (MQ=255)
                                                   gTGCGGGGGTTTTCTCTCTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCg                   >  1:4465944/1‑71 (MQ=255)
                                                         gggTTTTCTCTCTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGa                 <  1:1113653/67‑1 (MQ=255)
                                                            ttttCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTgcgc                          <  1:708275/55‑1 (MQ=255)
                                                            ttttCTCTCTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCaa                                  <  1:2124949/47‑1 (MQ=255)
                                                            ttttCTCTCTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCaa                                  <  1:3128091/47‑1 (MQ=255)
                                                            ttttCTCTCTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCaa                                  <  1:3003723/47‑1 (MQ=255)
                                                              ttCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCgg        >  1:3857988/1‑71 (MQ=255)
                                                               tctcTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGa                 >  1:4914605/1‑61 (MQ=40)
                                                                ctctctCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCaa                                  <  1:124957/43‑1 (MQ=255)
                                                                ctctctCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCaa                                  <  1:73328/43‑1 (MQ=255)
                                                                ctctTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACAt                               >  1:2605328/1‑46 (MQ=255)
                                                                ctctTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCAcc                    >  1:4710655/1‑57 (MQ=255)
                                                                ctctTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCAcc                    >  1:2319123/1‑57 (MQ=255)
                                                                 tctTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTgcg                           >  1:3196446/1‑49 (MQ=255)
                                                                   tttCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCAc                     >  1:5163547/1‑53 (MQ=255)
                                                                   tttCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCt      <  1:4219710/68‑1 (MQ=255)
                                                                    ttCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt  >  1:1874567/1‑71 (MQ=255)
                                                                    |                                                                      
GTCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGT  >  minE/821537‑821675

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: