Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,735 A→G 32.5% intergenic (‑272/+182) 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,7350AG32.5% 113.8 / 63.8 80intergenic (‑272/+182)ldrA/ldrCtoxic polypeptide, small/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base G (12/14);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTG  >  minE/821673‑821805
                                                              |                                                                      
ggTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCgg                                                                >  1:2642846/1‑71 (MQ=39)
      gATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCtttt                                                           >  1:1104757/1‑70 (MQ=40)
      gATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTc                                                          <  1:5258918/71‑1 (MQ=39)
      gATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTc                                                          <  1:5623243/71‑1 (MQ=39)
        tactacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCaa                                                        >  1:3998554/1‑71 (MQ=39)
        tactacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCaa                                                        >  1:3703521/1‑71 (MQ=39)
         actacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCaa                                                        >  1:4528678/1‑70 (MQ=40)
          ctacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAAc                                                       <  1:1392233/70‑1 (MQ=40)
             cGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTa                                                                        >  1:4662579/1‑50 (MQ=255)
             cGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTa                                                                        >  1:359793/1‑50 (MQ=255)
                   aaTTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTg                                               <  1:2788131/69‑1 (MQ=35)
                      tCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTa                                                                        <  1:3211241/41‑1 (MQ=38)
                      tCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTa                                                                        <  1:703983/41‑1 (MQ=38)
                       cATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCtttt                                                           >  1:4463368/1‑53 (MQ=35)
                       cATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCtttt                                                           >  1:2354029/1‑53 (MQ=35)
                       cATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACa                                                      >  1:3531416/1‑58 (MQ=35)
                       cATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGc                                                    >  1:636252/1‑60 (MQ=35)
                        aTCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTa                                                                        <  1:1923710/39‑1 (MQ=40)
                        aTCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGt                                        <  1:3024025/71‑1 (MQ=39)
                           ttcacttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGtt                                       >  1:3937129/1‑69 (MQ=40)
                           ttcacttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGt                                        <  1:4887246/68‑1 (MQ=40)
                           ttcacttcaTCCCTGAAATGTTTGCAATAAAGAGAACATTCCGGCTTTTCAACAGCTGTTGCAGTCGtt                                       >  1:320400/1‑69 (MQ=255)
                            tcacttcaTCCCTGAAATGTTTGCAATAAAGAGAACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTc                                     <  1:576134/70‑1 (MQ=255)
                             cacttcaTCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCa                                    >  1:12353/1‑70 (MQ=11)
                             cacttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTcc                                                                  >  1:797557/1‑40 (MQ=40)
                             cacttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGttt                                      <  1:1446410/68‑1 (MQ=40)
                             cacttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTc                                     <  1:766275/69‑1 (MQ=40)
                             cacttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTc                                     <  1:4797940/69‑1 (MQ=40)
                             cacttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTc                                     <  1:2007512/69‑1 (MQ=40)
                             cacttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCa                                    >  1:2046808/1‑70 (MQ=255)
                             cacttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCAt                                   >  1:2184858/1‑71 (MQ=255)
                              acttcaTCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCAt                                   <  1:5603295/70‑1 (MQ=11)
                              acttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCAt                                   <  1:4844769/70‑1 (MQ=255)
                              acttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCAt                                   <  1:2873495/70‑1 (MQ=255)
                              acttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCAt                                   <  1:523176/70‑1 (MQ=255)
                              acttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATg                                  <  1:2383303/71‑1 (MQ=255)
                              acttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATg                                  <  1:216736/71‑1 (MQ=255)
                               cttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTc                                     <  1:5234717/67‑1 (MQ=40)
                               cttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGa                                 <  1:1444423/71‑1 (MQ=255)
                                  caTCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGa                                 >  1:1055425/1‑68 (MQ=11)
                                  caTCCCTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTg                              >  1:3842299/1‑71 (MQ=14)
                                   aTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGc                             <  1:1585178/71‑1 (MQ=255)
                                    tCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCAt                                   <  1:4156454/64‑1 (MQ=255)
                                     cccTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGct                            >  1:5006544/1‑70 (MQ=255)
                                      ccTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGc                                                    >  1:5061367/1‑45 (MQ=38)
                                      ccTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGc                                                    >  1:83291/1‑45 (MQ=38)
                                       cTGAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTg                         <  1:4718052/71‑1 (MQ=17)
                                        tGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTc                                                          <  1:1754790/37‑1 (MQ=40)
                                         gAAATGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGga                       >  1:2975512/1‑71 (MQ=17)
                                         gAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGctc                           >  1:3716784/1‑67 (MQ=255)
                                         gAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGctc                           >  1:2800009/1‑67 (MQ=255)
                                          aaaTGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCAt                                   <  1:1047272/58‑1 (MQ=11)
                                          aaaTGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGgat                      <  1:976120/71‑1 (MQ=17)
                                          aaaTGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGgat                      <  1:2869499/71‑1 (MQ=17)
                                          aaaTGTTTGCAATAAAGAGTACATTCCGGCTTTTCaa                                                        >  1:5271687/1‑37 (MQ=40)
                                          aaaTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTgg                        <  1:4647999/69‑1 (MQ=255)
                                          aaaTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgat                      <  1:3996150/71‑1 (MQ=255)
                                          aaaTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgat                      <  1:3052047/71‑1 (MQ=255)
                                          aaaTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgat                      <  1:5206207/71‑1 (MQ=255)
                                          aaaTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGga                       >  1:807401/1‑70 (MQ=255)
                                           aaTGTTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGgatg                     <  1:3749695/71‑1 (MQ=18)
                                               tttGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGc                 >  1:385316/1‑71 (MQ=18)
                                                ttGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGgatg                     >  1:4706826/1‑66 (MQ=18)
                                                ttGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCt                >  1:1955198/1‑71 (MQ=18)
                                                ttGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGc                                              >  1:2712986/1‑41 (MQ=38)
                                                 tGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGgatgag                   >  1:4322892/1‑66 (MQ=18)
                                                 tGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGgatgag                   >  1:5250076/1‑66 (MQ=18)
                                                 tGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCt                <  1:2690654/70‑1 (MQ=18)
                                                 tGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTg                              >  1:4844364/1‑56 (MQ=255)
                                                    aaTGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTg                         <  1:1423380/58‑1 (MQ=17)
                                                    aaTGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTg                         <  1:4827439/58‑1 (MQ=17)
                                                     aTAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGc                                              <  1:4906049/36‑1 (MQ=40)
                                                      tgaaGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATg                                  <  1:4325493/45‑1 (MQ=14)
                                                      tgaaGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATg                                  <  1:1688084/45‑1 (MQ=14)
                                                       aaaGAGTACATTCCGGCTTTTCAACAGCTGTTGCAg                                            >  1:2203093/1‑36 (MQ=33)
                                                        aaGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGgatga                    >  1:504803/1‑59 (MQ=21)
                                                        aaGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTc        <  1:5614866/71‑1 (MQ=25)
                                                          gagTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCgg      <  1:3616340/71‑1 (MQ=37)
                                                          gagTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCgg      >  1:3947259/1‑71 (MQ=37)
                                                          gagTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCgg      <  1:329877/71‑1 (MQ=37)
                                                          gagTGCATTCCGGTTTTTAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCggg     >  1:1143634/1‑71 (MQ=255)
                                                           agTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCa                                    <  1:4530933/40‑1 (MQ=255)
                                                           agTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgatga                    <  1:2635013/56‑1 (MQ=255)
                                                             tACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGt                               >  1:1192865/1‑43 (MQ=255)
                                                             tACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGatgatg                  >  1:4154224/1‑56 (MQ=255)
                                                              aCATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTg  >  1:4311177/1‑71 (MQ=255)
                                                              |                                                                      
GGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTG  >  minE/821673‑821805

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: