Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,728 A→G 24.6% intergenic (‑265/+189) 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,7280AG24.6% 139.5 / 40.4 73intergenic (‑265/+189)ldrA/ldrCtoxic polypeptide, small/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base A (25/30);  new base G (9/9);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCT  >  minE/821658‑821791
                                                                      |                                                               
ggATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATa                                                                 <  1:304988/71‑1 (MQ=255)
 gATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATa                                                                 >  1:2825760/1‑70 (MQ=255)
     aTTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAag                                                              <  1:2725243/69‑1 (MQ=255)
               ggTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCgg                                                  >  1:2642846/1‑71 (MQ=39)
                     gATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCtttt                                             >  1:1104757/1‑70 (MQ=40)
                     gATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTc                                            <  1:5623243/71‑1 (MQ=39)
                     gATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTc                                            <  1:5258918/71‑1 (MQ=39)
                       tactacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCaa                                          >  1:3703521/1‑71 (MQ=39)
                       tactacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCaa                                          >  1:3998554/1‑71 (MQ=39)
                        actacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCaa                                          >  1:4528678/1‑70 (MQ=40)
                         ctacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAAc                                         <  1:1392233/70‑1 (MQ=40)
                            cGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTa                                                          >  1:359793/1‑50 (MQ=255)
                            cGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTa                                                          >  1:4662579/1‑50 (MQ=255)
                                 cAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATaaa                                                               >  1:656551/1‑40 (MQ=40)
                                 cAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATaaa                                                               >  1:4800916/1‑40 (MQ=40)
                                  aaTTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTg                                 <  1:2788131/69‑1 (MQ=35)
                                     tCATCTTCACTTCATCCCTGAAATGTTTGCAATaa                                                                >  1:3548618/1‑35 (MQ=37)
                                     tCATCTTCACTTCATCCCTGAAATGTTTGCAATaa                                                                >  1:377284/1‑35 (MQ=37)
                                     tCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTa                                                          <  1:703983/41‑1 (MQ=38)
                                     tCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTa                                                          <  1:3211241/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:523176/70‑1 (MQ=255)
                                             acttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCAt                     <  1:4844769/70‑1 (MQ=255)
                                             acttcaTCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCAt                     <  1:2873495/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:83291/1‑45 (MQ=38)
                                                     ccTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGc                                      >  1:5061367/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:2800009/1‑67 (MQ=255)
                                                        gAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGctc             >  1:3716784/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:3052047/71‑1 (MQ=255)
                                                         aaaTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgat        <  1:5206207/71‑1 (MQ=255)
                                                         aaaTGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgat        <  1:3996150/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)
                                                                      aaaGAGTACATTCCGGCTTTTCAACAGCTGTTGCAg                              >  1:2203093/1‑36 (MQ=33)
                                                                      |                                                               
GGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCT  >  minE/821658‑821791

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: