Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,479,154 G→A 23.8% A113A (GCG→GCA phoU → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,479,1540GA23.8% 97.5 / 22.5 75A113A (GCG→GCAphoUDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (19/38);  new base A (6/12);  total (25/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

TAAAACCATTGCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTGGAGTCGCTGGGCCGTCAT  >  minE/2479094‑2479217
                                                            |                                                               
tAAAACCATTGCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGgaga                                                          <  1:244354/68‑1 (MQ=255)
  aaaCCATTGCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGaaa                                                        >  1:1039069/1‑68 (MQ=255)
  aaaCCATTGCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAAt                                                       >  1:1349356/1‑69 (MQ=255)
  aaaCCATTGCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAAt                                                       >  1:1541149/1‑69 (MQ=255)
  aaaCCATTGCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAAt                                                       >  1:2481336/1‑69 (MQ=255)
  aaaCCATTGCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAAt                                                       >  1:2709044/1‑69 (MQ=255)
  aaaCCATTGCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAAt                                                       >  1:1115380/1‑69 (MQ=255)
          gCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTccc                                                 >  1:2990508/1‑67 (MQ=255)
          gCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTccc                                                 >  1:2988807/1‑67 (MQ=255)
          gCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTccc                                                 >  1:428341/1‑67 (MQ=255)
          gCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTccc                                                 >  1:547567/1‑67 (MQ=255)
          gCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTccc                                                 >  1:1775180/1‑67 (MQ=255)
          gCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATACTccc                                                 >  1:3062528/1‑67 (MQ=255)
              aGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCagcagc                                           <  1:3169869/69‑1 (MQ=255)
              aGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCagcagc                                           <  1:2643911/69‑1 (MQ=255)
              aGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCagcagc                                           <  1:276698/69‑1 (MQ=255)
              aGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCagcagc                                           <  1:604429/69‑1 (MQ=255)
              aGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCagcagc                                           <  1:1989831/69‑1 (MQ=255)
              aGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCagcagc                                           <  1:1946244/69‑1 (MQ=255)
              aGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCagcagc                                           <  1:1742999/69‑1 (MQ=255)
              aGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCagcagc                                           <  1:1287992/69‑1 (MQ=255)
              aGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCagcagc                                           <  1:1346193/69‑1 (MQ=255)
              aGCTGGAGCGTATTGGCGACGTGGAGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCagcagc                                           <  1:77141/69‑1 (MQ=255)
               gCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCagcagc                                           <  1:2654266/68‑1 (MQ=255)
                     gCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCACTTGAGAAAtt                                                      <  1:3308700/51‑1 (MQ=255)
                     gCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCACTGGAGAAAtt                                                      <  1:2776886/51‑1 (MQ=255)
                     gCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCACTGGAGAAAtt                                                      <  1:335393/51‑1 (MQ=255)
                     gCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCACTGGAGAAAtt                                                      <  1:684239/51‑1 (MQ=255)
                     gCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCACTGGAGAAAtt                                                      <  1:3277864/51‑1 (MQ=255)
                     gCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCACTGGAGAAAtt                                                      <  1:305730/51‑1 (MQ=255)
                     gCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCACTGGAGAAAtt                                                      <  1:2923977/51‑1 (MQ=255)
                     gCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCACTGGAGAAAtt                                                      <  1:701667/51‑1 (MQ=255)
                     gCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCACTGGAGAAAtt                                                      <  1:1025019/51‑1 (MQ=255)
                     gCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCACTGGAGAAAtt                                                      <  1:1606906/51‑1 (MQ=255)
                     gCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCACTGGAGAAAtt                                                      <  1:1802003/51‑1 (MQ=255)
                     gCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCACTGGAGAAAtt                                                      <  1:2029959/51‑1 (MQ=255)
                       gTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGt                                  <  1:3350322/69‑1 (MQ=255)
                            ggCGACGTTGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTg                             <  1:2234054/69‑1 (MQ=255)
                            ggCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGc                               >  1:387198/1‑67 (MQ=255)
                            ggCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCt                              >  1:1450014/1‑68 (MQ=255)
                            ggCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCt                              >  1:770540/1‑68 (MQ=255)
                            ggCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCt                              >  1:2322971/1‑68 (MQ=255)
                            ggCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCt                              >  1:3300391/1‑68 (MQ=255)
                            ggCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCt                              >  1:1352993/1‑68 (MQ=255)
                            ggCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTg                             <  1:582551/69‑1 (MQ=255)
                            ggCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTg                             <  1:3436101/69‑1 (MQ=255)
                                 cGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAg                        <  1:56639/69‑1 (MQ=255)
                                 cGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAg                        <  1:1621491/69‑1 (MQ=255)
                                 cGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAg                        <  1:2111993/69‑1 (MQ=255)
                                 cGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAg                        <  1:2833099/69‑1 (MQ=255)
                                  gTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAg                        <  1:2884828/68‑1 (MQ=255)
                                  gTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAg                        <  1:2702495/68‑1 (MQ=255)
                                         aCAAAATCTGCCGTACTGCACTGGAGAAATTCTCCCAGCAGCATCa                                       >  1:1305993/1‑46 (MQ=255)
                                         aCAAAATCTGCCGTACTGCACTGGAGAAATTCTCCCAGCAGCATCAGCCATTGCTGGTAAGTCTGGAg                 >  1:1894821/1‑68 (MQ=255)
                                         aCAAAATCTGCCGTACTGCACTGGAGAAATTCTCCCAGCAGCATCAGCCATTGCTGGTAAGTCTGGAGt                >  1:2615109/1‑69 (MQ=255)
                                         aCAAAATCTGCCGTACTGCACTGGAGAAATTCTCCCAGCAGCATCAGCCATTGCTGGTAAGTCTGGAGt                >  1:1627013/1‑69 (MQ=255)
                                         aCAAAATCTGCCGTACTGCACTGGAGAAATTCTCCCAGCAGCATCAGCCATTGCTGGTAAGTCTGGAGt                >  1:1899125/1‑69 (MQ=255)
                                         aCAAAATCTGCCGTACTGCACTGGAGAAATTCTCCCAGCAGCATCAGCCATTGCTGGTAAGTCTGGAGt                >  1:679757/1‑69 (MQ=255)
                                          cAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTGGAGTc               <  1:3005602/69‑1 (MQ=255)
                                          cAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTGGAGTc               <  1:624877/69‑1 (MQ=255)
                                          cAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTGGAGTc               <  1:1414510/69‑1 (MQ=255)
                                          cAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTGGAGTc               <  1:1815126/69‑1 (MQ=255)
                                                cTGCCGTACTCCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTg                             <  1:2530073/49‑1 (MQ=255)
                                                     gTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTGGAGTCGCTGGGc        <  1:1169420/65‑1 (MQ=255)
                                                      tACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTTGTAAg                        <  1:2924973/48‑1 (MQ=255)
                                                        cTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTg                    <  1:861149/50‑1 (MQ=255)
                                                        cTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTg                    <  1:558626/50‑1 (MQ=255)
                                                        cTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTg                    <  1:1202794/50‑1 (MQ=255)
                                                        cTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTg                    <  1:1231437/50‑1 (MQ=255)
                                                        cTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTg                    <  1:3350696/50‑1 (MQ=255)
                                                        cTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTg                    <  1:1376668/50‑1 (MQ=255)
                                                        cTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTg                    <  1:2379838/50‑1 (MQ=255)
                                                        cTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTg                    <  1:2188515/50‑1 (MQ=255)
                                                        cTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTg                    <  1:1643648/50‑1 (MQ=255)
                                                        cTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTGGAGTCGCTGGGCCGTCAt  >  1:3421268/1‑68 (MQ=255)
                                                            |                                                               
TAAAACCATTGCCGAGCTGGAGCGTATTGGCGACGTGGCGGACAAAATCTGCCGTACTGCGCTGGAGAAATTCTCCCAGCAGCATCAGCCGTTGCTGGTAAGTCTGGAGTCGCTGGGCCGTCAT  >  minE/2479094‑2479217

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: