Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 949,284 T→A 5.1% Y518F (TAT→TTT)  dos ← cAMP phosphodiesterase, heme‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE949,2840TA 5.1% 531.9 / 14.6 177Y518F (TAT→TTT) doscAMP phosphodiesterase, heme‑regulated
Reads supporting (aligned to +/- strand):  ref base T (83/85);  new base A (7/2);  total (90/87)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

Only 100 of 183 total aligned reads displayed.

TTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCC  >  minE/949220‑949351
                                                                |                                                                   
ttttACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATc                                                                 <  1:4853207/69‑1 (MQ=255)
ttttACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCca                                                               >  1:2404257/1‑71 (MQ=255)
  ttACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAAt                                                                  >  1:4494184/1‑66 (MQ=255)
         ttCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATtgt                                                      >  1:1304276/1‑71 (MQ=255)
         ttCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATtgt                                                      >  1:1889258/1‑71 (MQ=255)
          tCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCcattgcatt                                                        >  1:3801805/1‑68 (MQ=255)
          tCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCcattgcatt                                                        >  1:3259066/1‑68 (MQ=255)
          tCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCcattgcatt                                                        >  1:3921946/1‑68 (MQ=255)
          tCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATtgt                                                      <  1:4013923/70‑1 (MQ=255)
          tCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATtgt                                                      <  1:2610093/70‑1 (MQ=255)
           cGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATtgt                                                      <  1:4059801/69‑1 (MQ=255)
                 tCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              <  1:2204591/71‑1 (MQ=255)
                  cGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGt                                               >  1:4447356/1‑69 (MQ=255)
                  cGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              >  1:3428763/1‑70 (MQ=255)
                  cGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              >  1:3647636/1‑70 (MQ=255)
                   gCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCca                                                               >  1:907889/1‑52 (MQ=255)
                   gCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              >  1:3172090/1‑69 (MQ=255)
                   gCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              >  1:195857/1‑69 (MQ=255)
                   gCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGga                                            <  1:832786/71‑1 (MQ=255)
                     cGGGCTGAAGAACTGCCAGCCGTTACCGCCATTTTTGCGAATAAAATcc                                                                >  1:238329/1‑49 (MQ=255)
                      gggCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGagag                                         >  1:1837520/1‑71 (MQ=255)
                      gggCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGagag                                         >  1:2154142/1‑71 (MQ=255)
                      gggCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGagag                                         >  1:3543590/1‑71 (MQ=255)
                       ggCTGAAGAACTGCCAGCCGTTACCGCCATTTTTGCGAATAAAATCCATTGCATTGTGAGCAGTGGAGAGc                                        <  1:3418388/71‑1 (MQ=255)
                       ggCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              <  1:3951963/65‑1 (MQ=255)
                       ggCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              <  1:3426810/65‑1 (MQ=255)
                       ggCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGc                                        <  1:1665404/71‑1 (MQ=255)
                         cTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCcattt                                                            <  1:4188905/49‑2 (MQ=255)
                         cTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCcattg                                                            <  1:1379516/49‑1 (MQ=255)
                         cTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCcattg                                                            <  1:3221560/49‑1 (MQ=255)
                         cTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCcattg                                                            <  1:2110009/49‑1 (MQ=255)
                            aagaaCTGCCAGCCGTTACCGCCATTCTTGCGAAtata                                                                    >  1:1688614/1‑38 (MQ=255)
                            aagaaCTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCca                                                               >  1:1075228/1‑43 (MQ=255)
                            aagaaCTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCca                                                               >  1:1269457/1‑43 (MQ=255)
                            aagaaCTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTa                                   <  1:1390422/71‑1 (MQ=255)
                            aagaaCTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTa                                   <  1:3305071/71‑1 (MQ=255)
                                aCTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCcat                                                              >  1:3784309/1‑40 (MQ=255)
                                aCTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAg                                                   >  1:1476048/1‑51 (MQ=255)
                                aCTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAg                                                   >  1:94200/1‑51 (MQ=255)
                                aCTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcg                               <  1:3082458/71‑1 (MQ=255)
                                aCTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcg                               <  1:491823/71‑1 (MQ=255)
                                aCTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcg                               <  1:4043950/71‑1 (MQ=255)
                                aCTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcg                               <  1:2519010/71‑1 (MQ=255)
                                aCTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcg                               <  1:1811619/71‑1 (MQ=255)
                                aCTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcg                               <  1:806614/71‑1 (MQ=255)
                                 cTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgc                              <  1:99338/71‑1 (MQ=255)
                                 cTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgc                              <  1:2016355/71‑1 (MQ=255)
                                 cTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgc                              <  1:2022420/71‑1 (MQ=255)
                                 cTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgc                              <  1:4969797/71‑1 (MQ=255)
                                  tGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGgaga                                          <  1:632008/58‑1 (MQ=255)
                                  tGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGgaga                                          <  1:3123118/58‑1 (MQ=255)
                                  tGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGgaga                                          <  1:3111264/58‑1 (MQ=255)
                                    ccAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGt                           <  1:1161079/71‑1 (MQ=255)
                                    ccAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGt                           <  1:4568043/71‑1 (MQ=255)
                                       gCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              <  1:4024905/49‑1 (MQ=255)
                                       gCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              <  1:4689188/49‑1 (MQ=255)
                                       gCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              <  1:4519679/49‑1 (MQ=255)
                                       gCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              <  1:3886360/49‑1 (MQ=255)
                                       gCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              <  1:4675600/49‑1 (MQ=255)
                                       gCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              <  1:1115038/49‑1 (MQ=255)
                                       gCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              <  1:2003048/49‑1 (MQ=255)
                                       gCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              <  1:334784/49‑1 (MQ=255)
                                       gCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTg                                              <  1:3164259/49‑1 (MQ=255)
                                       gCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGtttt                        <  1:642695/71‑1 (MQ=255)
                                       gCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGtttt                        <  1:3473689/71‑1 (MQ=255)
                                       gCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGtttt                        <  1:2887414/71‑1 (MQ=255)
                                        ccGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGtttt                        <  1:2830747/70‑1 (MQ=255)
                                         cGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTa                      >  1:3778022/1‑71 (MQ=255)
                                                   aTTCTTGCGAATATAATCCATTGCATTGTGAGCAg                                                >  1:1807322/1‑35 (MQ=255)
                                                   aTTCTTGCGAATATAATCCATTGCATTGTGAGCAg                                                >  1:726236/1‑35 (MQ=255)
                                                   aTTCTTGCGAATATAATCCATTGCATTGTGAGCAg                                                >  1:5150759/1‑35 (MQ=255)
                                                   aTTCTTGCGAATATAATCCATTGCATTGTGAGCAg                                                >  1:2356825/1‑35 (MQ=255)
                                                   aTTCTTGCGAATATAATCCATTGCATTGTGAGCAg                                                >  1:482657/1‑35 (MQ=255)
                                                   aTTCTTGCGAATATAATCCATTGCATTGTGAGCAg                                                >  1:2773818/1‑35 (MQ=255)
                                                   aTTCTTGCGAATATAATCCATTGCATTGTGAGCAg                                                >  1:2791261/1‑35 (MQ=255)
                                                   aTTCTTGCGAATATAATCCATTGCATTGTGAGCAg                                                >  1:945000/1‑35 (MQ=255)
                                                   aTTCTTGCGAATATAATCCATTGCATTGTGAGCAg                                                >  1:3513873/1‑35 (MQ=255)
                                                   aTTCTTGCGAATATAATCCATTGCATTGTGAGCAg                                                >  1:35478/1‑35 (MQ=255)
                                                   aTTCTTGCGAATATAATCCATTGCATTGTGAGCAg                                                >  1:4346069/1‑35 (MQ=255)
                                                   aTTCTTGCGAATATAATCCATTGCATTGTGAGCAg                                                >  1:3615542/1‑35 (MQ=255)
                                                   aTTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGACGt            <  1:5090155/71‑1 (MQ=255)
                                                   aTTCTTACGAATATAATCCATTGCATTGTGAGCAg                                                >  1:669841/1‑35 (MQ=255)
                                                    ttCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcg                               <  1:4027638/51‑1 (MQ=255)
                                                    ttCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTa           <  1:2718/71‑1 (MQ=255)
                                                      cTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgcg                             >  1:2051324/1‑51 (MQ=255)
                                                      cTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgcg                             >  1:4104916/1‑51 (MQ=255)
                                                      cTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgcg                             >  1:4403292/1‑51 (MQ=255)
                                                      cTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGtttt                        >  1:4289661/1‑56 (MQ=255)
                                                      cTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGc         >  1:4862765/1‑71 (MQ=255)
                                                      cTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGc         >  1:2627983/1‑71 (MQ=255)
                                                      cTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGc         >  1:2934958/1‑71 (MQ=255)
                                                        tGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCt        >  1:5059332/1‑70 (MQ=255)
                                                         gCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCt        <  1:1528482/69‑1 (MQ=255)
                                                         gCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCt        <  1:545286/69‑1 (MQ=255)
                                                         gCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCt        <  1:4697918/69‑1 (MQ=255)
                                                         gCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCt        <  1:4361120/69‑1 (MQ=255)
                                                             atatAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGcc  <  1:1420626/71‑1 (MQ=255)
                                                              tataATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGttttt                       <  1:310320/49‑1 (MQ=255)
                                                              tataATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGttttt                       <  1:3686532/49‑1 (MQ=255)
                                                              tataATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGt                           <  1:2273558/45‑1 (MQ=255)
                                                                |                                                                   
TTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCC  >  minE/949220‑949351

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: