Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,150,016 G→A 12.5% intergenic (‑69/‑133) osmE ← / → nadE DNA‑binding transcriptional regulator/NAD synthetase, NH3/glutamine‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,150,0160GA12.5% 206.9 / 17.7 72intergenic (‑69/‑133)osmE/nadEDNA‑binding transcriptional regulator/NAD synthetase, NH3/glutamine‑dependent
Reads supporting (aligned to +/- strand):  ref base G (29/34);  new base A (4/5);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

ATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTCGC  >  minE/1149946‑1150083
                                                                      |                                                                   
atatTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTg                                                                     <  1:1650394/71‑1 (MQ=255)
 tatTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGa                                                                    <  1:4035270/71‑1 (MQ=255)
 tatTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGa                                                                    <  1:1958026/71‑1 (MQ=255)
 tatTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGa                                                                    <  1:1906565/71‑1 (MQ=255)
      cGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGa                                                                <  1:826042/70‑1 (MQ=255)
      cGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGa                                                                <  1:800124/70‑1 (MQ=255)
      cGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGa                                                                <  1:3056412/70‑1 (MQ=255)
      cGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTAAACGa                                                                <  1:768539/70‑1 (MQ=255)
      cGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTAAACGa                                                                <  1:1805650/70‑1 (MQ=255)
             ttGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTc                                                        >  1:3500672/1‑71 (MQ=255)
             ttGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTc                                                        >  1:1598951/1‑71 (MQ=255)
             ttGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTc                                                        >  1:2848332/1‑71 (MQ=255)
             ttGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTc                                                        >  1:4121582/1‑71 (MQ=255)
             ttGTTTATCAGCGTGTTAGATAAGACTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTc                                                        >  1:3452094/1‑71 (MQ=255)
               gTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCt                                                         <  1:2950792/68‑1 (MQ=255)
               gTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCt                                                         <  1:2230518/68‑1 (MQ=255)
               gTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTAAACGAAACGGCTcc                                                       >  1:2269093/1‑70 (MQ=255)
                   aTCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCtt                                                   <  1:1741849/70‑1 (MQ=255)
                   aTCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTAAACGAAACGGCTCCGCtt                                                   <  1:3322190/70‑1 (MQ=255)
                    tCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTc                                                 <  1:2105183/71‑1 (MQ=255)
                      aGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTAAACGAAACGGCTCCGCTTTCag                                               <  1:3321537/71‑1 (MQ=255)
                        cGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCagag                                             <  1:14268/71‑1 (MQ=255)
                        cGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCagag                                             <  1:1100035/71‑1 (MQ=255)
                         gtgtTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCagag                                             <  1:3155620/70‑1 (MQ=255)
                          tgtTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGa                                           >  1:3880746/1‑71 (MQ=255)
                          tgtTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGa                                           >  1:2751335/1‑71 (MQ=255)
                            ttAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGAtt                                         >  1:3388296/1‑71 (MQ=255)
                            ttAGATAAGCCTGGAATACATTGGGAGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGAt                                          <  1:588243/70‑1 (MQ=255)
                              aGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTcc                                       >  1:1683322/1‑71 (MQ=255)
                              aGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTAAACGAAACGGCTCCGCTTTCAGAGGATTcc                                       >  1:2033397/1‑71 (MQ=255)
                              aGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTAAACGAAACGGCTCCGCTTTCAGAGGATTcc                                       >  1:1931634/1‑71 (MQ=255)
                               gATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTcc                                       <  1:3646437/70‑1 (MQ=255)
                                 tAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGa                                                                <  1:436191/43‑1 (MQ=255)
                                 tAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCagag                                             >  1:1448937/1‑62 (MQ=255)
                                 tAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTg                                     <  1:1645037/70‑1 (MQ=255)
                                  aaGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGaaa                                                              >  1:3387741/1‑44 (MQ=255)
                                   aGCCTGGACTACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTa                                   <  1:761777/70‑1 (MQ=255)
                                   aGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGaaa                                                              >  1:1058952/1‑43 (MQ=255)
                                   aGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGaaa                                                              >  1:476593/1‑43 (MQ=255)
                                   aGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGaa                                                               >  1:1601678/1‑42 (MQ=255)
                                   aGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGt                                    >  1:2798072/1‑69 (MQ=255)
                                   aGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTAt                                  >  1:2802723/1‑71 (MQ=255)
                                   aGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTAt                                  >  1:2331211/1‑71 (MQ=255)
                                      cTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGAc                               >  1:1373791/1‑71 (MQ=255)
                                      cTGGAATACATTGGGCGCTTTTTCAAGCCCGTAAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGa                                <  1:4107948/70‑1 (MQ=255)
                                       tGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCaga                                              >  1:3652661/1‑55 (MQ=255)
                                          aaTACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGttt                           >  1:2648958/1‑71 (MQ=255)
                                            tACATTGGGCGCTTTTTCAAGCCCGTAAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGtttt                          >  1:3816861/1‑70 (MQ=255)
                                              cATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGAt                                          <  1:1134811/52‑1 (MQ=255)
                                              cATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTaa                        <  1:2371755/70‑1 (MQ=255)
                                                ttGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCt                                      <  1:1308857/54‑1 (MQ=255)
                                                ttGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTATATGACg                              >  1:2103543/1‑62 (MQ=255)
                                                      gcTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTc               <  1:2617879/71‑1 (MQ=255)
                                                      gcTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTc               <  1:1296087/71‑1 (MQ=255)
                                                       cTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCt                                      <  1:3459444/47‑1 (MQ=255)
                                                        tttttCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTg                                     <  1:3763264/47‑1 (MQ=255)
                                                        cttttCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAg             <  1:16824/70‑1 (MQ=255)
                                                         ttttCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTg                                     >  1:941649/1‑46 (MQ=255)
                                                         ttttCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGt                                    >  1:2851489/1‑47 (MQ=255)
                                                          tttCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCt                                      <  1:1230409/44‑1 (MQ=255)
                                                          tttCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGc            <  1:3626994/70‑1 (MQ=255)
                                                          tttCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGc            <  1:151137/70‑1 (MQ=255)
                                                            tCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGccc          <  1:3924459/70‑1 (MQ=255)
                                                            tCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCg         >  1:3042423/1‑71 (MQ=255)
                                                            tCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCg         >  1:1101593/1‑71 (MQ=255)
                                                              aaGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCAt                 >  1:2766189/1‑61 (MQ=255)
                                                               aGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTg      >  1:2293589/1‑71 (MQ=255)
                                                                gCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCAGCTGt     >  1:761604/1‑71 (MQ=255)
                                                                 cccGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTc    <  1:1310523/71‑1 (MQ=255)
                                                                  ccGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTc    <  1:1282833/70‑1 (MQ=255)
                                                                   cGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTCg   >  1:2747054/1‑70 (MQ=255)
                                                                   cGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTCGc  <  1:2079487/71‑1 (MQ=255)
                                                                     tGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCAtt                <  1:3657891/55‑1 (MQ=255)
                                                                      |                                                                   
ATATTCCGTCCTCTTGTTTATCAGCGTGTTAGATAAGCCTGGAATACATTGGGCGCTTTTTCAAGCCCGTGAACGAAACGGCTCCGCTTTCAGAGGATTCCTGTATGACGTTTTAACCACCATTCAGCCCGCTGTCGC  >  minE/1149946‑1150083

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: