Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,688 T→A 21.0% intergenic (‑225/+229) 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,6880TA21.0% 127.6 / 31.1 62intergenic (‑225/+229)ldrA/ldrCtoxic polypeptide, small/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base T (19/30);  new base A (4/9);  total (23/39)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAAC  >  minE/821621‑821752
                                                                   |                                                                
caccaGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGtct                                                                <  1:3694624/70‑1 (MQ=255)
     gggATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTc                                                          <  1:3247632/71‑1 (MQ=255)
     gggATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTc                                                          <  1:1960772/71‑1 (MQ=255)
       gATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCAt                                                        >  1:3753573/1‑71 (MQ=255)
        aTAAAGCCCCCGCAACATTGCGCCTCACCGGATAACGCCGGCTTGGTGTGGATACTACGACTCAATTCATc                                                       >  1:420335/1‑71 (MQ=11)
             gCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcac                                                  >  1:5298653/1‑71 (MQ=255)
                 ccGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcactt                                                <  1:5222463/69‑1 (MQ=255)
                 ccGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATctt                                                     <  1:5470837/64‑1 (MQ=255)
                 ccGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATctt                                                     <  1:357918/64‑1 (MQ=255)
                  cGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCAt                                             >  1:851332/1‑71 (MQ=255)
                    cAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATcc                                           >  1:1706426/1‑71 (MQ=255)
                    cAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATcc                                           <  1:4355330/71‑1 (MQ=255)
                     aaCATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTAACTTCATccc                                          <  1:2057318/71‑1 (MQ=255)
                      aCATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGTATACTACGt                                                                  >  1:4644459/1‑46 (MQ=255)
                       cATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCt                                         >  1:4058878/1‑70 (MQ=255)
                       cATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCt                                         >  1:5251838/1‑70 (MQ=255)
                       cATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCt                                         >  1:3173459/1‑70 (MQ=255)
                        aTTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCaa                                                             <  1:3422698/49‑1 (MQ=255)
                             gcCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGt                                  <  1:1405894/71‑1 (MQ=255)
                             gcCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGt                                  <  1:5009579/71‑1 (MQ=255)
                                  aCCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                              <  1:5444792/70‑1 (MQ=255)
                                  aCCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCa                             <  1:858252/71‑1 (MQ=255)
                                    cGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCaa                            <  1:3770626/70‑1 (MQ=255)
                                    cGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCaa                            <  1:2101632/70‑1 (MQ=255)
                                     ggATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATa                          <  1:304988/71‑1 (MQ=255)
                                      gATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGaa                                      <  1:1160941/58‑1 (MQ=255)
                                      gATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGaa                                      <  1:1308997/58‑1 (MQ=255)
                                      gATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGa                                       <  1:5187130/57‑1 (MQ=255)
                                      gATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGa                                       <  1:233630/57‑1 (MQ=255)
                                      gATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATa                          >  1:2825760/1‑70 (MQ=255)
                                        tAATTCCGGCTTGGTGTGGATACTACGTCTCAAtt                                                           >  1:1965953/1‑35 (MQ=255)
                                        tAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAAt                                    >  1:2051320/1‑58 (MQ=255)
                                          aTTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcacttc                                               <  1:2722744/45‑1 (MQ=255)
                                          aTTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcacttc                                               <  1:1247612/45‑1 (MQ=255)
                                          aTTCCGGCTTGGTGTGGATACTACGTCTCAATTCATcttcacttc                                               <  1:147074/45‑1 (MQ=255)
                                          aTTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAag                       <  1:2725243/69‑1 (MQ=255)
                                           ttCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGttt                                <  1:1147776/59‑1 (MQ=40)
                                           ttCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGttt                                <  1:3641226/59‑1 (MQ=40)
                                           ttCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCaa                            <  1:4671516/63‑1 (MQ=39)
                                              cGGCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGttt                                <  1:5044334/56‑1 (MQ=255)
                                                gCTTGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAAtg                          <  1:1363072/60‑2 (MQ=255)
                                                   tGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGaaa                                     >  1:2147106/1‑46 (MQ=255)
                                                   tGGTGTGGATACTACGACTCAATTCATCTTCACTTCATCCCTGaaa                                     >  1:2633534/1‑46 (MQ=255)
                                                    ggTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCgg           >  1:2642846/1‑71 (MQ=39)
                                                      tgtgGATACTACGACTCAATTCATCTTCACTTCATCCCTGa                                       <  1:807495/41‑1 (MQ=255)
                                                      tgtgGATACTACGACTCAATTCATCTTCACTTCATCCCTGa                                       <  1:4514396/41‑1 (MQ=255)
                                                      tgtgGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                              <  1:2035275/50‑1 (MQ=255)
                                                      tgtgGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                              <  1:3590591/50‑1 (MQ=255)
                                                      tgtgGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                              <  1:5193936/50‑1 (MQ=255)
                                                        tgGATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATg                                   <  1:2183216/43‑1 (MQ=255)
                                                          gATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCtttt      >  1:1104757/1‑70 (MQ=40)
                                                          gATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTc     <  1:5258918/71‑1 (MQ=39)
                                                          gATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTc     <  1:5623243/71‑1 (MQ=39)
                                                          gATACTACGACTCAATTCATCTTCACTTCATCCCTGAAATGt                                  >  1:4679074/1‑42 (MQ=255)
                                                            tactacGTCTCAATTCATCTTCACTTCATCCCTGAAATGttt                                >  1:3600496/1‑42 (MQ=38)
                                                            tactacGTCTCAATTCATCTTCACTTCATCCCTGAAATGttt                                >  1:831210/1‑42 (MQ=38)
                                                            tactacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCaa   >  1:3998554/1‑71 (MQ=39)
                                                            tactacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCaa   >  1:3703521/1‑71 (MQ=39)
                                                             actacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCaa   >  1:4528678/1‑70 (MQ=40)
                                                              ctacGTCTCAATTCATCTTCACTTCATCCCTGAAATGtt                                 <  1:731949/39‑1 (MQ=40)
                                                              ctacGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAAc  <  1:1392233/70‑1 (MQ=40)
                                                              ctacGACTCAATTCATCTTCACTTCATCCCTGAAATGTTTGc                              <  1:1785173/42‑1 (MQ=255)
                                                                 cGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTa                   >  1:4662579/1‑50 (MQ=255)
                                                                 cGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTa                   >  1:359793/1‑50 (MQ=255)
                                                                   |                                                                
CACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAAC  >  minE/821621‑821752

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: