Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,294,653 T→C 31.3% I198V (ATT→GTT)  nac ← DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,294,6530TC31.3% 70.3 / 40.7 83I198V (ATT→GTT) nacDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (36/21);  new base C (14/12);  total (50/33)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

CAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATCAACGCTTT  >  minE/1294587‑1294715
                                                                  |                                                              
cAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAg                                                              <  1:732322/69‑1 (MQ=255)
cAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAg                                                              <  1:3423058/69‑1 (MQ=255)
cAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAg                                                              <  1:3430730/69‑1 (MQ=255)
cAATTTCACCAATAACTTTTGCCGTGAGACGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAg                                                              <  1:3300280/69‑1 (MQ=255)
 aaTTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGc                                                             >  1:2606110/1‑69 (MQ=255)
 aaTTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGc                                                             >  1:2929786/1‑69 (MQ=255)
   tttCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCAc                                                           <  1:2904942/69‑1 (MQ=255)
   tttCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCAc                                                           <  1:242047/69‑1 (MQ=255)
   tttCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCAc                                                           <  1:2462302/69‑1 (MQ=255)
       aCCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGt                                                        <  1:2770705/68‑1 (MQ=255)
       aCCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGt                                                        <  1:1560295/68‑1 (MQ=255)
          aaTAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATc                                                    <  1:3008200/69‑1 (MQ=255)
          aaTAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATc                                                    <  1:2582921/69‑1 (MQ=255)
           aTAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATc                                                    <  1:3279490/68‑1 (MQ=255)
              aCTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTa                                                       >  1:2392970/1‑62 (MQ=255)
              aCTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTa                                                       >  1:2288677/1‑62 (MQ=255)
              aCTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTa                                                       >  1:1292873/1‑62 (MQ=255)
                ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:315232/67‑1 (MQ=255)
                ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:1837619/67‑1 (MQ=255)
                ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:80109/67‑1 (MQ=255)
                ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:622801/67‑1 (MQ=255)
                ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:1148902/67‑1 (MQ=255)
                ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:508264/67‑1 (MQ=255)
                ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:3314151/67‑1 (MQ=255)
                ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:429097/67‑1 (MQ=255)
                ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:1239962/67‑1 (MQ=255)
                ttttGCCGTGAGCAGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTCACAGCGCTGTAATCACTg                                                <  1:193745/67‑1 (MQ=255)
                ttttGCCGAGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:1413046/67‑1 (MQ=255)
                ctttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:2061566/66‑1 (MQ=255)
                        tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa                                      >  1:1410566/1‑67 (MQ=255)
                        tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa                                      >  1:2658558/1‑67 (MQ=255)
                        tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa                                      >  1:509537/1‑67 (MQ=255)
                        tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa                                      >  1:2829861/1‑67 (MQ=255)
                        tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa                                      >  1:3129163/1‑67 (MQ=255)
                        tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa                                      >  1:1540754/1‑67 (MQ=255)
                        tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa                                      >  1:1667411/1‑67 (MQ=255)
                                      gAAAAAGCCTCATCAACACGAAGTCTAATAGCTCTGTAATCACTGGGGAGAAAGAGGTTCATTTgcgc                         <  1:200157/68‑1 (MQ=255)
                                      gAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTgcgc                         <  1:581660/68‑1 (MQ=255)
                                                   caACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTcaca            >  1:326691/1‑68 (MQ=255)
                                                   caACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTcaca            >  1:979945/1‑68 (MQ=255)
                                                   caACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTcaca            >  1:135880/1‑68 (MQ=255)
                                                   caACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTcaca            >  1:1094098/1‑68 (MQ=255)
                                                   caACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTcaca            >  1:1141663/1‑68 (MQ=255)
                                                   caACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTcaca            >  1:224224/1‑68 (MQ=255)
                                                   caACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGATCATTTGCGCAATAGCATTcac             >  1:2908018/1‑67 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGgaga                                          >  1:1741398/1‑37 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGgaga                                          >  1:2781195/1‑37 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGgaga                                          >  1:3147351/1‑37 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGgaga                                          >  1:1399485/1‑37 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGgaga                                          >  1:2288255/1‑37 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGgaga                                          >  1:1179890/1‑37 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAgg                                    >  1:2883717/1‑43 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCCCATc          >  1:1026622/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          <  1:514258/69‑1 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:1175499/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:1088258/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:932475/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:3385205/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:960593/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:2322847/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:3308906/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          <  1:1815171/69‑1 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:1932477/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:3074455/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          <  1:2440204/69‑1 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:2808529/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          <  1:2321251/69‑1 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:2877775/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:2986692/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:3055528/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          <  1:1938604/69‑1 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          <  1:3074761/69‑1 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:3085988/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:311959/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:1990888/1‑69 (MQ=255)
                                                    aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          <  1:3226331/69‑1 (MQ=255)
                                                            gTCTAACAGCTCTGTAATCACTGGGGAGAAAAAGGTTCATTTGCGCAATAGCATTAACATCAACGCttt  >  1:1194068/1‑69 (MQ=255)
                                                            gTCTAACAGCGCTGTAATCACTGGGGAGAAAAAGGTTCATTTGCGCAATAGCATTAACATCAACGCttt  >  1:561007/1‑69 (MQ=255)
                                                            gTCTAACAGCGCTGTAATCACTGGGGAGAAAAAGGTTCATTTGCGCAATAGCATTAACATCAACGCttt  >  1:612716/1‑69 (MQ=255)
                                                            gTCTAACAGCGCTGTAATCACTGGGGAGAAAAAGGTTCATTTGCGCAATAGCATTAACATCAACGCttt  >  1:1431482/1‑69 (MQ=255)
                                                            gTCTAACAGCGCTGTAATCACTGGGGAGAAAAAGGTTCATTTGCGCAATAGCATTAACATCAACGCttt  >  1:2281485/1‑69 (MQ=255)
                                                            gTCTAACAGCGCTGTAATCACTGGGGAGAAAAAGGTTCATTTGCGCAATAGCATTAACATCAACGCttt  >  1:2058155/1‑69 (MQ=255)
                                                            gTCTAACAGCGCTGTAATCACTGGGGAGAAAAAGGTTCATTTGCGCAATAGCATTAACATCAACGCttt  >  1:3227566/1‑69 (MQ=255)
                                                                  |                                                              
CAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATCAACGCTTT  >  minE/1294587‑1294715

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: