Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,103,998 A→C 31.7% intergenic (‑246/‑3) agaR ← / → kbaZ DNA‑binding transcriptional dual regulator/tagatose 6‑phosphate aldolase 1, kbaZ subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,103,9980AC31.7% 118.0 / 88.7 123intergenic (‑246/‑3)agaR/kbaZDNA‑binding transcriptional dual regulator/tagatose 6‑phosphate aldolase 1, kbaZ subunit
Reads supporting (aligned to +/- strand):  ref base A (28/56);  new base C (14/25);  total (42/81)
Fisher's exact test for biased strand distribution p-value = 8.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

TGACGCAGTATCAACTGAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAATTTA  >  minE/2103936‑2104059
                                                              |                                                             
tGACGCAGTATCAACTGAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGa                                                         <  1:3133266/69‑1 (MQ=255)
          tCAACTGAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCt                                                  >  1:1042804/1‑66 (MQ=255)
          tCAACTGAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCt                                                  >  1:1489849/1‑66 (MQ=255)
          tCAACTGAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCt                                                  >  1:1821138/1‑66 (MQ=255)
          tCAACTGAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCt                                                  >  1:2851656/1‑66 (MQ=255)
          tCAACTGAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCt                                                  >  1:3377138/1‑66 (MQ=255)
          tCAACTGAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCt                                                  >  1:109529/1‑66 (MQ=255)
          tCAACTGAAACAAAACGAAAGATTAAAATCGCAGTAATCTGAACTGGAGACGAAAGTGa                                                         >  1:2828610/1‑59 (MQ=255)
           cAACTGAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGAc                                               >  1:2666853/1‑68 (MQ=255)
           cAACTGAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGAc                                               >  1:1405270/1‑68 (MQ=255)
           cAACTGAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGAc                                               >  1:1626172/1‑68 (MQ=255)
            aaCTGAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACa                                              >  1:2120162/1‑68 (MQ=255)
                gAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAAt                                         <  1:1129018/69‑1 (MQ=255)
                gAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAAt                                         <  1:3082352/69‑1 (MQ=255)
                gAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAAt                                         <  1:1296346/69‑1 (MQ=255)
                gAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAAt                                         <  1:2113233/69‑1 (MQ=255)
                gAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAAt                                         <  1:828775/69‑1 (MQ=255)
                gAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTCACAGAAAt                                         <  1:405094/69‑1 (MQ=255)
                  aaCAAAACGAAAGATTAATATCGCAGTAATCGGAACTGGAGAGGCAAg                                                            >  1:2652903/1‑48 (MQ=255)
                  aaCAAAACGAAAGATTAATATCGCAGTAATCGGAACTGGAGAGGCAAGTGAAACATCTGACAGaaa                                          >  1:2550901/1‑66 (MQ=255)
                  aaCAAAACGAAAGATTAATATCGCAGTAATCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATg                                        >  1:1499408/1‑68 (MQ=255)
                  aaCAAAACGAAAGATTAATATCGCAGTAATCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATg                                        >  1:3318424/1‑68 (MQ=255)
                  aaCAAAACGAAAGATTAATATCGCAGTAATCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATg                                        >  1:364888/1‑68 (MQ=255)
                  aaCAAAACGAAAGATTAATATCGCAGTAATCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATg                                        >  1:274795/1‑68 (MQ=255)
                  aaCAAAACGAAAGATTAATATCGCAGTAATCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATg                                        >  1:2629987/1‑68 (MQ=255)
                  aaCAAAACGAAAGATTAATATCGCAGTAATCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATg                                        >  1:544669/1‑68 (MQ=255)
                  aaCAAAACGAAAGATTAATATCGCAGTAATCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATg                                        >  1:1925025/1‑68 (MQ=255)
                  aaCAAAACGAAAGATTAATATCGCAGTAATCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATg                                        >  1:81524/1‑68 (MQ=255)
                  aaCAAAACGAAAGATTAATATCGCAGTAATCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATg                                        >  1:863459/1‑68 (MQ=255)
                  aaCAAAACGAAAGATTAATATCGCAGTAATCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATg                                        >  1:1581709/1‑68 (MQ=255)
                  aaCAAAACGAAAGATTAATATCGCAGTAATCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATg                                        >  1:855663/1‑68 (MQ=255)
                   aCAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATg                                        <  1:1390433/67‑1 (MQ=255)
                   aCAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATg                                        <  1:3278962/67‑1 (MQ=255)
                   aCAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATg                                        <  1:834923/67‑1 (MQ=255)
                   aCAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATg                                        <  1:2250933/67‑1 (MQ=255)
                   aCAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGt                                      <  1:392353/69‑1 (MQ=255)
                   aCAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGt                                      <  1:3129498/69‑1 (MQ=255)
                   aCAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGt                                      <  1:678456/69‑1 (MQ=255)
                   aCAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGt                                      <  1:2432733/69‑1 (MQ=255)
                   aCAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGt                                      <  1:362874/69‑1 (MQ=255)
                   aCAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGt                                      <  1:682364/69‑1 (MQ=255)
                    cAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGt                                      <  1:2994852/68‑1 (MQ=255)
                     aaaaCGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTga                                    <  1:2562034/69‑1 (MQ=255)
                      aaaCGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTga                                    <  1:1356194/68‑1 (MQ=255)
                            aaGATTAATGTCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGa                                                <  1:280758/50‑1 (MQ=255)
                            aaGATTAATGTCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGa                                                <  1:683162/50‑1 (MQ=255)
                            aaGATTAATGTCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGa                                                <  1:362911/50‑1 (MQ=255)
                            aaGATTAATGTCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGa                                                <  1:1017040/50‑1 (MQ=255)
                            aaGATTAATGTCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGa                                                <  1:1704920/50‑1 (MQ=255)
                            aaGATTAATGTCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGa                                                <  1:1677422/50‑1 (MQ=255)
                            aaGATTAATGTCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGa                                                <  1:815648/50‑1 (MQ=255)
                            aaGATTAATGTCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGa                                                <  1:3341640/50‑1 (MQ=255)
                            aaGATTAATGTCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGa                                                <  1:2053313/50‑1 (MQ=255)
                            aaGATTAATGTCGCAGTAATCTGAACTGGACAGGAAAGTGAAACATCTGa                                                <  1:2259303/50‑1 (MQ=255)
                              gATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAGCATCTGACAGAAATGGTGAGACAGcaca                           <  1:2672048/69‑1 (MQ=255)
                              gATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGcaca                           <  1:2862496/69‑1 (MQ=255)
                                      tCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATg                                        <  1:2241919/48‑1 (MQ=255)
                                      tCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGt                                      <  1:441046/50‑1 (MQ=255)
                                              aTCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACaaa           >  1:2795075/1‑69 (MQ=255)
                                              aTCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACaaa           >  1:18997/1‑69 (MQ=255)
                                              aTCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACaaa           >  1:2029277/1‑69 (MQ=255)
                                              aTCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACaaa           >  1:132500/1‑69 (MQ=255)
                                              aTCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACaaa           >  1:2505921/1‑69 (MQ=255)
                                              aTCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACaaa           >  1:234459/1‑69 (MQ=255)
                                              aTCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACaaa           >  1:2059466/1‑69 (MQ=255)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:397192/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:342350/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:1253804/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:590052/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:880799/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:3310104/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:1527776/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:1776031/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:184008/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:1954829/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:2028029/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:1303072/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:2893400/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:2919334/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:3319798/44‑1 (MQ=38)
                                              aTCGGAACTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:1289147/44‑1 (MQ=38)
                                                  gAAGTGGAGATGCAAGTGAAACTTCTGTCAGAAATGGTga                                    <  1:347860/40‑1 (MQ=255)
                                                    aCTGGAGAGGCAAGTGAAACATCTGACAGAAATGGTga                                    <  1:1260856/38‑1 (MQ=255)
                                                     cTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAAtt    <  1:1324421/69‑1 (MQ=255)
                                                     cTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAAtt    <  1:1202077/69‑1 (MQ=255)
                                                     cTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAAtt    <  1:1428049/69‑1 (MQ=255)
                                                     cTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAAtt    <  1:55691/69‑1 (MQ=255)
                                                     cTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAAtt    <  1:740295/69‑1 (MQ=255)
                                                     cTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAAtt    <  1:341260/69‑1 (MQ=255)
                                                      tGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAAtt    <  1:155451/68‑1 (MQ=255)
                                                      tGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAAtt    <  1:1360059/68‑1 (MQ=255)
                                                      tGGAGAGGAAAGTGAAACATCGGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAAtt    <  1:378978/68‑1 (MQ=255)
                                                       tgagaGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAATTTa  >  1:2186585/2‑69 (MQ=255)
                                                       ggAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAATTTa  >  1:362242/1‑69 (MQ=255)
                                                       ggAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAATTTa  >  1:1775014/1‑69 (MQ=255)
                                                       ggAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAATTTa  >  1:1043362/1‑69 (MQ=255)
                                                       ggAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAATTTa  >  1:1653384/1‑69 (MQ=255)
                                                       ggAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAATTTa  >  1:2756367/1‑69 (MQ=255)
                                                       ggAGAGGAAAGAGAAACATCTGACAGAAATGGTGAGACAGCACaa                          >  1:291496/1‑45 (MQ=255)
                                                              aaaGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACa             <  1:47814/51‑1 (MQ=255)
                                                              |                                                             
TGACGCAGTATCAACTGAAACAAAACGAAAGATTAATATCGCAGTAATCTGAACTGGAGAGGAAAGTGAAACATCTGACAGAAATGGTGAGACAGCACAAAGCGGGCAAAACAAATGGAATTTA  >  minE/2103936‑2104059

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: