Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,826 A→T 6.9% intergenic (‑363/+91) 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,8260AT 6.9% 163.9 / 4.2 58intergenic (‑363/+91)ldrA/ldrCtoxic polypeptide, small/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base T (2/2);  total (30/28)
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.38e-01

GCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGG  >  minE/821760‑821894
                                                                  |                                                                    
gCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGAt                                                                   <  1:1542625/70‑1 (MQ=255)
 cAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATa                                                                  <  1:3177636/70‑1 (MQ=255)
 cAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATa                                                                  <  1:122937/70‑1 (MQ=255)
    tCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATa                                                                  <  1:827126/67‑1 (MQ=255)
     cGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGa                                                              >  1:3876598/1‑70 (MQ=255)
      gTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATaaa                                                                <  1:3553379/67‑1 (MQ=255)
      gTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATaaa                                                                <  1:3558366/67‑1 (MQ=255)
      gTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGAc                                                             >  1:1534259/1‑70 (MQ=255)
      gTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACa                                                            <  1:3605144/71‑1 (MQ=255)
       tttCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACaa                                                           <  1:1550834/71‑1 (MQ=255)
        ttCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACaa                                                           >  1:805594/1‑70 (MQ=255)
        ttCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACaa                                                           >  1:2815770/1‑70 (MQ=255)
        ttCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACaa                                                           >  1:1166661/1‑70 (MQ=255)
         tCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATa                                                                  >  1:1257833/1‑62 (MQ=255)
         tCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACaaa                                                          >  1:3073055/1‑70 (MQ=255)
          cATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACaaaa                                                         <  1:2844917/70‑1 (MQ=255)
            tGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATg                                                      <  1:1786145/71‑1 (MQ=255)
              aGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGc                                                    >  1:1938505/1‑71 (MQ=255)
                     tGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATa                                                                  <  1:1876050/50‑1 (MQ=255)
                     tGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAg                                                               <  1:1144989/53‑1 (MQ=255)
                     tGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGTCAAAAATGGCCGCAGGc                                             >  1:1091721/1‑71 (MQ=255)
                     tGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACaaa                                                          <  1:3680061/58‑1 (MQ=255)
                     tGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGcc                                                   >  1:2308534/1‑65 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGAt                                                                   <  1:4142411/47‑1 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGAt                                                                   <  1:67543/47‑1 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACaaaaa                                                        >  1:3037820/1‑58 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACaaaaa                                                        >  1:2747211/1‑58 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCt                                            <  1:67060/70‑1 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCt                                            <  1:286879/70‑1 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTg                                           >  1:311712/1‑71 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTg                                           >  1:612083/1‑71 (MQ=255)
                       gatgatGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTg                                           >  1:306837/1‑71 (MQ=255)
                        atgatgCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGt                                          <  1:940345/71‑1 (MQ=255)
                           atgCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACaaa                                                          <  1:2707590/52‑1 (MQ=255)
                           atgCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTAc                                       <  1:2900802/71‑1 (MQ=255)
                              cTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAAATAAAGACAAAAATGGCCGCAGGCTGTTAcccc                                    <  1:1868178/71‑1 (MQ=255)
                                tCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATaaa                                                                >  1:2391684/1‑41 (MQ=255)
                                 ccAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGc                                 >  1:2042186/1‑71 (MQ=255)
                                    gCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACaaa                                                          <  1:2643939/43‑1 (MQ=255)
                                       cGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCgg                           >  1:2433141/1‑71 (MQ=255)
                                        gggTTGCCAATATTTACTTGTGGAAGAGATAAAGACaaaaa                                                        >  1:1257041/1‑41 (MQ=255)
                                        gggTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGc                                             <  1:2663213/52‑1 (MQ=255)
                                        gggTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCt                                   >  1:2502319/1‑62 (MQ=255)
                                          gTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGttt                        <  1:872903/71‑1 (MQ=255)
                                          gTTGCCAATATTTACTTGTGGAAGAGATAAAGACaaa                                                          <  1:1555046/37‑1 (MQ=255)
                                             gCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCg                      >  1:3548500/1‑70 (MQ=255)
                                               cAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGgc                   >  1:1494606/1‑71 (MQ=255)
                                                 atatTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGgcg                  >  1:195728/1‑70 (MQ=255)
                                                 atatTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGgcg                  >  1:688340/1‑70 (MQ=255)
                                                 atatTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGgcg                  >  1:1805901/1‑70 (MQ=255)
                                                  tatTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCa                <  1:2654123/71‑1 (MQ=255)
                                                   atTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCa                                                >  1:1974344/1‑38 (MQ=255)
                                                   atTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCa                                                >  1:490210/1‑38 (MQ=255)
                                                   atTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCa                >  1:3641521/1‑70 (MQ=255)
                                                    tttACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTa                                        >  1:2595143/1‑45 (MQ=255)
                                                     ttACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCggg                    >  1:990995/1‑64 (MQ=255)
                                                      tACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCatat             <  1:2571420/70‑1 (MQ=255)
                                                           gtgGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCa        <  1:487707/70‑1 (MQ=255)
                                                                  agaTAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACgg  >  1:1245081/1‑69 (MQ=255)
                                                                  |                                                                    
GCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGG  >  minE/821760‑821894

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: